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ABSTRACT 

We present measurements of the Hubble diagram for 103 Type la supcrnovae (SNe) ■with rcdshifts 
0.04 < z < 0.42, discovered during the first season (Fall 2005) of the Sloan Digital Sky Survey-II 
(SDSS-II) Supernova Survey. These data fill in the redshift "desert" between low- and high-redshift 
SN la surveys. Within the framework of the mlcs2k2 light-curve fitting method, we use the SDSS-II 
SN sample to infer the mean reddening parameter for host galaxies, Ry ~ 2.18 ± 0.14stat ± 0.48syst, 
and find that the intrinsic distribution of host-galaxy extinction is well fit by an exponential function, 
P{Av) = exp(— Ay/rv), with ry = 0.334 ± 0.088 mag. We combine the SDSS-II measurements 
with new distance estimates for published SN data from the ESSENCE survey, the Supernova Legacy 
Survey (SNLS), the Hubble Space Telescope (HST), and a compilation of nearby SN la measurements. 
A new feature in our analysis is the use of detailed Monte Carlo simulations of all surveys to account 
for selection biases, including those from spectroscopic targeting. Combining the SN Hubble diagram 
with measurements of baryon acoustic oscillations from the SDSS Luminous Red Galaxy sample and 
with cosmic microwave background temperature anisotropy measurements from WMAP, we estimate 
the cosmological parameters w and ^m, assuming a spatially flat cosmological model (FwCDM) with 
constant dark energy equation of state parameter, w. We also consider constraints upon J^m and Q/^ 
for a cosmological constant model (ACDM) with w = —1 and non-zero spatial curvature. For the 
FwCDM model and the combined sample of 288 SNe la, we find w = -0.76 ± 0.07(stat)±0.11(syst), 
= 0.307 ± 0.019(stat)±0.023(syst) using mlcs2k2 and w = -0.96 ± 0.06(stat) ± 0.12(syst), 
Om = 0.265 ± 0.016(stat) ± 0.025(syst) using the SALT~-ii fitter. We trace the discrepancy between 
these results to a difference in the rest-frame UV model combined with a different luminosity correction 
from color variations; these differences mostly affect the distance estimates for the SNLS and HST 
supcrnovae. We present detailed discussions of systematic errors for both light-curve methods and 
find that they both show data-model discrepancies in rest-frame [/-band. For the SALT-ii approach, 
we also see strong evidence for rcdshift-dcpendence of the color-luminosity parameter (/3) . Restricting 
the analysis to the 136 SNe la in the Near by -I- SDSS-II samples, wc find much better agreement 
between the two analysis methods but with larger uncertainties: w — —0.92 ± 0.13(stat) ^0 33 (syst) 
for mlcs2k2 and w = -0.92 ± O.ll(stat) ^015 (syst) for SALT-ii. 
Subject headings: supernova cosmology: cosmological parameters 
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1. INTRODUCTION 

Ten years ago, measurements of the Hubble diagram 
of Type la supernovae (SNe) provided the first direct 
evidence for cosmic acceleration (Riess et al. 1998; Perl- 
mutter et al. 1999). In the intervening decade, dedicated 
SN surveys have brought tremendous improvements in 
both the quantity and quality of SN la data, and SNe la 
remain the method of choice for precise relative distance 
determination over cosmological scales (e.g., Lcibundgut 
2001; Filippenko 2005). We now have in hand large, ho- 
mogeneously selected samples of SNe la with relatively 
dense time-sampling in multiple passbands at redshifts 
z > 0.3, most recently from the ESSENCE project (Mik- 
naitis et al. 2007; Wood-Vasey et al. 2007) and Supernova 
Legacy (SNLS) Survey (Astier et al. 2006), augmented by 
smaller samples from the Hubble Space Telescope (HST) 
that extend to higher redshift (Garnavich et al. 1998; 
Knop et al. 2003; Riess et al. 2004, 2007). These data 
have confirmed and sharpened the evidence for acceler- 
ated expansion. Cosmic acceleration is most commonly 
attributed to a new energy-density component known 
as dark energy (for a review, see Frieman, Turner, & 
Huterer (2008)). The recent SN measurements, in com- 
bination with measurements of the baryon acoustic os- 
cillation (BAO) feature in galaxy clustering and of the 
cosmic microwave background (CMB) anisotropy, have 
provided increasingly precise constraints on the density, 
r^DE, and equation of state parameter, w, of dark energy. 

Despite these advances, a number of concerns rc- 
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main about the robustness of current SN cosmology 
constraints. The SN la Hubble diagram is constructed 
from combining low- and high-rcdshift SN samples that 
have been observed with a variety of telescopes, instru- 
ments, and photometric passbands. Photometric off- 
sets between these samples are highly degenerate with 
changes in cosmological parameters, and these offsets 
could be hidden in part because there is a gap or "redshift 
desert" between the low-redshift (z < 0.1) SNe, found 
with small-aperture, wide-field telescopes, and the high- 
redshift {z > 0.3) SNe discovered by large- aperture tele- 
scopes with relatively narrow fields. In addition, the low- 
redshift SN measurements that are used both to anchor 
the Hubble diagram and to train SN distance estima- 
tors were themselves compiled from combinations of sev- 
eral surveys using different telescopes, instruments, and 
selection criteria. Increasing the robustness of the cos- 
mological results calls for larger supernova samples with 
continuous redshift coverage of the Hubble diagram; it 
also necessitates high-quality data, with homogeneously 
selected, densely sampled, multi-band SN light curves 
and well-understood photometric calibration. 

The Sloan Digital Sky Survey-II Supernova Survey 
(SDSS-II SN Survey) (Frieman et al. 2008), one of the 
three components of the SDSS-II project, was designed 
to address both the paucity of SN la data at interme- 
diate redshifts and the systematic limitations of previ- 
ous SN la samples, thereby leading to more robust con- 
straints upon the properties of the dark energy. Over the 
course of three three- month seasons, the SDSS-II SN Sur- 
vey discovered and measured well-sampled, multi-band 
light curves for roughly 500 spectroscopically confirmed 
SNe la in the redshift range 0.01 <z< 0.45. This data 
set fills in the redshift desert and for the first time in- 
cludes both low- and high-redshift SN measurements in 
a single survey. The survey takes advantage of the ex- 
tensive database of reference images, object catalogs, and 
photometric calibration previously obtained by the SDSS 
(for a description of the SDSS, see York et al. (2000)). 

In this paper, we present the Hubble diagram based 
on spectroscopically confirmed SNe la from the first full 
season (FaU 2005) of the SDSS-II SN Survey. To derive 
cosmological results, we include information from BAO 
(Eisenstein et al. 2005) and CMB measurements (Ko- 
matsu et al. 2008) , and we also combine our data with our 
own analysis of public SN la data sets at lower and higher 
redshifts. We fit the SN la light curves with two mod- 
els, mlcs2k2 (Jha et al. 2007) and SALT-ii (Guy et al. 
2007). We use the publicly available SALT-ii software 
with minor modifications, but we have made a number 
of improvements to the implementation of the mlcs2k2 
method, as described in §5. 

Two companion papers explore related analyses with 
the same SN data sets. Lampeitl et al. (2009) com- 
bine the SDSS-II SN data with different BAO constraints 
and with measurements of redshift-space distortions and 
of the Integrated Sachs- Wolfe effect to derive joint con- 
straints on dark energy from low-redshift (z < 0.4) mea- 
surements only; they also explore the consistency of the 
SN and BAO distance scales. SoUerman et al. (2009) use 
SN, BAO, and CMB measurements to constrain cosmo- 
logical models with a time- varying dark energy equation 
of state parameter as well as more exotic models for cos- 
mic acceleration. In all three papers, we use a consistent 
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analysis of the SN data. Differences in cosmological in- 
ferences are attributable to differences in (a) the SN data 
included, (b) the other cosmological data sets included, 
and (c) the cosmological model space considered. 

The outline of the paper is as follows. In §2, we 
briefly describe the operation and data processing for 
the SDSS-II SN Survey, which have been more exten- 
sively described in Sako et al. (2008). In §3, we sum- 
marize the spectroscopic analysis leading to final red- 
shift and SN type determinations (Zheng et al. 2008) 
and the photometric analysis leading to final supernova 
flux measurements (Holtzman et al. 2008) for SDSS-II 
SNe. In § 4, we present the SN samples and selection 
criteria applied to the light-curve data. We describe and 
compare the mlcs2k2 and SALT-ii methods in § 5. In 
§6 we describe detailed Monte Carlo simulations for the 
SDSS-II SN Survey and other SN data sets that we use 
to determine survey efficiencies and their dependences 
on SN luminosity, extinction, and redshift. Modeling of 
the survey efhciencies is needed to correct for selection 
biases that affect SN distance estimates. In §7 we use 
a larger spectroscopic -|-photometric SDSS-II SN sample 
to determine host-galaxy dust properties that are used in 
the mlcs2k2 fits. In particular, we present new measure- 
ments of the mean dust parameter, Ry, and of the extinc- 
tion (Av) distribution. In §8 we describe the cosmologi- 
cal likelihood analysis, which combines the SN la Hubble 
diagram with BAO and CMB measurements. In §9, we 
present a detailed study of systematic errors, showing 
how uncertainties in model parameters and in calibra- 
tions impact the results. In §10 we discuss the super- 
nova Hubble diagrams using the mlcs2k2 and SALT-ii 
fitters and derive constraints on cosmological parame- 
ters. We provide a detailed comparison of the MLCS2k2 
and SALT-II results in § 11, and we conclude in § 12. 
Appendices provide details on the methods for warping 
the SN la spectral template for K-corrcctions, model- 
ing the filter passbands for the nearby SN la sample, 
determining the magnitudes of the primary photometric 
standard stars, extracting the distribution of host-galaxy 
dust extinction from the SDSS-II sample, and estimat- 
ing error contours that include systematic uncertainties. 
They also include discussion of the scatter in the SDSS- 
II Hubble diagram and of the translation of the SALT-ii 
model into the mlcs2k2 framework. 

2. SDSS-II SUPERNOVA SURVEY 

The scientific goals, operation, and basic data process- 
ing for the SDSS-II SN Survey are described in Frieman 
et al. (2008), and details of the SN search algorithms 
and spectroscopic observations are given in Sako et al. 
(2008). Here we provide a brief summary of the Fall 
2005 campaign, in order to set the context for the data 
analysis. 

The SDSS-II Supernova Survey primary instrument 
was the SDSS CCD camera (Gunn et al. 1998) mounted 
on a dedicated 2.5-m telescope (Gunn et al. 2006) at 
Apache Point Observatory (APO), New Mexico. The 
camera obtains, nearly simultaneously, images in five 
broad optical bands: ugriz (Fukugita et al. 1996). The 
camera was used in time-delay-and-intcgrate (TDI, or 
drift scan) mode, which provides efficient sky coverage. 
The Supernova Survey scanned at the normal (sidereal) 
SDSS survey rate, which yielded 55-s integrated expo- 



sures in each passband; the instrument covered the sky 
at a rate of approximately 20 square degrees per hour. 

On most of the usable observing nights in the period 1 
September through 30 November 2005, the SDSS-II SN 
Survey scanned a region (designated stripe 82) centered 
on the celestial equator in the Southern Galactic hemi- 
sphere that is 2.5° wide and runs between right ascen- 
sions of 20'" and 4'^'', covering a total area of 300 sq. deg. 
Due to gaps between the CCD columns, on a given night 
slightly more than half of the declination range of the 
stripe was imaged; on succeeding nights, the survey al- 
ternated between the northern (N) and southern (S) dec- 
lination strips of stripe 82 (see Stoughton et al. (2002) 
for a description of the SDSS observing geometry). Ac- 
counting for CCD gaps, bad weather, nearly full Moon, 
and other observing programs, a given region was im- 
aged on average every four to five nights under a variety 
of conditions. This relatively high cadence enabled us to 
obtain well-sampled light curves, typically starting well 
before peak light. 

At the end of each night of imaging, the SN data 
were processed using a dedicated 20-CPU computing 
cluster at APO. Images were processed through the 
PHOTO photometric reduction pipeline to produce cor- 
rected w, g, r, z, z frames (Lupton et al. 2001; Ivezic et al. 
2004), each with an astrometric solution (Pier et al. 
2003), point-spread- function map, and zero-point. A co- 
added template image, consisting of typically 8 stacked 
images taken in previous years, was matched to the 
new image and subtracted from it. Subtracted gri im- 
ages were searched for pixel clusters with an excess flux 
(roughly 3a) above the noise in the subtracted image, 
and a position and total PSF flux were assigned for each 
significant detection. We positionally matched detec- 
tions in multiple passbands: objects are detections in 
at least two of the three gri passbands with a displace- 
ment of less than 0.8" between detections in each filter. 
This displacement cut was chosen to ensure high effi- 
ciency for objects with low signal-to-noise. The g and 
r exposures of a given object were taken five minutes 
apart, enabling many fast asteroids to be removed by the 
0.8" requirement. Finally, a catalog of 10^ previously de- 
tected variables (mainly stars and AGN) and 4 million 
stars (r < 21.5) was used to reject detections within 1" of 
any object in the catalog; nearly 40,000 such detections 
were automatically vetoed during the Fall 2005 survey. 

During the season, 20" x 20" cut-outs of the resulting 
~ 140, 000 object images were visually scanned by hu- 
mans, '^^ typically within 24 hours of when the data were 
obtained. The human scanning was done to eliminate 
objects that were clearly not supernovae, such as un- 
subtracted diffraction spikes, other subtraction artifacts, 
and obvious asteroids. To monitor the software pipelines 
and human scanning efficiency, 'fake' supernovae were 
inserted on top of galaxies in the images during process- 
ing. Approximately 11,400 of the objects were tagged by 
a scanner as a possible supernova candidate. Nearly 60% 
of the candidates appeared only once during the survey; 
most of these are likely slow-moving solar system ob- 
jects. After a night of observations, each candidate light 

During the 2006 season we implemented more aggressive soft- 
ware cuts that reduced the number of objects scanned by over an 
order of magnitude. 
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curve (in g, r, i) was updated and compared with a set 
of supernova light curve templates that include SNe la 
as a function of redshift, intrinsic luminosity, and ex- 
tinction, as well as non-la SN types. Light curves that 
matched best to a SN la template (at any reasonable 
redshift, luminosity, and extinction) were preferentially 
scheduled for spectroscopic follow-up observations. Can- 
didates with r-band magnitude r < 20.5 were given high- 
est priority for follow-up, regardless of photometric SN 
type; for SNe la, this magnitude cut corresponds roughly 
to redshifts z < 0.15. For fainter SN la candidates, spec- 
troscopic priority was given to candidates with the best 
chance of acquiring a useful spectrum. In order of im- 
portance, the prioritization criteria were: (i) SN is well- 
separated (> 1") from the core of its host galaxy, (ii) 
reasonable SN/galaxy brightness contrast based on vi- 
sual inspection, and (iii) SN host-galaxy is relatively red 
(early- type). In most detection in at least two 

epochs was required before a spectrum was obtained. 

Spectra of supernova candidates and, where possible, 
their host galaxies were obtained in Sept. -Dec. 2005 with 
a number of telescopes (Frieman et al. 2008; Zheng et al. 
2008): the Hobby-Eberly 9.2-meter at McDonald Ob- 
servatory, the Astrophysical Research Consortium 3.5- 
meter at APO, the Subaru 8.2-meter on Mauna Kea, the 
Hiltner 2.4-mcter at MDM Observatory, the 4.2-meter 
William Herschel Telescope on La Palma, and the Keck 
10-meter on Mauna Kea. Approximately 90% of the 
SN la candidates that were spectroscopically observed 
were confirmed as SNe la. 

As noted below (§3.1), 146 spectroscopically observed 
candidates from 2005 were classified as definitive or pos- 
sible SNe la based on analysis of their spectra. For these 
candidates, there are a total of more than 2000 photomet- 
ric epochs, where each epoch corresponds to a measure- 
ment (not necessarily a detection) in the ugriz passbands 
within a time window of —15 days to -1-60 days relative to 
peak brightness in the supernova rest-frame. About half 
of the epochs were recorded in "photometric" conditions, 
defined as no moon, PSF less than 1.7", and no clouds 
as indicated by the SDSS cloud camera, which moni- 
tors the sky at lO/^m (Hogg et al. 2001). Another 30% 
of the measurements were recorded in non-photometric 
(but moonless) conditions. The remaining 20% of the 
measurements were taken with the moon above the hori- 
zon. 

3. SDSS SN SPECTROSCOPIC AND PHOTOMETRIC 
REDUCTION 

For each supernova candidate found during the sur- 
vey, the on-mountain software pipeline described in § 2 
delivered preliminary photometric measurements. Sim- 
ilarly, spectroscopic observations were reduced in near- 
real time so that estimates of SN type and redshift could 
be made. Although these initial measurements were suf- 
ficient for discovering and confirming SNe, for the final 
analysis and sample selection we require more accurate 
photometry (Holtzman et al. 2008) and a more uniform 
spectroscopic analysis (Zheng et al. 2008). This section 
briefly describes these techniques. 

3.1. Supernova Typing and Redshift Determination 

After the finish of the Fall 2005 season, all of the su- 
pernova spectra were processed with IRAF (Tody 1993). 



Classification of the reduced SN spectra was aided by 
the IRAF package rvsao.xcsao (Tonry & Davis 1979), 
which cross-correlates the spectra with libraries of SN 
spectral templates and searches for significant peaks. De- 
tails of this analysis are described in Zheng et al. (2008). 
About half of the supernova spectra had an excellent 
template match, while the other half required more hu- 
man judgment for the SN typing. Based on this analysis, 
130 candidates were classified as confirmed SNe la and 
16 candidates were classified as probable SNe la. 

For 29 of these 146 candidates, we have used the SDSS 
host-galaxy spectroscopic redshift as reported in the 
SDSS DR4 database; typical redshift uncertainties are 
l-2xl0~^. For SN 2005hj, a host-galaxy spectroscopic 
redshift and its uncertainty were obtained by Quimby 
et al. (2007). For 82 of the candidates that do not have a 
host spectroscopic redshift in the DR4 database, we use 
the redshift from host-galaxy spectral features obtained 
with our own spectroscopic observations. The redshift 
precision in those cases is estimated to be 0.0005, the rms 
difference between our host-galaxy redshifts and those 
measured by the SDSS spectroscopic survey (DR4) for 
a sample in which both redshifts are available. For the 
remaining 34 candidates, our redshift estimate is based 
on spectroscopic features of the supernovae, with an esti- 
mated uncertainty of 0.005, the rms spread between the 
SN redshifts and host-galaxy redshifts. In summary, 77% 
of the spectroscopically confirmed and probable SNe la 
have spectroscopic redshifts determined from host-galaxy 
features, while the rest have redshifts based on SN spec- 
tral features. The redshifts are determined in the helio- 
centric frame and then transformed to the CMB frame 
as described in § 8. 

The redshift distribution for the 130 confirmed SNe la 
from the 2005 season is shown below in Fig. 2e. The rel- 
ative deficit of confirmed SNe at redshifts between 0.15 
and 0.25 is due to the finite spectroscopic resources that 
were available for the Fall 2005 campaign and to the 
relative priorities given to low- and high-redshift candi- 
dates for the different telescopes (Sako et al. 2008). Sub- 
sequently, host-galaxy redshifts have been obtained for 
most of the "missing" SN la candidates with SN la-like 
light curves in this redshift range. These photometri- 
cally identified (but spectroscopically unconfirmed) can- 
didates with host-galaxy redshifts are used in the deter- 
mination of host-galaxy dust properties (§ 7), but we do 
not include them in the Hubble diagram for this analy- 
sis. Compared to the Fall 2005 season, spectroscopic ob- 
servations during the 2006 and 2007 seasons were more 
complete around redshifts z ^ 0.2. 

3.2. Supernova Photometry 

To achieve precise and reliable SN photometry, we de- 
veloped a new technique called "Scene Model Photome- 
try" (SMP) that optimizes the determination of super- 
nova and host-galaxy fluxes. This method and the Fall 
2005 SN photometry results are described in detail in 
Holtzman et al. (2008). 

The basic approach of SMP is to simultaneously model 
the ensemble of survey images covering a SN location as 
a time- varying point source (the SN) and sky background 
plus time-independent galaxy background and nearby 
calibration stars, all convolved with a time- varying PSF. 
The calibration stars arc taken from the SDSS catalog 
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for stripe 82 produced by Ivezic et al. (2007). The fit- 
ted parameters are supernova position, supernova flux 
for each epoch and passband, and the host galaxy inten- 
sity distribution in each passband. The galaxy model for 
each passband is a 20 x 20 grid (with a grid-scale set by 
the CCD pixel scale, 0.4" x 0.4") in sky coordinates, and 
each of the 400 x 5 = 2000 galaxy intensities is an inde- 
pendent fit parameter. As there is no pixel re-sampling 
or image convolution, the procedure yields correct sta- 
tistical error estimates. Holtzman et al. (2008) describes 
the rigorous tests that were carried out to validate the 
accuracy of SMP photometry and of the error estimates. 

Although we have obtained additional imaging on 
other telescopes for a subsample of the confirmed SNe la, 
only photometry from the SDSS 2.5 m telescope is used 
in this analysis. Fig. 1 shows four representative SDSS-II 
SN la light curves processed through SMP and provides 
an indication of the typical sampling cadence and signal- 
to-noise as a function of redshift. 

The fluxes and magnitudes returned by SMP are in the 
native SDSS system (Ivezic et al. 2007). The SDSS pho- 
tometric system is nominally on the AB system, but the 
native flux in each filter differs from that of a true AB 
system by a small amount. AB-magnitudes are obtained 
by adding the AB-offsets in Table 1 to the native magni- 
tudes. The offsets are determined by comparing photo- 
metric measurements of the HST standard solar analogs 
P3330E, P177D, and P041C with synthetic magnitudes 
based on the published HST spectra (Bohlin 2006) and 
SDSS filter bandpasses. Since the standard stars are too 
bright to be measured directly with the SDSS 2.5 m tele- 
scope, the measurements are taken with the 0.5-meter 
SDSS Photometric Telescope (the PT) and transformed 
to the native system of the SDSS telescope. The tech- 
nique of transferring the PT magnitudes to the native 
SDSS system is identical to that used to obtain the SDSS 
photometric calibration (Tucker et al. 2006). The uncer- 
tainty in the AB offsets is estimated to be 0.003, 0.004, 
0.004, 0.007, 0.010 mag (for u,g,r,i,z) based on the in- 
ternal consistency of the three standard solar analogs. 
The uncertainties given in Table 1 are larger, since they 
also account for the ~ 10 A uncertainties in the cen- 
tral wavelengths (given in the same Table) of the SDSS 
filters. 



TABLE 1 

AB OFFSETS AND CENTRAL WAVELENGTH UNCERTAINTIES FOR THE 
SDSS FILTERS. 



SDSS filter 


AB offset (mag) and 
its uncertainty^ 


uncertainty (A) on 
central wavelength 


u 


-0.037 ± 0.014 


8 


9 


-f 0.024 ± 0.009 


7 


r 


+0.005 ± 0.009 


16 


i 


+0.018 ±0.009 


25 


z 


+0.016 + 0.010 


38 



Errors account for uncertainties in the central wavelengths of 
the SDSS filters. 



4. SUPERNOVA SAMPLE SELECTION 

In this section, we describe the light-curve selection 
criteria used to define the SN la samples. To mini- 



mize systematic errors associated with analysis meth- 
ods and assumptions, we perform a nearly uniform anal- 
ysis on data from SDSS-II, the published data from 
ESSENCE (Wood-Vasey et al. (2007); hereafter WV07), 
SNLS (Astier et al. 2006), HST (Riess et al. 2007), and a 
nearby SN la sample collected over a decade from several 
surveys and a number of telescopes (Jha, Riess, & Kirsh- 
ner (2007) ; hereafter JRK07) . Although these data sam- 
ples are analyzed in a homogeneous fashion, we present 
more details about the SDSS-II analysis since these data 
are presented here for the first time and, more impor- 
tantly, because we use the SDSS-II sample in § 7 to make 
inferences about the SN la population that we apply to 
all the data samples. 

Light curves with good time sampling and good signal- 
to-noise are needed to yield reliable distance estimates. 
We therefore apply stringent selection cuts to all five pho- 
tometric data samples used in this analysis. The cuts are 
also chosen to define samples whose selection functions 
can be reliably modeled with the Monte Carlo simula- 
tions described in §6. In future analyses the cuts will be 
further refined based on studies with simulated samples. 

We first present the selection cuts we have applied and 
then discuss briefly the rationale for each of them. Defln- 
ing Tiost as the rest-frame time, such that Trost = corre- 
sponds to peak brightness in rest-frame _B-band accord- 
ing to mlcs2k2, we select for inclusion in the cosmology 
analysis SN la light curves that satisfy the following cri- 
teria: 

1. For SDSS-II, ESSENCE, SNLS, and HST, at least 
one measurement is required before peak brightness 
(Ticst < days); for the nearby sample, at least 
one measurement is required with Tj-est < +5 days. 
The requirement on the nearby sample is relaxed, 
because nearly half the sample would be rejected 
by the more stringent cut of Trcst < days. 

2. at least one measurement with Tiost > +10 days. 

3. at least five measurements with — 15 < Tj-ost < +60 
days. 

4. at least one measurement with signal-to-noise ra- 
tio (SNR) above 5 for: each of SDSS g, r, and 
i\ both SNLS r and i (no requirement on z); 
HST F814W.WFPC2 and at least one other HST 
passband. For the ESSENCE sample, we adopt 
the cuts from WV07: at least one measurement at 
Trost < +4 days that has SNR> 5, at least one 
measurement at Trcst > +9 days that has SNR> 5, 
and at least 8 total measurements with SNR> 5. 
Since the nearby SN la sample includes only events 
with high SNR, no SNR requirement is needed for 
that sample. 

5. T'fit > 0.001, where Pflt is the mlcs2k2 hght-curve 
fit probability based on the per degree of free- 
dom (see §5.1). 

6. z > Zmin = 0.02, which only affects the nearby 
SN la sample. 

For all the data samples we only include unambiguous 
spectroscopically confirmed SNe la; in particular, for the 
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Fig. 1.— Light curves for four SDSS-II SNc la at different redshifts: SN 2005fr at z = 0.09, SN 2005fb at z = 0.18, SN 2005fr at 
z = 0.29, and SN 2005gq at z = 0.39. The passbands are SDSS g (top), r (middle), and i (bottom). Points are the SMP flux measurements 
(flux = 10(11-0"'"), where m is the SN magnitude) with ±1 cr photometric errors indicated. Solid curves show the best- fit mlcs2k2 model 
fits (see §5.1), and dashed curves give the ±1 cr error bands on the model fits. The Modified Julian Date (MJD) under each set of light 
curves is the fitted time of peak brightness for rest-frame B-band. 



SDSS-II sample, we do not include the 16 spectroscopi- 
cally probable SNe la (see §3.1). Moreover, for the SDSS- 
II sample, we use only g,r,i photometry in the analysis, 
and we reject ~ 4% of the epochs for which the Scene 
Model Photometry pipeline (§3.2) did not return a reli- 
able flux estimate. 

For the first three requirements in the list above, a 
"measurement" corresponds to a recorded photometric 
measurement in a single passband and can have any 
signal-to-noise value, i.e., a significant detection is not 
necessary. These requirements collectively ensure that 
the time sampling of the light curve is sufficient to yield 
a robust light-curve model fit, with coverage before and 
after peak light so that the epoch of peak light can be reli- 
ably estimated. To illustrate the motivation for requiring 
a measurement before peak light (which was not explic- 
itly required in either the Astier et al. (2006) or WV07 
analyses), consider SN gl33 in the ESSENCE sample, 
which has no measurements before peak and is there- 
fore rejected in our analysis. Compared to the published 
values in WV07, our mlcs2k2 fitted time of maximum 
brightness is 10 days earlier, and our fitted distance mod- 
ulus is 0.4 magnitudes smaller. The fourth requirement, 
on SNR, similarly puts a fioor on the quality of the light- 
curve data. The fifth requirement, on mlcs2k2 light- 
curve fit probability Vfn , is designed to remove obviously 
peculiar SNe in an objective fashion. This cut removes 
the previously identified peculiar SNe la in the SDSS- 
II sample, 2005hk (Phillips et al. 2007; Chornock et al. 
2006), 2005gj (Aldering et al. 2006; Prieto et al. 2007), 
and SDSS-II SN 7017 (which is similar to 2005gj). In 
the nearby sample, it rejects the following peculiar SNe: 
1992bg, 1995bd, 1998de, 1999aa, 1999gd, 2001ay, 2001bt, 
2002bf, and 2002cx. 

The sixth selection criterion, corresponding to czmin = 
6000 km/sec, removes objects from the nearby sample 



for which the typical galaxy peculiar velocity, Vpoc ''^ 300 
km/scc, is a non-negligible fraction of the Hubble reces- 
sion velocity. In principle, this cut on redshift could be 
replaced by a redshift- and position-dependent weight- 
ing covariance factor that includes the effects of both 
random and correlated peculiar velocities (Hui & Greene 
2006; Cooray & Caldwell 2006). In this analysis, we fol- 
low recent practice and simply impose a lower redshift 
bound, but this approach raises the issue of how to select 
z,„i„. Astier et al. (2006) and WV07 used z^in = 0.015. 
However, using MLCS2k2, JRK07 found that the Hubble 
parameter inferred from the lowest-redshift SNe, with 
z < 0.025, is systematically higher than that obtained 
using more distant (0.025 < z < 0.1) nearby SNc, consis- 
tent with an earlier result of Zehavi et al. (1998). JRK07 
also noted that varying z^in from 0.008 to 0.027 changes 
the dark energy equation of state by Sw ^ 0.2 for the 
nearby SN la sample in combination with a simulated 
ESSENCE sample. As a consequence, Riess et al. (2004, 
2007) used z^in = 0.023 (czmin = 7000 km/sec), i.e., they 
only included SNe beyond the so-called "Hubble bubble" . 
On the other hand, Conley et al. (2007) found that the 
Hubble bubble is not significant when the SALT-ii fitter 
is used. As discussed in §9.1, we find that the best-fit 
value of w is sensitive to the choice of Zmin whether we 
use mlcs2k2 or SALT-ii. Varying Zmin, we find that 
•Zmin = 0.02 corresponds to the middle of the range of 
w variations for the MLCS2k2 method. For the SALT-ii 
method, w varies rapidly with ZaUn near Zmin ~ 0.015, 
and is more stable when Zmin ^ 0.02. On this basis, we 
choose Zniin = 0.02 for both light-curve fitting methods 
and include the effects of varying Zj^in in the systematic 
error budget. 

For the SDSS-II sample of 130 spectroscopically con- 
firmed SNe la from the Fall 2005 season, 103 satisfy 
these selection criteria. The cut-rejection statistics are 
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as follows: 3 are photometrically peculiar SNc la that 
fail the Vat requirement; 9 have no measurement be- 
fore peak brightness — most of these were discovered early 
in the survey season; 11 have no measurement with 
Trest > +10 days — most of these were discovered late 
in the survey season or were at the high-redshift end of 
the distribution; and 4 SNe la in the high-redshift tail, 
z ~ 0.4, fail the SNR requirement. 

With the selection criteria defined above, the number 
of SN la events used for fitting is shown in Table 2 for 
each sample; a total of 288 SNe la are included in the 
fiducial analysis (in systematic error tests, e.g., varying 
Zmin, this number fluctuates by a small amount). Ta- 
ble 2 also shows the average number of measurements 
per SN la for each sample, where a measurement is an 
observation in a single passband in the rest-frame time 
interval —15 to -1-60 days. The average number of mea- 
surements is about 50 for both the nearby and SDSS-II 
samples, in the twenties for ESSENCE and SNLS, and 
11 for HST. We note that our selection requirements 
are more restrictive than those applied in previous anal- 
yses. WV07 included 60 out of 105 spectroscopically 
confirmed ESSENCE SNe la for their MLCS2k2 analy- 
sis'^, while our cuts select 56. WV07 selected 45 SNe la 
from the nearby sample {z > 0.015), while we include 
33 {z > 0.02); the difference is mainly due to the differ- 
ent redshift cuts. Astier et al. (2006) included 71 SNLS 
SNe la, while we retain 62 from the same sample. 



TABLE 2 

Redshift range, number of SNe passing selection cuts, and 
mean number of measurements for each sn sample. 



training the models using pre-existing data, and different 



sample 


redshift 






(obs passbands) 


range 






Nearby (UBVRI) 


0.02 - 0.10 


33 


52 


SDSS-II {gri) 


0.04 - 0.42 


103 


48 


ESSENCE {RI) 


0.16 - 0.69 


56 


21 


SNLS {griz) 


0.25 - 1.01 


62 


27 


HST (FllOW, F160W, 


0.21 - 1.55 


34 


11 


F606W, F775W, F850LP) 









''Number of SNe la passing cuts. 

''Average number of measurements per SN la, in the interval 
-15 < Trcst < -1-60 days. 



Figure 2 shows distributions in the SDSS-II sample — 
before selection cuts — for some of the variables used in 
sample selection, as well as the SDSS-II SN la redshift 
distribution before and after selection cuts are applied. 
Figure 3 shows the redshift distribution for all five sam- 
ples, along with the average of the maximum observed 
signal-to-noise as a function of redshift. 

5. LIGHT-CURVE ANALYSIS 

In this section, we describe our methods of analyzing 
supernova light curves and extracting distance estimates. 
The two light-curve fitting methods we employ, mlcs2k2 
and SALT-ii, reflect different assumptions about the na- 
ture of color variations in SNe la, different approaches to 

39 WV07 include 60 ESSENCE SNe la for the analysis that 
includes the SNLS sample; for their analysis of Nearby-|-ESSENCE 
(excluding SNLS) they require z < 0.67, resulting in 57 ESSENCE 
SNe la. 
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Fig. 2. — For the spectroscopically confirmed SN la sample from 
the SDSS-II 2005 season, distributions are shown for: (a) number 
of gri measurements with — 15 < Trcst < +60 days as a function 
of day in the survey season when the SN reached peak luminosity. 
Vertical arrows show the start (Sept. 1) and end dates (Nov. 30) 
of the survey season. SNe that lie in the overlap region of strips 
82N and 82S (solid dots) tend to have more measurements; (b) 
logjQ of maximum g-band SNR versus redshift; (c) time of first 
measurement relative to peak light in rest-frame B; (d) time of 
last measurement (not necessarily detection) relative to peak light 
- the pile-up near 60 days is from SNe that have measurements 
past 60 days; (e) redshifts before (130, thin line) and after (103, 
thick line) selection cuts are applied. The arrows in panels (b), (c), 
and (d) indicate the selection cuts. 



nearby (JRK07) 
SDSS-II (this paper) 
ESSENCE (WV07) 
SNLS (Astier06) 
HST (Riess07) 





40 


Ln 


35 


o 




o 


30 


CD 


25 


Q. 


20 




CO 


0) 


15 


'i_ 






10 


LU 






^ . >l I II 



O nearby (JRK07) 

• SDSS-II (this paper) 

■ ESSENCE (WV07) 

□ SNLS (Astier06) 

^.0 HST(Riess07) 



O 



b 



><ioo 

£ 80 
60 
40 
20 


0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

redshift 

Fig. 3. — Top panel: summed redshift distribution for the five 
SN la samples indicated in the legend. Bottom panel: maximum 
observed signal-to-noise ratio (among all passbands) as a function 
of redshift, averaged in bins of width Az = 0.05; error bars indicate 
the rms spread within each bin. All selection requirements have 
been applied. 
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ways of determining model parameters. 

5.1. MLCS2k2 Fitting Method 

The Multicolor Light Curve Shape method, known as 
MLCS2k2 in its current incarnation (JRK07), has been 
in use for more than a decade; the original MLCS ver- 
sion (Riess et al. 1998) was used by the High-z Supernova 
Team in the discovery of cosmic acceleration. For each 
supernova, MLCS2k2 returns an estimated distance mod- 
ulus and its uncertainty; the redshift and distance modu- 
lus for each SN are inputs to the cosmology fit discussed 
in§ 8. 

MLCS2k2 describes the variation among SN la light 
curves with a single parameter (A). Excess color vari- 
ations relative to the one-parameter model are assumed 
to be the result of extinction by dust in the host galaxy 
and in the Milky way. The mlcs2k2 model magnitude 
is given by 

+ K£ + K}f,+f^ + X^i, , (1) 

where e is an epoch index that runs over the observa- 
tions, / are observer-frame filter indices. /' = UBVRI 
are the rest-frame filters for which the model is de- 
fined, A is the mlcs2k2 shape-luminosity parameter 
that accounts for the correlation between peak luminos- 
ity and the shape/duration of the light curve, Xhost is the 
host-galaxy extinction, Xmw is the Milky Way extinc- 
tion, Kffi is the K~correction between rest-frame and 
observer-frame filters, and /i is the distance modulus, 
which satisfies ^ = Slog^oCi^L/lO pc), where is the 
luminosity distance. We use this model for SN epochs in 
the rest- frame time range — 15 < Tjcst < +60 days rel- 
ative to rest-frame B-band maximum. Observer-frame 
passbands are included that satisfy 3200 < A//(l -I- z) < 
9500 A, where A/ is the mean wavelength of the filter 
passband, and z is the redshift of the SN la. To ac- 
count for larger model uncertainties in the restframe-UV 
region, a K-correction uncertainty of 0.0006 x (3500 — 
A/) mag is added in quadrature to the model error for 
A/ < 3500 A. 

In the mlcs2k2 model, the shape-luminosity param- 
eter A describes the intrinsic SN color dependence on 
brightness, and Xhost describes SN color variations from 
reddening (extinction) by dust in the host galaxy, which 
is assumed to behave in a manner similar to dust in 
the Milky Way. In particular, the extinction is de- 
scribed by the parametrization of Cardelli, Clayton, & 

Mathis (1989) (hereafter CCM89), = C'^' [a^' + 

/Rv)Av^ where Ay is the extinction in magnitudes 
in l/-band, = 1, &^ = 0, and the relative extinc- 
tion in other passbands is determined by the parame- 
ter the ratio of F-band extinction to color-excess, 
Rv = Av/E{B - V). For the Milky Way, the value of 
Rv averaged over a number of lines of sight is Ry = 3.1; 
this global value has been adopted in previous SN anal- 
yses using mlcs2k2. For the galaxies that host SNe la, 
we instead adopt Ry — 2.18 ± 0.50, as derived in § 7.2 
from the SDSS-II SN data. 

The coefficients M'^'-^ , p'^'^ , and q^'^ are model vec- 
tors that have been evaluated using nearly 100 well- 



observed low-redshift SNe as a training set. M'^'f is 
the absolute magnitude for a SN la with A = 0. As- 
suming a Hubble parameter h — i/o/lOO km/sec/Mpc 
= 0.65, the resulting absolute magnitudes at peak bright- 
ness are -20.00, -19.54, -19.46, -19.45, -19.18 mag for 
U,B,V, R, /, respectively. The p and q vectors translate 
the shape-luminosity parameter A into a change in the 
SN la absolute magnitude. The p'^'f values (at peak 
brightness) vary among passbands from 0.6 to 0.8, and 
the q^'-^ vary from 0.1 to 0.9; therefore, intrinsically faint 
(bright) SNe have positive (negative) values of A. 

We use model vectors based on the procedure outlined 
in JRK07, but with two notable differences. First, the 
vectors have been re-evaluated based on our determina- 
tion of the dust parameter, Ry = 2.18. Since most of the 
nearby objects used in the training have low extinction, 
retraining with a different value of Ry has little effect 
on the vectors and therefore on the cosmological results. 
Note that the insensitivity of the MLCS2k2 training to 
the value of Ry does not imply that the estimated dis- 
tances for high-redshift SNe are insensitive to the value 
of Ry, especially since the latter samples include highly 
extinguished SNe. The impact of Ry on the cosmology 
results is presented in §9. 

The second change in the model vectors from JRK07 
involves adjustments to the M^'-^ that were developed 
during the course of the WV07 analysis of the ESSENCE 
data. For the model training with the nearby SN la sam- 
ple, it was assumed that the observed Ay distribution 
has the functional form of an exponential distribution 
convolved with a Gaussian centered at Ay = 0. However, 
the mlcs2k2 training process resulted in a convolution 
Gaussian that is not centered at zero; adjustments were 
made in the vectors such that the Gaussian is centered 
at zero. The main caveat in this procedure is that the 
selection efficiency for the nearby sample is small and un- 
known, and therefore it is not straightforward to model 
the observed Ay distribution in terms of an underlying 
population. The M*^'^ adjustments for UBVRI depend 
only on the passband and are independent of epoch. For 
the model vectors determined with Ry = 2.2, the mag- 
nitude adjustments relative to the values in JRK07 are 

gj^UBVBJ ^ _^o.050, +0.020, 0.0, -0.002, -0.033. (2) 

K-corrections transform the mlcs2k2 SN rest-frame 
Landolt-system magnitudes to the magnitudes of a red- 
shifted SN in an observed passband. K-corrections are 
computed following the prescription of Nugent, Kim, & 
Perlmutter (2002), which requires a SN spectrum at each 
epoch, the spectrum of a reference star, and the reference 
star magnitude in each passband. As explained in Ap- 
pendix A, we use a single template spectrum for each 
SN epoch and warp it to match the colors of the SN 
model. Since the Landolt photometry is not associated 
with a precisely defined set of filters, we use the stan- 
dard UBVRI and Bx filters defined by Bessell (1990) 
and apply a color transformation to obtain photometry 
in the Landolt system. That procedure is detailed in Ap- 
pendix B. In place of the traditional primary reference 
star Vega, we choose BD+17"4708 (Oke & Gunn 1983) 
as our primary reference because it has been measured 
by Landolt, it has a precise HST STIS spectrum (Bohlin 
2007) and it is the primary reference for SDSS photom- 
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etry. Wc have also carried out the analysis with Vega 
as the primary reference and include the difference as a 
systematic error. The primary magnitudes for each fil- 
ter system are given in Table 22 of Appendix C for both 
BD+17 and Vega. 

A light-curve fit determines the likelihood function £ 
of the observed magnitudes or fluxes as a function of four 
model parameters for each SN la: (i) time of peak lumi- 
nosity in rest-frame i3-band, to, (ii) shape-luminosity pa- 
rameter, A, (iii) host-galaxy extinction at central wave- 
length of rest-frame F-band, Ay, and (iv) the distance 
modulus, fi. The redshift (z) is accurately determined 
from the spectroscopic analysis, so it is not included as a 
fit parameter; the redshift uncertainty is included in the 
cosmology analysis (§ 8). For each supernova, the log of 
the posterior probability Ppost, or statistic, is given 
by 

= -2 In Ppost (to, A, Ay, Ai|data) 
= -21n£(data|to, A,^y,/x) - 21nPprior(-z, Ay, A) , 

(3) 

where PpHor is a Bayesian prior (see below), and the log- 
likelihood is given by 

-21nZ: . (4) 

[ j '''i.stat + ''^i, model J 

Here the index i runs over all measured epochs and 
observer-frame passbands, and F^'^^'^ is the observed flux 
for measurement i. The statistical measurement uncer- 
tainty, (Tstat , is estimated from the Scene Model Photom- 
etry as described in § 3.2. For the model uncertainty, 
fmodei, we usc thc diagonal elements of the mlcs2k2 co- 
variance matrix, which are estimated from the spread in 
the training sample of SNe. For example, at the epoch 
of peak brightness (to) these model errors are 0.11, 0.07, 
0.08, 0.10, 0.11 mag for U, B, V, R, I, respectively; thc un- 
certainties increase monotonically with time away from 
to- As explained below in the list of modifications, we 
do not use the off-diagonal mlcs2k2 correlations in this 
analysis. 

Since Ay is a physical parameter that is always posi- 
tive, and since it is not well constrained if thc peak signal- 
to-noise is low or if the observations do not span a large 
wavelength range, the mlcs2k2 fit includes a Bayesian 
prior on the extinction. The prior forbids negative values 
of Ay and encodes information about the distribution of 
extinction in SN host galaxies as well as the selection 
efficiency of the survey. Since there is degeneracy be- 
tween the inferred values of Ay and fx, the prior leads to 
reduced scatter in the Hubble diagram. For the nearby 
SN sample, which has high peak signal-to-noise for all 
objects, the prior has no impact on the Hubble scatter; 
for the other samples we employ, the prior reduces the 
Hubble scatter by a factor of 1.3-2. For this analysis, 
thc prior is defined to be 

Pprior{z,Ay,A)^P{Ay)P{A) X 

escarch(2. Ay, A)ecuts(z, Ay, A), (5) 

where P{Ay) and P(A) arc thc underlying SN la pop- 
ulation distributions of Ay and A, and wc assume that 
these distributions are independent of redshift. We de- 
termine them from SDSS-H SN data in § 7. For SNe 



passing the selection cuts, the parameter A is typically 
precisely determined by thc light-curve fit, so thc prior 
on A does not have a significant impact on thc inferred 
parameters. The functions esearch and ecuts are survey- 
dependent efficiency factors associated with the survey 
selection functions and with the sample selection cuts. 
The efficiencies are determined from Monte Carlo simu- 
lations in conjunction with the observed data distribu- 
tions for each survey in § 6. Tests with high-statistics 
simulations have verified that the prior in Eq. 5 leads to 
unbiased results for cosmological parameters. 

In the mlcs2k2 fit, the estimated value and uncer- 
tainty for each model parameter, e.g., the distance modu- 
lus fi, arc obtained by marginalizing the posterior (Eq. 3) 
over the three other parameters and taking the mean and 
rms of the resulting one-dimensional probability distri- 
bution. In the marginalization integrals, we usc 11 bins 
in each of thc parameters; this choice is dictated by the 
computational time required for the large number of sys- 
tematics tests (see § 9.1). We have compared results with 
11 and 15 integration bins and find excellent agreement. 

To implement the MLCS2k2 method, we have written 
a new version of the fitting package with several modifi- 
cations from JRK07: 

1. We fit in calibrated flux instead of magnitudes. In 
previous analyses using mlcs2k2, the fits were car- 
ried out using magnitudes, and data with SNR< 5 
were typically excluded in order to avoid ill-defined 
magnitudes associated with negative fiux measure- 
ments. The SNR cut results in a biased determi- 
nation of the shape-luminosity parameter A and 
therefore of the distance modulus fi. Fitting in flux 
enables a proper treatment of errors for all mea- 
surements and results in a negligible bias in A and 
fi, as determined from a simulation. This change is 
crucial for our analysis, since ^ 40% of the SDSS-II 
SN measurements (with — 15 < Trcst < +60 days) 
have SNR< 5. 

2. We have made two improvements to the treatment 
of K-corrections. First, we use the updated SN la 
spectral templates from Hsiao et al. (2007), which 
result in better consistency between the data and 
the best-fit mlcs2k2 model for observer-frame fil- 
ters that map onto rest-frame i?-band. Second, we 
have improved the spectral warping used for K- 
corrections as explained in Appendix A. 

3. The mlcs2k2 model includes off-diagonal covari- 
ances in thc model magnitudes to account for 
brightness correlations between different epochs 
and passbands; in this analysis, wc ignore the off- 
diagonal covariances for two reasons. The primary 
reason is that the MLCS2k2 model covariances ap- 
pear to display unphysical behavior. The correla- 
tion coefficient pij = cov{i, j) / (Ji<Tj between epochs 
i and j decreases discontinuously from unity at ti = 
tj (Fig. 4): the correlation between epochs sepa- 
rated by only one day is weak, 0.2 < pt,t+i < 0.8, 
and thus docs not penalize (via x^) random varia- 
tions of 0.1 mag over one-day time-scales. The 
observed smoothness of high-quality SN la light- 
curve data rules out such large intrinsic fiuctua- 
tions, suggesting that random instrumental noise 
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may have been included in the model covariance 
matrix. The impact of the off-diagonal covariances 
on determination of the cosmological parameters 
{w and flu) from the SN data is much smaller 
than the statistical uncertainties. Second, there 
is a subtle limitation when measurements at the 
same epoch in two observer- frame passbands /i , /2 
are matched onto the same rest-frame filter /' us- 
ing Arest = Aobs/(l + z) for each passband. In the 
MLCS2k2 model, there is an artificial 100% cor- 
relation between the two rest-frame model magni- 
tudes. This feature arises for the observed ugriz 
filters used by SDSS-II and SNLS, but does not ap- 
pear for the Bessell filters used in the nearby and 
ESSENCE samples. 

4. We have extensively modified the prior (Eq. 5). 
The mlcs2k2 prior in JRK07 is intended to re- 
flect the true distribution of Ay. In analyzing 
the ESSENCE data, WV07 used a different Ay 
prior and multiplied it by a simulated efficiency 
that depends upon extinction, intrinsic luminosity, 
and redshift. In our analysis we use more detailed 
Monte Carlo simulations (§ 6) of each data sample 
to estimate the survey efficiencies that are incor- 
porated into the priors, and we use the SDSS-II 
SN data sample to determine the underlying Ay 
distribution. 

5. In mlcs2k2, the reddening parameter Ry is 
treated as a fixed global parameter. In JRK07 and 
WV07, Ry was set to the average Milky Way value 
of 3.1. In our analysis, we use Ry = 2.18 ± 0.50 as 
empirically determined from the SDSS-II SN sam- 
pie (§ 7). 

Some example fits for SDSS-II SN light curves using 
the modified version of mlcs2k2 are shown above in 
Fig. 1. Figure 5 shows the average fractional residuals 
between the MLCS2k2 model light curves and the data for 
the SNe in each survey and for each rest-frame UBVR 
passband. The overall data-model agreement is good, ex- 
cept for some late-time epochs and J7-band. The [/-band 
residuals are discussed later in more detail (§ 10.1.3). 
Fig. 6 shows the fit parameters Ay and A vs. redshift 
for SNe in the different surveys. The impact of the prior 
requiring Ay > is immediately evident in the top-left 
panel. If we split each SN sample at its median red- 
shift, the average Ay for the lower-redshift SNe is larger 
than for the higher-redshift SNe; this Ay-difference is 
0.1 mag for the nearby sample, and ~ 0.05 mag for 
the other SN samples. The prior discussed above ac- 
counts for this redshift-dependent shift. The right panel 
in Fig. 6 shows the fitted Ay vs. redshift using a flat 
prior, P{Ay) = P{A) = 1 in Eq. 5. Although there are 
many SNe with Ay < 0, in § 7.3 we show that an un- 
derlying extinction distribution with Ay > 0, combined 
with measurement uncertainties, is consistent with the 
"negative- Ay" distribution obtained from fitting with a 
flat prior. Since A is well constrained by the light-curve 
fits, the A distribution with a flat prior is very similar to 
that using the nominal prior. 

We have checked the results of the modified mlcs2k2 
fitter with the distance estimates derived by WV07 for 



the ESSENCE and nearby SN samples. For this compar- 
ison, we use the WV07 extinction prior and efficiency, as 
described in their Eqs. 2 and 3. The WV07 efficiency 
function accounts for missing SNe at high redshift and 
for the bias arising from using only measurements with 
SNR > 5; for comparison, we therefore use the same 
SNR cut. The modified fitter is run in a mode that repli- 
cates the original mlcs2k2 fitter, with two exceptions: 
First, as noted above, we fit in flux instead of magnitude. 
Second, the K-corrections use the average spectral tem- 
plate of Hsiao et al. (2007), while WV07 used a library 
of spectra and interpolated K-corrections to the desired 
epoch. To compare our "nearby-|-ESSENCE" analysis 
with WV07, we fit light curves for the 45 nearby SNe la 
(0.015 < z< 0.1) and 57 ESSENCE SNe la analyzed by 
WV07. The rms scatter between our fitted distance mod- 
uli (fi) and those from WV07 is 0.03 mag and 0.05 mag 
for the nearby and ESSENCE samples, respectively. Our 
marginalized value for the dark energy equation of state 
parameter w, using the SDSS BAO prior (see § 8), agrees 
to within 0.01 with the result of WV07. 

We stress that this comparison with WV07 is a consis- 
tency check of our version of MLCS2k2 relative to previ- 
ous versions. When we analyze the present SN samples 
with mlcs2k2 (§ 10), our different prior and mlcs2k2 
model parameter values result in cosmological parameter 
estimates that differ significantly from those of WV07, as 
discussed in § 10.1.4. 
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Fig. 4. — For the mlcs2k2 model, correlation coefficient PAt,o 
between B-band epoch at peak brightness (to) and time At = 
t — to, where PAt,o = cov(Af, 0)/(TAto"0, as a function of At, and 
correlation coefficient pt— tio.tio between epoch at 10 days past 
peak {tio) and time t — tio. The spikes at and 10 days correspond 
to the requirements po,0 = 1 a^nd PiO.lO = 1- 



5.2. SALT-II Fitting Method 

The SALT-ii light curve fitting method (Guy et al. 
2007) has been developed by the SNLS collaboration. 
The SALT-II model employs a two-dimensional surface in 
time and wavelength that describes the temporal evolu- 
tion of the rest-frame spectral energy distribution (SED) 
for SNe la. The temporal resolution of the model is 1 day, 
and the wavelength resolution is 10 A, allowing accurate 
synthesis of model fluxes to compare with photometric 
data. The model is created from a combination of photo- 
metric light curves and hundreds of SN la spectra. When 
there are measurement gaps in the spectral surface, the 
unmeasured regions of the SED are determined from in- 
terpolations of the measured regions. The photometric 
data is mostly from the nearby sample (JRK07) but also 
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includes higher-redshift data (z > 0.1) to better con- 
strain the rest-frame ultraviolet behavior of the model. 
For a complete list of SN light curves and spectra used 
for training, see Table 2 in Guy at al. (2007). 

In SALT-ii, the rest-frame iiux at wavelength A and 
time f (t = at B-band maximum) is modeled by 



dF„ 



dX 



- (i , A) = xo X [Mo (t , A) + x-i X Ml (i , A)] 
X exp[c X CL{\)] . 



(6) 

Mo{t,X), Mi(i, A), and CL{\) are determined from the 
training process described in Guy et al. (2007). The 
Mo surface represents the average spectral sequence, and 
is very similar to the sequence of average spectral tem- 
plates (Hsiao et al. 2007) that we use for the mlcs2k2 
K-corrections. Mi is the first moment of variability 
about this average, accounting for the well-known corre- 
lation of both peak brightness and color with light-curve 
shape, and xi is the stretch parameter, the analog of 
the mlcs2k2 A parameter. CL{X) is the mean color 
correction term, and c is a measure of SN la color. Al- 
though the color variation is not explicitly attributed to 
dust extinction, in the optical region CL(A) is reason- 
ably well approximated by the CCM89 extinction law 
with Rv ~ 2. In the UV region, CL{X) exceeds the 
CCM89 extinction by about 0.07 mag. 

The spectral-time surfaces are defined for rest-frame 
times —20 < Ticst < +50 days relative to the time 
of maximum brightness, and for rest-frame wavelengths 
that span 2000 to 9200 A. We use the salt-ii spec- 
tral surfaces obtained from retraining the model using 
Bessell-filter shifts based on HST standards, as discussed 
in Appendix B (Table 21), but otherwise using the same 
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Fig. 5. — Data-model fractional residuals as a function of rest- 
frame epoch in 5-day bins for mlcs2k2 light-curve fits. The rest- 
frame passband and SN sample are indicated on each plot. Mea- 
surements with SNR< 6 are excluded, and error bars indicate the 
rms spread. For SNLS, the residuals are shown only for SNe with 
z < 0.5 as explained later in §10.2.4. F^ata (Fmodcl) is the SN 
flux from the data (best-fit mlcs2k2 model). Vertical dashed lines 
indicate epoch of peak brightness (Ti-cst = 0); horizontal dashed 
lines indicate Fjata = Fmodd- 
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VS. redshift, for the different SN la samples indicated on the plot. 
Right panels: fitted A vs. redshift. Upper panels are from fit with 
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technique and data as described in Guy et al. (2007). The 
UBVRI magnitudes for the primary reference Vega are 
taken from Fukugita et al. (1996): these are 0.02, 0.03, 
0.03, 0.03, 0.024, respectively, and are slightly different 
from those used to crosscheck the mlcs2k2 method. Al- 
though the wavelength coverage of the spectral surface 
is rather broad, the SALT-ii model includes only those 
observer-frame passbands for which 2900 < A//(l-|-2;) < 
7000 A, where A/ is the mean wavelength of the filter 
and z is the SN la redshift. 

To compare with photometric SN data, the observer- 
frame flux in passband / is calculated as 



FLit) = i^ + ^) 



dX' 



X'^{t,X')Tf{X'{l + z)) 



(7) 

where T^{X) defines the transmission curve of observer- 
frame passband /. For the SDSS-II, ESSENCE, SNLS, 
and HST samples, T-f{X) is provided by each survey. For 
the nearby sample, T^{X) is given by the Bessell (1990) 
U BVRI filter response curves, with wavelength shifts as 
described in Appendix B and listed in Table 21. 

The model uncertainty accounts for the covariance be- 
tween Mo{t, A) and Mi{t, A) at the same epoch and wave- 
length. Although spectral covariances between different 
epochs and wavelengths are not considered, the model 
does account for covariances between integrated fluxes 
at different epochs within the same filter. Each SN la 
light curve is fit separately using Eqs. 6 and 7 to deter- 
mine the parameters xq, xi, and c. However, the salt- 
ii light-curve fit does not yield an independent distance 
modulus estimate for each SN. As discussed in § 8.2, the 
distance moduli are determined as part of a global fit to 
an ensemble of SN light curves in which cosmological pa- 
rameters and global SN properties are also determined. 
The SALT-II fits do not include informative priors on the 
fit parameters or the effects of selection efficiencies. We 
correct the SALT-ii results for selection biases using a 
Monte Carlo simulation (see §6 and 9.2). 
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In most cases, the SALT-ii light-curve fits arc quahta- 
tively very similar to the mlcs2k2 fits on a per-object 
basis. The average rest-frame light curve residuals for the 
SALT-II fits are shown in Fig. 7 for each survey for filters 
UBVR; note that the [/-band residuals for the nearby 
SNe show some discrepancy, as will be discussed later. 
Fig. 8 shows the fitted values for the color parameter c 
and stretch parameter xi versus redshift for SNe in the 
different surveys. 

As a crosscheck on our use of the pubhc SALT-ii code, 
we have compared our fits of the 71 SNe la from Astier 
et al. (2006) to fits done by the developer (J. Guy, private 
communication) and find good agreement. The mean 
difference in the color (c) is 0.003 ± 0.003, with an rms 
dispersion of 0.02, and the mean difference in the shape- 
luminosity parameter (xi) is — 0.014±0.021, with an rms 
dispersion of 0.18. The slight differences are attributable 
to the use of different versions of the code and of the 
error (dispersion) map. 
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Fig. 7. — Data-model fractional residuals as a function of rest- 
frame epoch in 5-day bins, for SALT-II light-curve fits. The rest- 
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5.3. Comparison of mlcs2k2 and SALT-ii Fitters 

We end this section by briefly comparing and contrast- 
ing the SALT-II and mlcs2k2 methods. The mlcs2k2 
rest-frame model for the intrinsic SN brightness is de- 
fined in discrete UBVRI passbands corresponding to 
the Landolt system. For each SN light curve, a com- 
posite SN spectrum is warped based on the model fit to 
the observed SN colors at each epoch, and the warped 
spectrum is used to perform the K-corrections needed 
to transform the rest-frame model to the observer-frame 
fluxes. The SALT-ii model uses a composite SN spectrum 
that depends on both the epoch and intrinsic luminosity 
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Fig. 8. — Left panel: SALT-ll fitted color values (c) vs. redshift, 
for the SN la samples indicated on the plot. Right panel: fitted 
stretch parameter values, xi, vs. redshift. 



as well as an epoch-independent color term. This spec- 
trum is used to model the rest-frame fluxes. Both the 
SALT-II and mlcs2k2 models are trained using nearby 
SN la data, but SALT-ii training also includes higher- 
redshift data that reduces the dependence on the nearby 
SN sample and provides better constraints on the rest- 
frame ultraviolet regions of the spectrum. 

The SALT-II parameters xi and c are analogous to the 
mlcs2k2 parameters A and Ay- The parameters xi and 
A are essentially equivalent in describing the correlation 
between SN light-curve shape and brightness, but c and 
Ay have different meanings. mlcs2k2 assumes that all 
intrinsic SN color variations are captured in the model by 
the light-curve shape-luminosity correlation and that any 
additional observed color variation is due to reddening by 
host-galaxy dust. The color term (c) in SALT-ii describes 
the excess color (red or blue) of a SN relative to that 
of a fiducial SN with fixed stretch parameter xi. The 
excess color could be from host-galaxy extinction, from 
variations in SN color that are independent of xi , or from 
other effects, and SALT-ii does not attempt to separate 
these effects. SALT-ii uses c to reduce the scatter in the 
Hubble diagram in a manner analogous to the use of xi. 
The global SALT-ii parameter /3, defined below in § 8.2, is 
the analog of the global MLCS2k2 dust parameter Rb = 
Rv + 1; one expects /3 ~ Rb excess color variation 
is purely due to host-galaxy extinction. The SALT~ii f3 
parameter is determined from the global fit to the Hubble 
diagram for the entire SN la sample under analysis; we 
determine the mlcs2k2 Ry parameter by modeling the 
observed colors of a specific subset of the SN data (§ 7.2). 

Concerning correlations among model parameters, 
mlcs2k2 and SALT-ii treat different aspects. The 
mlcs2k2 model includes covariances between different 
epochs and passbands, but we have excluded the off- 
diagonal covariances as explained in § 5.1. The SALT- 
II model includes covariances between integrated fluxes 
at different epochs within the same passband, but co- 
variances between passbands are not considered. SALT-ii 
also includes the covariance between the spectral surfaces 
Mo{t, A) and Mi{t, A) at each epoch and wavelength bin 
(Eq. 6), but it does not include covariances between dif- 
ferent epochs and passbands. 

Within the mlcs2k2 framework, each light-curve fit 
yields an estimated distance modulus along with its es- 
timated error, independent of cosmological assumptions. 
By contrast, in SALT-ii the distance modulus estimate 
for a given supernova is based on a global fit to the en- 
semble of supernovae within a parametrized cosmological 
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model (see § 8.2). A result of this global minimization 
in SALT-ii is that a distance modulus bias in a particu- 
lar redshift range, such as could arise from including a 
poorly calibrated filter, will induce a bias over the entire 
redshift range of the sample. For the determination of 
cosmological parameters, this tends to reduce the sen- 
sitivity to systematic problems and hence can lead to 
smaller systematic uncertainties. However, this reduced 
sensitivity can also make biases more difficult to identify. 
An explicit example of this is described in § 10.2.4. 

Fitting with MLCS2k2 usually incorporates a Bayesian 
prior (Eqs. 3-5) that reduces the scatter in the Hubble 
diagram by incorporating information about the under- 
lying Ay distribution and the survey efficiencies. The 
prior, and the resulting Hubble scatter, do not depend 
on cosmological parameters. Because of the assumption 
that excess color variation is due to extinction by dust, 
the prior excludes values of Ay < 0. In mlcs2k2, SNe 
with very blue apparent colors (bluer than the template) 
are assigned Ay — 0, and the data-model color discrep- 
ancy is attributed to fluctuations. In SALT-ii, apparently 
blue SNe are assigned negative colors (c < 0) that result 
in larger luminosities and distance moduli compared to 
mlcs2k2. 

Within the SALT-ii framework, scatter in the Hubble 
diagram is explicitly minimized by simultaneously ad- 
justing global SN parameters along with the cosmologi- 
cal parameters; this minimization is described in § 8.2. 
In contrast to MLCS2k2, the SALT-II Hubble scatter de- 
pends on the cosmological parameters, and there is no 
mechanism to account for the survey efficiency directly 
in the fits. To correct for biases related to the survey 
efficiencies (§ 8.2), we use the Monte Carlo simulations 
described in § 6. 

The MLCS2k2 and SALT-ii light curve fit-residuals can 
be visually compared in Figures 5 and 7; the data and 
models are consistent for rest-frame passbands BVR, but 
there are discrepancies for [/-band in both cases. We 
address this issue in more detail in §10.1.3 and 10.2.4, and 
we compare the mlcs2k2 and SALT-ii results explicitly 
in § 11. 

6. MONTE CARLO SIMULATION: DETERMINING THE 
SELECTION EFFICIENCY 

All surveys suffer from incompleteness and selection ef- 
fects of various kinds. Supernovae that are intrinsically 
subluminous or highly extinguished by dust have less 
chance of being included in a flux-limited sample than 
more typical SNe. In addition, with limited spectroscopic 
resources, higher priority may be given to SN candidates 
with the best chances of yielding reliable identiflcations, 
e.g., by focusing on events that appear well separated 
from the host galaxy or for which the host has either low 
surface brightness or early-type colors and morphology 
that suggest low dust content. These selection effects 
become more pronounced at the high-redshift end of a 
survey, where only the brightest, unextinguished SNe will 
satisfy selection cuts. If SN la brightness were a perfectly 
standardizable distance indicator, such selection effects 
would not be an issue for cosmological analysis. However, 
intrinsic variations in SN brightness, photometric errors, 
and uncertainties in estimating host-galaxy dust extinc- 
tion lead to signiflcant uncertainties and possible biases 
in distance estimates, particularly for SNe observed with 



low signal-to-noise. In order to extract unbiased cos- 
mological parameter estimates, biases must either be re- 
duced to an acceptably small level by the analysis pro- 
cedure or else a correction scheme must be adopted. 

We have developed detailed Monte Carlo simulations 
of the different SN surveys in order to determine the sur- 
vey selection (or efficiency) functions and their impacts 
on SN distance estimates for both mlcs2k2 and SALT-ii. 
The simulations also enable us to verify the estimates of 
systematic errors due to uncertainties in the light-curve 
model parameters. The simulated efficiency is a major 
component in the mlcs2k2 fit prior discussed above in 
§ 5.1. Determining the host-galaxy extinction depen- 
dence of the efficiency is critical for the mlcs2k2 method, 
because the extinction is often poorly determined from 
the data. For the SALT-II method, the simulation and ef- 
ficiency play no direct role in the fitting, but they enable 
us to estimate and correct for biases in the cosmological 
parameters as desribed in § 8.2. 

Ideally, survey simulations would be based on artificial 
SNe la embedded into survey images, as was done during 
the SDSS-II SN survey to monitor the efficiency of the 
search pipelines (see § 2). We do not have access to the 
images for the other surveys, and full image-level sim- 
ulations would require a large amount of computing to 
perform the many variations that arc needed for the anal- 
ysis. We have instead developed a fast light-curve simu- 
lation''° that is based upon actual survey conditions and 
that therefore accounts for non-photometric conditions 
and varying time intervals between observations due to 
bad weather. At each survey epoch and sky location, 
the simulation uses the measured point spread function 
(PSF), zero point, CCD gain, and sky background to de- 
termine the noise and to convert the simulated model 
magnitudes into CCD counts. The simulation also incor- 
porates a model for host-galaxy light and dust extinction. 
We have obtained the necessary observational informa- 
tion for the SDSS-II, ESSENCE, SNLS, and HST sur- 
veys to carry out these detailed simulations. The nearby 
SN la sample is a heterogeneous sample collected over 
many years by different observers and telescopes, and 
we do not have the information needed to make detailed 
simulations of this sample. 

Here we describe the simulation within the context of 
the mlcs2k2 light-curve model and comment on the dif- 
ferences needed to simulate light curves in the SALT-ii 
model. We select a random SN redshift from a power- 
law distribution, dN/dz ~ (1 z)f^, with /3 = 1.5 ± 0.6, 
as determined by our recent analysis of the SN la rate 
(Dilday et al. 2008). A SN la luminosity parameter A 
and host-galaxy extinction Ay are selected from under- 
lying distributions that we have inferred from our data 
(§ 7.3). The mlcs2k2 model is used to convert A into 
rest-frame UBVRI magnitudes. These generated SN la 
magnitudes are increased according to the selected Ay 
and the CCM89 extinction law using Ry = 2.18, as de- 
termined in § 7.2. The reddened UBVRI magnitudes 
are K-corrected into observer-frame magnitudes. A ran- 
dom sky coordinate is selected from the survey area, and 
Milky Way extinction is applied based on the maps of 

The simulation, along with the light curve fitters described 
in § 5, are publicly available in a software package called SNANA: 
http : //www. sdss . org/supernova/SNANA .html 
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SchlegeL Finkbcincr, & Davis (1998). A random date for 
peak brightness is selected from the smvcy time frame, 
and all observed epochs at the selected sky coordinate 
are identified from the actual survey observations. For 
each observation epoch, the measured survey zero point 
is used to convert the simulated magnitude into a simu- 
lated flux. For simulations based on the SALT-ii model, 
the mlcs2k2 parameters A and Ay are simply replaced 
by the corresponding SALT-ii parameters {xi,c), drawn 
from empirical distributions. 

The simulated noise for each epoch and filter includes 
Poisson fluctuations from the SN la (signal) flux, sky 
background, CCD read noise, and host-galaxy back- 
ground. The signal noise is based on the number of 
CCD photoelectrons calculated from the simulated flux. 
The sky background is computed from the measured 
sky background per pixel, which is summed over an ef- 
fective aperture based on the measured PSF at that 
survey epoch and sky coordinate. For SN redshifts 
•ZSN < 0.4, noise from the host galaxy is simulated by 
associating the SN with a host from the SDSS galaxy 
photometric redshift catalog (Oyaizu et al. 2008), ran- 
domly selected such that Zgai ^ zsn- From the SDSS 
DR5 (Adelman-McCarthy et al. 2007) photoPrimary 
database (Stoughton et al. 2002), we use the fitted ex- 
ponential surface brightness profile in j'-band as a prob- 
ability distribution from which the SN position within 
the galaxy is randomly selected, i.e., we assume that 
the SN la rate within a galaxy is proportional to the 
local r-band luminosity. The host-galaxy background is 
computed by integrating the exponential galaxy model 
within the same effective aperture that is used for the 
sky noise. The exponential profile is not appropriate for 
early-type galaxies, but this model is meant only as an 
estimate of the range of host-galaxy background light 
expected. The host-galaxy noise exceeds the sky noise 
for only ~ 10% of the simulated SNe la with zsn < 0.4. 
For redshifts greater than 0.4, the lack of simulated host- 
galaxy noise is not significant, because the sky noise is 
dominant at these higher redshifts 

There remain two important aspects of the simulation 
that are less well-defined and therefore more difficult to 
model: (i) intrinsic variations in SN la properties, be- 
yond the shape-luminosity correlation, that lead to (so 
far) irreducible scatter in the Hubble diagram; and (ii) 
search-related inefficiencies beyond those due to photo- 
metric signal-to-noise and selection cuts, e.g., those asso- 
ciated with spectroscopic selection. Below, we describe 
our modeling of these features in the simulation. 

6.1. Simulating Variations of Intrinsic SN Brightness 

Using the MLCS2K2-based simulation described above, 
the resulting scatter in the Hubble diagram for SDSS-II 
SNe at z < 0.15 is only 0.06 magnitudes, well below the 
observed scatter of ~ 0.15 mag. To make the model more 
realistic, we introduce intrinsic fluctuations in the sim- 
ulated luminosity. The models for intrinsic fluctuations 
described below are empirically determined to match the 
observed Hubble scatter and are not based on a physical 
model. 

We have implemented two models of intrinsic SN vari- 
ations. The default method we use for mlcs2k2, called 
"color-smearing," introduces an independent fluctuation 
in each passband, and the fluctuation is the same for 



all epochs within each passband. A random number rfi 
from a unit-variance Gaussian distribution is chosen for 
each rest- frame passband /'. A magnitude fluctuation, 
Srufi ~ rfia'j.-,, is added to the generated magnitude at 

all epochs, where Cj, is the magnitude uncertainty at 
peak brightness given by the mlcs2k2 model in pass- 
band /'. In this method, the intrinsic model colors are 
randomly varied by typically ~ 0.1 mag. Since the sim- 
ulated color variations are the same at all epochs, this 
model does not respect the Lira law (Phillips et al. 1999), 
the empirical observation that intrinsic SN la colors have 
smaller variations at epochs later than about 2 months 
after explosion; this defficiency has a negligible impact 
on our analysis because our requirements of good light 
curve coverage make our simulated efficiencies insensitive 
to the magnitudes at such late epochs. 

The second model of intrinsic variation, which we use 
for the SALT-II method, and as a crosscheck for the 
mlcs2k2 method, is called "coherent luminosity smear- 
ing:" a coherent random magnitude shift, typically ^ 
0.15 magnitudes, is added to all epochs and passbands. 
In the coherent smearing method, the intrinsic model 
colors are not varied. 

A caveat in our implementation of intrinsic luminosity 
variations is that the mlcs2k2 simulation and fitter use 
different models of intrinsic fluctuations and covariances. 
Although the mlcs2k2 model includes a full covariance 
matrix, we argued in § 5.1 that these covariances do not 
accurately reflect intrinsic correlations. Since we use only 
the diagonal elements of the mlcs2k2 covariance ma- 
trix in the fitter, a literal translation for the simulation 
would be to implement a random intrinsic fluctuation of 
~ 0.1 mag independently for each epoch and passband. 
Since the observed smoothness of high-quality SN la light 
curve data rules out such large intrinsic epoch-to-epoch 
fluctuations, we have chosen the above methods to simu- 
late intrinsically smooth light curves. In future training 
of SN la light curves, it will be desirable to extract a 
model of intrinsic fluctuations and covariances that can 
be used consistently in both the light-curve fitter and the 
simulation. 

6.2. Simulation of Survey Search Efficiency 

The final step in the simulation is to model losses re- 
lated to the SN search. In particular, we need to account 
for SNe la that would have passed the light-curve selec- 
tion cuts of § 4 had they been identified, but that were 
missed due to inefficiencies in the SN search. 

Search-related losses come from the following sources: 
(a) the image-differencing pipeline can fail to detect ob- 
jects with very low signal-to-noise as well as objects with 
nearby artifacts such as a diffraction spike; (b) humans 
tasked with evaluating objects detected in subtracted im- 
ages may not correctly identify them as possible SN can- 
didates; (c) software used in spectroscopic targeting to 
fit light curves and photometrically classify SN candi- 
dates identified by humans may not correctly classify 
all SNe la; and (d) due to limited resources for spec- 
troscopic observations, not all photometrically identified 
SN la candidates will be targeted spectroscopically or re- 
sult in a spectrum with sufficient signal-to-noise to con- 
firm the SN type and determine its redshift. 

We define the overall survey efficiency, or survey selec- 
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tion function, as Csurvcy = Cscarch X Ecuts, whcre the search 
efficiency is further decomposed as egcarch = Csubtr x Cspoc- 
Here, egubtr describes the net search efficiency of the 
image-subtraction pipehne corresponding to step (a) 
above. The term egpoc describes the combination of steps 
(b), (c), and (d), which depends in part on human judg- 
ment for each SN, and ecuts is the factor associated with 
the final selection cuts described in § 4. These decompo- 
sitions of efficiency components are convenient because, if 
sufficient information about the search is available, then 
ecuts and egubtr can be reliably simulated. By contrast, it 
is usually impossible to directly simulate egpoc because it 
involves complex decision-making under varying circum- 
stances by many people involved in spectroscopic obser- 
vations. However, if the spectroscopic efficiency is nearly 
100% below some redshift for a given survey, then we can 
model Espcc at higher redshifts by comparing observed 
distributions of SNe properties (including redshift) to 
simulated distributions expected for a spectroscopically 
complete survey. In summary, our philosophy is to make 
as detailed a model of the survey efficiency components 
as the data and survey information allow, and to use 
data-simulation comparisons to constrain the other com- 
ponents. 

The different components of the efficiency arc likely to 
be correlated, so that the overall efficiency is not really a 
simple product as defined above. As an extreme exam- 
ple, if both the search and selection cuts required only 
that the maximum signal-to-noise ratio was greater than 
5, then Csearch = Ecuts, and the combined efficiency would 
be the same (not the product of the two). As discussed 
below for the SDSS-II survey, we can simulate the com- 
bined effect of the image-subtraction efficiency (esubtr) 
and selection cuts (ccuts)- To simplify notation we write 
the combined efficiency as a product, fgubtr x fcuts, but 
it should be understood that this product really refers 
to the combined efficiency taking into account all cor- 
relations. The espoc term is treated differently in that 
it is defined as an independent efficiency function that 
multiplies the other efficiency terms. 

For the SDSS-II, ESSENCE, SNLS, and HST samples, 
we have obtained detailed observing conditions, and the 
simulation accurately describes the effects of sample se- 
lection criteria, Ccuts- For the SDSS-II, we determine 
Csubtr using information from the fake SNe la that were 
inserted into the images during the survey (see below). 
For the non-SDSS samples, we do not have access to 
the pixel-level data, so we set Cgubtr = 1 and absorb 
the image-subtraction pipeline efficiency into the spectro- 
scopic efficiency, egpec = Csearch- For the nearby sample, 
we do not even have the information needed to determine 
the impact of the selection cuts; we therefore absorb all 
sources of inefficiency into the spectroscopic efficiency, 

^spec — ^survey- 

As noted above, for the SDSS-II, fake SNe la are used 
to infer Csubtr as a function of signal-to-noise ratio (SNR), 
as shown in Fig. 7 of Dilday et al. (2008). These efll- 
ciency curves are used by the light-curve simulation to 
probabilistically determine which measurements in which 
passbands result in detections; a single-epoch detection 
in any two of the three {gri) passbands is considered to be 
an object. An object found at two or more epochs results 
in a supernova candidate and is considered to be discov- 
ered by the image-subtraction pipeline. For the SDSS-II, 



the image-subtraction pipeline efficiency (egubtr) is com- 
plete up to a redshift of z '-^ 0.2 and then drops gradually 
to about 60% at z ~ 0.4. 

For the HST simulation, we use the single-epoch search 
efficiency as a function of magnitude, as described in 
Strolger et al. (2005) and Strolger & Riess (2006). The 
search was done with the F850LP_ACS passband (mean 
wavelength: 9070 A); it is fully efficient down to mag ~ 
23, and the efficiency drops to half at mag ~ 26. In- 
corporating this single-epoch efficiency profile into our 
simulation and requiring any one detection to discover a 
supernova, the simulated efficiency for finding SNe la is 
100% up to z ~ 1, and drops to about 65% at z = 1.6. 

6.2.1. Modeling Spectroscopic Efficiency 

Although selection effects associated with the pho- 
tometry can be simulated based on available informa- 
tion, modeling the spectroscopic efficiency for each sur- 
vey is more of a challenge. The SN redshift distribu- 
tions N(z) for the non-HST data and simulations are 
shown in Fig. 9. With the exception of the nearby SN 
sample, the simulations include losses from selection re- 
quirements (§ 4); for SDSS-II, losses from the image- 
subtraction pipeline are also included. For the nearby 
sample, the simulation is based on SDSS-II observing 
conditions and is chosen to be 100% efficient, since we do 
not have the information to model the effects of photo- 
metric selection requirements in that case; the simulated 
redshift distribution in narrow bins of redshift therefore 
scales as N{z) (ignoring evolution of the SN la rate 
at low redshift). 

Below some cutoff redshift, Zcut, the spectroscopic ef- 
ficiency is assumed to be ~ 100%; for z > Zcut the data- 
simulation discrepancy in f^iz) is taken to be due to 
spectroscopic inefficiency. For the SDSS-II and SNLS 
samples, the cutoff redshifts are estimated to be 0.15 
and 0.65, respectively, based on the levels of complete- 
ness given in their SN la rate measurements (Dilday et al. 
2008; Neill et al. 2006). For the ESSENCE sample, we es- 
timate Zcut ~ 0.45 based on visual inspection of the red- 
shift histogram in Fig. 9. For the nearby (z < 0.1) sam- 
ple, there is no redshift below which the spectroscopic 
completeness is assumed to be 100% (see discussion be- 
low). All four of these samples show significant spec- 
troscopic inefficiency at the upper ends of their redshift 
ranges. 

The data-simulation redshift comparison for the HST 
sample is shown in Fig. 10. Although the uncertainty in 
the SN la rate at these high redshifts is large, there is 
no evidence for significant spectroscopic inefficiency: the 
data-simulation redshift comparison is consistent with 
the claim that there are no significant losses for z < 1.4 
(Riess et al. 2007). We therefore assume that Cgpcc — 1 
for the HST survey, and this sample is not included in 
the efficiency discussion below. 

The spectroscopic efficiency as a function of redshift is 
defined to be espec = esurvey/(esubtr x Ccuts), which cor- 
responds to the ratios of the data and simulation his- 
tograms in Fig. 9. For each sample, this efficiency can 
be fitted by an exponential function, 

espcc(z) = Co exp[-(z - 2cut)/Ci] for z > Zcut (8) 
espGc(2:) = l for z < Zcut , 
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Fig . 9 . — Comparison of redshift distributions for data (dots) and 
simulations (histograms) for the nearby (Jha et al. 2007), SDSS-II 
(2005 data only), ESSENCE (Wood-Vascy et al. 2007), and SNLS 
(Astior et al. 2006) samples. The nearby sample is shown on a 
logarithmic scale, while the other vertical scales are linear. The 
simulations include losses from selection requirements (§ 4), and, 
for the SDSS-II, losses from the image-subtraction pipeline. The 
shaded regions show uncertainties in the simulated distributions 
due to uncertainty in the SN la rate vs. redshift (/3 = 1.5 it 0.6 
for dN/dz ~ (1 + z)f^). Except for the nearby sample, the sim- 
ulated distributions are scaled such that the integrated number 
of simulated SNe left of the cutoff redshifts (indicated by vertical 
arrows) match the data. Vertical error bars show the statistical 
uncertainty; horizontal bars show the bin-size. 
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Fig. 10. — Redshift distribution for HST sample (Riess et al. 
2007); data are shown by dots and the simulation by a his- 
togram. The shaded region shows the uncertainty from the 
redshift-dependence of the SN la rate (see Fig. 9 caption), but 
with the rate uncertainty doubled for z > 1. 

where i^o are determined from the fit. The dis- 

continuity at Zcut is motivated by the SDSS-II survey, 
for which our spectroscopic observation strategy targeted 
z < 0.15 candidates with very high priority. We have no 
evidence that the other surveys should have discontinu- 
ities in their redshift distributions, but the model above 
provides a reasonably accurate representation for the red- 
shift distributions of the ESSENCE & SNLS samples as 
well. This functional form is adopted for computational 
convenience in generating large Monte Carlo samples. 



Although the redshift dependence of the spectroscopic 
efficiency has now been estimated, we know that espcc(-z) 
depends in reality on a variety of factors, not simply red- 
shift, and we must model its dependence on those factors 
in order to properly model the selection function and as- 
sociated biases. As noted above, since spectroscopic tar- 
geting is a complex process, the underlying mechanism 
determining egpoc is not easily characterized. The sim- 
plest possibility would be that egpcc depends purely on 
redshift: in this case, e^^^^ = espoc(^) would be given by 
Eq. 8, and there would be no impact of spectroscopic 
efficiency on the survey bias, since redshift is precisely 
measured for each SN in the samples we consider. How- 
ever, this model is not a priori likely: it is more probable 



that 



-spec 



depends explicitly on both redshift and appar- 



ent SN la brightness, the well-known Malmquist bias. 

To consider this alternative, we define a "magnitude 
dimming" parameter A^dim as the difference between the 
simulated rcst-framc magnitude and the magnitude of 
the brightest possible SN la at peak light in rest-frame V 
band (for MLCS2k2) or in rest-frame B band (for SALT- 
ii). For the mlcs2k2 and SALT-ii models, A^dim is given 

by 



A^dim(MLCS2k2): 



Xdim(SALT-II): 



-9°'^(A^-A?, ^ 

:/3(c 



Arcf) 



Crcf) 



a{xi 



Xl. 



(9) 
of). (10) 



For mlcs2k2. Ay is the host-galaxy extinction in V- 
band, p"'^ and g"'^ are model parameters at the epoch 
of peak brightness in F-band (see Eq. 1), and Arcf = 
—0.3 corresponds to nearly the most intrinsically lumi- 
nous SN la in the training sample. For the SDSS-II, 
ESSENCE, and SNLS samples, 9 = 1; for the nearby 
sample, we set = 0, because the A distribution for the 
nearby sample appears unbiased relative to that of the 
underlying SN la population (see §7.3) while its Ay dis- 
tribution is clearly biased. For SALT-II, Xi and c describe 
the light-curve shape and color for each SN la (see § 5.2), 
2;i.ref = 2.6 and Cref — —0.26 correspond to the brightest 
SN la, and a, (3 are global SN la parameters determined 
in the cosmology fit (see § 8.2). 

Using the magnitude dimming parameter, we model 
the spectroscopic inefficiency due to apparent magnitude 
as 

.M _ _r -Mdnjm{z)] , (11) 



spec 



exp[ 



where m(z) is an exponential slope function determined 
by numerically solving 



dMdim Nsiuiz, Mdim) exp[- M dim /m{z)] 



DATA 



(^) 



(12) 



Here, iVsiM(^7 -A^dmi) is the number of simulated events 
in a two-dimensional bin of redshift and simulated A^dim, 
and N-dataIz) is the number of data events in the red- 
shift bin. In practice, m{z) is evaluated in discrete red- 
shift bins for z > Zcut and fit to an exponential function 
of redshift. 



m{z) = moexp[-(z - Zcnt)/mi] + ?7i2 



(13) 



where mo, mi, and m2 are determined in the fit, and 



m(z) 



oo for z < 



Since A^dim is defined for 
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y-band in the mlcs2k2 model and for B-band in the 
SALT-ii model, the associated mo, mi parameters are 
different. The parameter 7712 is non-zero only for the 
nearby sample. Qualitatively, Eq. 11 says that SNe that 
are intrinsically faint (large positive A) or extinguished 
(large Ay) will be under- represented at the high-redshift 
end of a sample. 

In principle, we now have two models for the spectro- 
scopic efficiency: one, denoted e^^^^, and explicitly given 
in Eq. 8, depends solely on redshift; the other depends 
explicitly on apparent brightness and implicitly on red- 
shift and is also constrained to match Eq. 8. While the 
latter model corresponds to expectations for Malmquist 
bias, the "redshift-only" model efp^,^ may be better suited 
to modeling spectroscopic selection that assigns lower 
priority to SN candidates near the cores of host galax- 
ies, since SN-host angular separation tends to decrease 
with redshift. Since both models satisfy Eq. 8, these two 
efficiency models can be linearly combined into a more 
general model: 



Nearby SN la DATA 



espoc(z, A^dim) = (1 - ^7M)espcc + ^^M^] 



M 



(14) 



To break this model degeneracy and determine the co- 
efficient Aj^, we compare the mean fitted extinction 
Ay(z) versus redshift for the data and the simulation 
after applying all efficiency factors, and we minimize the 
for the data-simulation difference. For the SDSS-II, 
Am ~ 0.8 ± 0.2, indicating that most of the search inef- 
ficiency is related to SN la apparent brightness. For the 
ESSENCE survey. Am ± 0.2, indicating that most 
of the search inefficiency is simply a function of redshift. 
For the SNLS sample, the best-fit Am = 1; however, 
since the data-simulation for Am = is only 0.3 
larger than the minimum -^M = 1) there is es- 

sentially no information on Am- The large uncertainty 
on Am is in part due to the relatively small spectro- 
scopic inefficiency for the SNLS. For the nearby sample, 
there is no redshift range for which the sample is 100% 
efficient, and we therefore need an additional constraint 
to determine the efficiency parametrization. Noting that 
the mean fitted extinction drops rapidly with redshift for 
the nearby sample (see Fig 11), we assume that Am = 1- 
Fig. 11 shows that the simulated efficiency works well in 
reproducing the strong extinction gradient in the nearby 
sample. 

As an illustration, the inferred efficiencies for the 
SDSS-II SN sample are shown in Fig. 12 as a function 
of Av, for different values of the redshift and shape- 
luminosity parameter A. The high quality of the simu- 
lation for the SDSS-II, ESSENCE, SNLS, and HST sam- 
ples is illustrated in Figs. 13-14, where we compare the 
observed redshift and flux distributions to simulations 
that include the efficiency functions derived above. The 
redshift comparisons have excellent /dof as a result of 
how the spectroscopic efficiency is determined. For the 
flux comparisons, the x^ /dof vary from 1 to a few; the 
larger x^ values are due to a few notable discrepancies 
in some of the flux bins. 

We have also compared the distributions of the number 
of epochs, earliest and latest epochs, noise, and peak 
colors, and find good data-simulation agreement in these 
distributions as well. 
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Fig. 11. — Left panel: mlcs2k2 fitted extinction Av versus red- 
shift for the nearby SN la sample (Jha et al. 2007). Right panel: 
fitted Ay distribution for the simulation using the parametrized ef- 
ficiency function of Eqs. 11 and 14 with Am = 1. Light-curve fits 
wore made using the mlcs2k2 prior of Eq. 5 and the exponential 
Ay distribution of §7.2. 
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Fig. 12. — Estimated efficiency versus extinction for the SDSS- 
II SN sample. Each panel corresponds to a different combination 
of redshift (z = 0.1, 0.3) and intrinsic SN la brightness (A = 
—0.3,-1-0.5 bright, faint). The curves correspond to different 
stages of the ofHciency: after the image-subtraction pipeline, egubtr 
(dotted); after spectroscopic confirmation, Egearch = ^subtr x ^spoc 
(dashed); and after selection cuts, Csurvey = ^search ^ ^cuts (solid). 
For z = 0.1, the spectroscopic efficiency is 100%, so the search 
(dashed) and image-subtraction (dotted) curves are the same. In 
all cases, ecuts ~ 0.8 at low extinction; this loss is mainly due to the 
requirement of good light-curve coverage, i.e., SNe that peak very 
early or very late in the observing season do not have adequately 
sampled light curves to satisfy the selection criteria. 

7. ESTIMATION OF HOST-GALAXY DUST PROPERTIES 

The Monte Carlo simulations used in the previous sec- 
tion to model survey selection functions rely on knowl- 
edge of host-galaxy dust properties, in particular on the 
underlying distribution of extinction, P(Ay), and on the 
mean reddening law parameter, i?y . Moreover, MLCS2k2 
distance estimates rely directly on our knowledge of these 
two quantities, which are assumed to be independent 
of redshift. In this section, we describe how we deter- 
mine these dust properties from the SDSS-II SN sample. 
These global dust properties are used in the mlcs2k2 
fitting prior for all SN samples. 

7.1. SDSS-II "Dust" Sample 
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Fig. 13. — Comparison of redshift distributions for data (dots) 
and simulations (histograms), for the SDSS-II, ESSENCE, and 
SNLS samples after applying the selection requirements in § 4. The 
simulations include all known effects, including the spectroscopic 
inefficiency. Each simulated distribution is scaled such that the 
total number of SNe matches the data. The SDSS data-simulation 
discrepancy for z > 0.4 is an artifact of our simple modeling of 
Espcc (Eq. 13). 

The results of § 6.2.1 indicate that current SN sam- 
ples suffer from significant spectroscopic selection effects 
(see Fig. 9). Since spectroscopic selection is likely to 
be biased against highly extinguished SNe, use of purely 
spectroscopic SN samples may lead to biased estimates 
of the distribution of host-galaxy dust properties and 
thereby to potentially biased distance estimates when a 
dust-distribution prior is applied. To address this issue, 
we use a nearly complete set of SDSS-II SNe la to mea- 
sure dust properties. For SDSS-II events with SN la-likc 
light curves that were not spcctroscopically confirmed as 
SNe la, we have embarked upon a program to obtain 
host-galaxy spectra and measure spectroscopic redshifts; 
we call these photometric SN la candidates. Based on 
distributions of the mlcs2k2 fit parameters as well as 
visual inspection of the light curve fits, we find that the 
requirement of a good SN la light-curve fit (see cut 5 in 
§ 4) to a well-sampled light curve is a good substitute 
in identifying SNe la when a confirming SN spectrum is 
lacking. 

To identify photometric SN la candidates, we start 
with all 4100 candidate events that were detected by the 
on-mountain frame-subtraction pipeline on two or more 
epochs in the 2005 observing season (see §2). We pro- 
cess these candidate light curves through Scene Model 
Photometry (§ 3.2), fit them with the MLCS2k2 method, 
and prioritize them for host-galaxy spectroscopy based 
on the quality of the light curve and the fit. We have 
obtained host-galaxy redshifts for the majority of candi- 
dates for which the host-galaxy r-band magnitude sat- 
isfies r < 20 and for a large subsample of fainter hosts 
as well. Adding this "host-z" photometric SN la sample 
to the spcctroscopically confirmed and probable SN la 
sample, the combined sample appears to be nearly spcc- 
troscopically complete to z ~ 0.3, when compared with 
the simulated sample. For z < 0.3, after the selection 
cuts of §4 are applied, the combined sample comprises 
81 confirmed SNe la, 6 probable SNe la, and 73 host-z 
SNe la. We refer to these 160 SNe la as the SDSS-11 
dust sample. To illustrate the importance of including 
the host-z subset, we note that the average fitted extinc- 
tion {Ay) is about 0.2 for the spcctroscopically confirmed 
sample and almost 0.4 for the host-z sample: ignoring the 
host-z subset would clearly lead to biased results for the 
distribution of host-galaxy dust properties. On the other 
hand, giving highest priority to bright host galaxies for 
the host-z follow-up program may preferentially select 
hosts of more extinguished SNe. 

To illustrate our understanding of the efficiency. Fig. 15 



shows the redshift distribution for the dust sample, com- 
pared with the simulation that includes known losses 
from the image-subtraction pipeline and selection cuts 
but docs not include losses related to spectroscopic ob- 
servations. This comparison shows that the dust sample 
is indeed nearly complete for redshifts z < 0.3. We can 
obtain an independent estimate of the dust-sample com- 
pleteness by counting the number of photometric SN la 
candidates that pass our selection cuts, that have a pho- 
tometric redshift (based on the host galaxy or the SN 
light curve) less than Zphot = 0.3, and that have not 
yet been targeted for host-galaxy spectroscopy. There 
are nearly 40 such events, indicating a completeness, af- 
ter selection cuts, of about 80% for the dust sample at 
z < 0.3. 

7.2. Determination of Dust Reddening Parameter Ry 

Although Rv varies along sight-lines through the 
Milky Way and likely varies with local environment 
within galaxies, we follow standard practice in treating 
it as a global parameter, because most current SN data 
are not adequate to determine it on an object-by-object 
basis. After briefly reviewing previous determinations of 
Rv, we describe the method we have developed to deter- 
mine Rv and its results. 

Measurements of dust properties in the Milky Way 
have favored an average value of Rv ~ 3.1 (Fitzpatrick 
& Massa 2007) along sight lines with moderate to sub- 
stantial extinction, and this has been used as a canon- 
ical value in the literature. However, studies of stellar 
colors in the direction of several thin cirrus clouds in the 
Milky Way indicate a value of Rv ^ 2 for those relatively 
low-extinction environments (Szomoru & Guhathakurta 
1999). For galaxies that host SNe la, two methods have 
been commonly used to determine Rv- The first in- 
fers an average or global Rv value based on statisti- 
cal averages of optical light-curve colors. This method 
has been applied to the nearby SN sample (Phillips 
et al. 1999; Altavilla et al. 2004; Reindl et al. 2005; 
Riess et al. 1996; Nobih & Goobar 2008), resulting in 
Rv values in the range of 2 — 3, somewhat lower than 
the Milky Way average value. The second method uses 
SNe la that have densely sampled, high signal-to-noise 
optical and NIR photometry, for which Rv can be esti- 
mated for individual events: this method has resulted in 
~ 1-2 (Krisciunas et al. 2007, 2006; Ehas-Rosaet al. 
2006; Wang et al. 2008), significantly below the canoni- 
cal Milky Way value. Calculations and simulations have 
shown that multiple scattering by circumstcUar dust can 
lead to lower values for the reddening parameter inferred 
from SNe la, Rv ~ 1.5-2.5 (Wang 2005; Goobar 2008). 

Within the framework of the MLCS2k2 light-curve 
model, we have used the SDSS-II dust sample and a 
variant of the average color method to make a new de- 
termination of Rv- Previous measurements with this 
method were based on samples with large and unknown 
selection effects; our determination of Rv is based upon 
a SN la sample with a well-understood selection effi- 
ciency. We assume that reddening is due to dust ex- 
tinction with a wavelength dependence described by the 
CCM89 model and that the mlcs2k2 model parameters 
{M,p,q in Eq. 1) accurately describe the SN la bright- 
ness and colors. The measurement of Rv is based upon 
comparing the average colors as a function of SN epoch 
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Fig. 14. — Comparison of fiux distributions for data (dots) and simulations (histograms) for SDSS-II, ESSENCE, SNLS, and HST. The 
observer-frame passband is indicated on each pfot (HST filter symbols are =NIC-F110W, 1 =NIC-F160W, 2 =ACS-F606W, 4 =ACS- 
F775W, 6 =ACS-F850LP). Each simulated distribution is scaled to have the same total number of entries as the data; bins with zero entries 
arc not plotted. The calibrated flux is io~'''^'™''g+ll. As in Fig. 13, the simulations include all the survey inefficiencies. 



^60 



iS 50 



40 
30 
20 

10 




• SDSS data 

(confirmed + host-z) ^ 

- Sim 



m 



simultaneously in SDSS-II, the colors are determined di- 
rectly from the data without the need to interpolate in 
time. 

We compare these color-versus-epoch measurements 
with those of a grid of simulated samples, where each 
sample is generated with a different value of Ry and 
Ay , and where Ay describes the generated Ay distribu- 
tion, with PiAy) = eyip{- Ay /Ay). The Ry and Ay 
grid sizes are 0.2 and 0.05, respectively. The simulated 
and data samples are subject to the same selection cuts 
and fit with mlcs2k2 in the same way. We determine the 
best-fit values of Ry and Ay by minimizing the following 
statistic between the data and the grid of simulations, 
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Fig. 15. — Redshift distribution for the SDSS-II dust sample 
(dots), which includes confirmed, probable, and host-2 photomet- 
ric SN la events, and for the simulation (histogram), which in- 
cludes losses from the image-subtraction pipeline and selection 
cuts. Shaded region reflects simulated uncertainty in the SN la 
rate as explained in Fig. 9 caption. Vertical arrow indicates z < 0.3 
requirement to select the dust sample used for the determination 
of host-galaxy dust properties. 



of the SDSS-II data to those of the simulation. By us- 
ing the nearly complete z < 0.3 dust sample, we mini- 
mize potential selection bias against extinguished SNe la; 
this sample also has a selection efficiency that is well de- 
scribed by the simulation. 

For the dust sample, we compute three mean observed 
SN la colors, g — r, r — i, and g — i, in one-day bins in 
rest-frame epoch, as shown in Fig. 16. Although the third 
color (g — i) is redundant in most cases, it provides infor- 
mation in the few cases in which the r-band measurement 
is not available. The rest-frame time is relative to the 
time of peak brightness in i?-band as determined by the 
mlcs2k2 light-curve fit (see § 5.1). While the light-curve 
fits include all observations regardless of signal-to-noise, 
the estimates of the mean observed colors include only 
those observations for which the signal-to-noise ratio is 
greater than 4. Since the r, i measurements are taken 
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where the epoch- index e runs over 1-day bins with — 5 < 
T^cst < 30 days, the data averages O^"*"^ are taken over 
all dust-sample SNe and epochs e surviving the cuts 
above, and the measured colors in the simulated sam- 
ples depend upon the input values of Ry and Ay. Each 
color uncertainty, e.g., CT(g_r)^, is estimated as rms/VlVe, 
where is the number of color measurements (summed 
over all SNe) at epoch e. The second term in brackets 
is a weighting to account for the mlcs2k2 model un- 
certainty; is the minimum model uncertainty at 

the epoch of peak brightness, and model 
uncertainty at epoch e. The model- uncertainty ratio is 
therefore unity at Trcst = and decreases as the model 
uncertainty increases for epochs away from peak bright- 
ness, so that epochs with large errors are downweighted. 
We have tested this method with 100 simulated mock 
data samples as described in § 6; the input values of Ry 
and Ay are recovered, and the statistical uncertainties, 
although they are smaller than the grid sizes, match the 
spread in recovered values. 
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Using the SDSS-II data and simulations, wc find 

i?v = 2.18±0.14,tat±0.48,y,t (16) 

Av = 0.358 ± 0.026,tat ± 0.068syst (17) 

and a correlation coefficient of 0.17. The relatively small 
correlation coefficient confirms that the method is inde- 
pendently sensitive to Ry and Ay. Figure 16 compares 
the average observed colors with those of the simulation 
using the best-fit values in Eq. 16. The values are 
somewhat larger than expected, particularly for g ~ r. 
There is also a notable data-simulation discrepancy for 
epochs past about 10 days for the g — r and g — i colors, 
although these late-time epochs carry less weight in the 
due to the increasing model errors. 

To estimate the systematic uncertainties in this mea- 
surement, we have varied aspects of the procedure and 
determined their effects on the recovered Ry and Ay. 
In particular, we lowered the redshift cutoff for the dust 
sample to z < 0.25 (nominal is 0.3), varied the simulated 
redshift distribution of the underlying SN la population, 
dN/dz = a(l 4- z)^, within the correlated Icr errors on 
a and (3 from Dilday et al. (2008), varied the minimum 
epoch from —10 to days (nominal is —5 days), varied 
the maximum epoch from 4-25 to 4-35 days (nominal is 
4-30 days), varied the minimum signal-to- noise between 
2 and 8 (nominal is 4), excluded each of the three colors 
g — r, g — i, and r — i individually from the minimiza- 
tion, and ignored the MLCS2k2 model-uncertainty terms 
in the statistic. We also ran the simulation without 
intrinsic color fluctuations (see §6.1). The changes due 
to each of these variations to the nominal results for Ry 
and Ay were added in quadrature to obtain an estimate 
of the total systematic uncertainty. Table 3 summarizes 
the contributions to the systematic uncertainties. The 
largest source of uncertainty comes from reducing the 
redshift-range from 0.3 to 0.25, which changes Ry by 
0.36 ±0.17, and the error reflects the uncorrelated uncer- 
tainty. Although this .Ry-shift is marginally consistent 
with a random fluctuation, we have included the shift as 
a systematic error. 

As a crosscheck of how the inferred Ry depends on the 
assumed exponential form of the extinction distribution, 
P(Ay), we have repeated the procedure using different 
distributions for Ay in the simulation. In particular, 
we use a flat, truncated Ay distribution, P{Ay) = 1 
for Ay < 2 Ay, and P{Ay) = when Ay > 2 Ay, this 
results in a negligible change in the inferred values of Ry 
and Ay. Thus, the inferred value of Ry appears to be 
insensitive to the assumed Ay distribution. 

7.2.1. The SALT-ll Approach to Measuring Ry 

The SALT-ii analog of Ry (called f3, see Eq. 29 below) 
is determined by a very different method that involves 
minimizing the residual scatter in the Hubble diagram. 
For comparison, we have tried a similar procedure within 
the MLCS2k2 framework, selecting the global value of Ry 
that minimizes the scatter in the Hubble diagram. For 
this test, we select spectroscopically confirmed SNe la 
from the SDSS-II sample that satisfy the light-curve cri- 
teria of § 4 and that have redshifts z < 0.15; the resulting 
sample includes 24 SNe la. We concentrate on this low- 
redshift sample because it is nearly complete and because 
the measurement uncertainty on each distance modulus 
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Fig. 16. — Mean SN la colors (g — r), {g — i), and (r — i) vs. rest- 
frame epoch, for the SDSS-II SN dust sample. Data are shown by 
the filled circles; the histogram overlay shows the simulation using 
the best- fit values of Ry and Ay- The data-simulation (for 36 
degrees of freedom) is indicated on each panel. 

estimate is well below the intrinsic scatter. The resulting 
Hubble scatter- minimized Ry value is about 1.7, roughly 
0.5 below our nominal result in Eq. 16. This result holds 
whether we use the default prior (Eq. 5) or a flat prior 
in the light-curve fits. When we apply this procedure 
to the larger SDSS-II dust sample (§ 7.1) that extends 
to z = 0.3, the Hubble scatter-minimized Ry value is 
slightly smaller. Hicken et al. (2009b) have applied the 
same approach to the nearby CfA3 sample (Hicken et al. 
2009a) and find that Ry = 1.7 minimizes the Hubble 
scatter, consistent with our results for the SDSS-II sam- 
ples. 

To test this approach to extracting Ry , we have imple- 
mented it on the simulated SDSS-II SN sample (§ 6). We 
generate a set of light curves with fixed Ry and fit them 
with mlcs2k2 in the same way that the data are fit. For 
the default color-smearing model of intrinsic SN luminos- 
ity variation (see § 6.1), the Hubble scatter-minimized 
Ry extracted from this process is biased low by 0.5 with 

TABLE 3 

Summary of uncertainties for the determination of Ry 

AND Ay. 



source 


criRv) 


^{Ay) 


statistical 


0.14 


0.026 


redshift range 0.25 — 0.3 


0.36 


0.031 


assumed dN/dz : /3 = 0.9 - 2.1 


0.11 


0.004 


epoch range 


0.25 


0.045 


signal to noise > 2, 8 


0.11 


0.011 


exclude one color 


0.08 


0.030 


ignore mlcs2k2 model uncertainty 


0.05 


0.018 


remove color smearing model 


0.14 


0.017 


'I'otal Systematic 


0.48 


0.068 


Total Uncertainty 


0.50 


0.072 
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respect to the input value. This bias is consistent with 
the difference we see between the seatter-mininiized Ry 
and the Ry we infer from the mean colors of the SDSS- 
II dust sample. For the alternative "coherent luminosity 
smearing" model, however, which results in no intrinsic 
color variations, the scatter-minimized Ry is unbiased. 

While the result above is suggestive, an important 
caveat is that we have not evaluated the uncorrelated sys- 
tematic uncertainty on the difference between our nom- 
inal Rv extraction and the scatter-minimized Ry- Our 
study suggests, however, that the SALT-ii (3 parameter 
could be biased low if the color-smearing model is a rea- 
sonable description of intrinsic SN la luminosity varia- 
tion. 

7.3. Determination of the Underlying Distributions of 
Ay and A 

The mlcs2k2 simulation and the prior used in the 
MLCS2k2 fitting method require knowledge of the un- 
derlying distribution for the l/-band extinction due to 
host-galaxy dust, P{Av). The distribution is also needed 
for the shape-luminosity parameter A, although the lat- 
ter has less impact on the results, because A is better- 
determined by the light-curve data for each SN. We de- 
termine these distributions using the Bayesian unfolding 
method of D'Agostini (1995), which essentially uses un- 
derlying trial distributions in the simulation to make pre- 
dictions for the observed distributions of fitted Ay and 
A. To limit selection biases in this procedure, we use the 
SDSS-II dust sample of § 7.1. 

We fit the SDSS-II light curves with mlcs2k2 using a 
flat prior on Ay, i.e., the fitted Ay value is allowed to be 
negative. In the fit, the extinction as a function of wave- 
length is given by the CCM89 model with Ry = 2.18, as 
derived in § 7.2. The underlying distributions for Ay and 
A are determined such that when they are input into the 
simulation, the fitted distributions from the simulated 
light curves match the distributions from the mlcs2k2 
fits to the SDSS-II data. Technical descriptions of this 
procedure and of the assumed underlying distributions 
are given in Appendix D. 




Fig. 17. — Fitted Ay (left) and A (right) distributions, where 
the fits use the mlcs2k2 model with a flat Ay prior. The mea- 
sured SDSS-II dust sample distributions are shown by dots; the 
simulations are shown by histograms. The underlying (generated) 
Ay distribution is shown by the dashed curve in the left panel. 

The results for the data and Monte Carlo simulation 
are shown in Fig. 17. Although the generated Ay dis- 
tribution includes only non-negative values of Ay, there 
is a tail of fitted negative values that is well-described 
by the simulation and arises from photometric errors 



and intrinsic SN color variations. Our procedure de- 
termines the underlying Ay distribution without assum- 
ing a functional form for P(Ay), but it turns out that 
this distribution is well described by an exponential, 
P{Av) = exp(-Ay/Tv), with 



TV = [0.334 ± 0.072] 
= 0.334 ±0.088 . 



0.10 X {Ry - 2.18) 



(18) 
(19) 



The uncertainty in ry includes statistical and system- 
atic uncertainties as described in Appendix D. The /in- 
dependence is shown explicitly in Eq. 18, along with the 
internal measurement error. Eq. 19 shows the total un- 
certainty including the uncertainty on Ry. The agree- 
ment between ry and Ay (from § 7.2) is well within the 
expected dispersion based on simulated tests. 
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Fig. 18. — (a) Underlying Ay distributions, P{Ay): dotted is 
from JR,K07, solid is the glos distribution from WV07, and the 
hatched region shows the distribution used in this paper, derived 
from the SDSS-II SN dust sample; (b) Underlying A distribution, 
P(A), used in the mlcs2k2 prior, also derived from the SDSS-II 
dust sample. 



The underlying Ay distribution, including the uncer- 
tainties, is shown as the hatched region in Fig. 18a; we 
use this P{Av) as part of the MLCS2k2 prior in Eq. 
5. Our inferred Ay distribution is marginally consistent 
with that of JRK07, who derived a prior of exponential 
form from the nearby sample, with = 0.46 (dashed 

curve in Fig. 18a). Our Ay distribution is also some- 
what consistent with the "galactic line-of-sight" (glos) 
prior, based on theoretical considerations (Hatano et al. 
1998; Commins 2004; Riello & Patat 2005), that was 
used by WV07 for the ESSENCE analysis (solid curve 
in Fig. 18a). 

The underlying distribution of A is described by an 
asymmetric Gaussian with mean Aq = —0.24, and Gaus- 
sian widths of o-_ = 0.24 for A < and a+ = 0.48 for 
A > (see Appendix D for uncertainties). As shown in 
the right panel of Fig. 17, when input into the simulation 
this distribution leads to a distribution of fitted A that 
is in good agreement with the observed distribution of 
fitted A from the SDSS-II dust sample. The underlying 
A distribution, P{A), is shown in Fig. 18b. In addition 
to the asymmetric Gaussian, we have added a tail to 
the A distribution at positive A, allowing for underlu- 
minous SNe where the underlying distribution is poorly 
measured due to small-number statistics; when used in 
the mlcs2k2 prior of Eq. 5, this tail ensures that the 
fitter is not heavily biased against underluminous SNe. 
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Since A is well determined from the light-curve shape, 
the fitted value of A has little sensitivity to the func- 
tional form of -P(A). To check the effect of the arbitrary 
tail in P(A), we have run the fits with amplitude of the 
tail region multiplied by half and by two; in both cases 
the rms- variation in fitted A is 0.01. 

To check that the derived extinction distribution re- 
flects that of the global SN population rather than that 
of a (possibly biased) SDSS-II sub-sample, we compare 
the fitted Ay distributions for the data and for simu- 
lations generated with the underlying Ay and A distri- 
butions derived above for each SN sample in Fig. 19. 
Here, the fits are carried out with a flat Ay prior, al- 
lowing Ay < 0. The overall agreement is good for all 
samples. This test illustrates that the inferred negative 
extinction values (when an uninformative prior is used) 
are consistent with being artifacts of the combination of 
low signal-to-noise and intrinsic color fluctuations, since 
the simulation is generated with Ay > 0. Note that the 
procedure used to determine the underlying Ay distri- 
bution ensures that the data and simulation agree well 
for positive Ay for the SDSS-II sample, but parameters 
have not been adjusted to match the distribution for the 
other samples, or for negative Ay. 




0.5 

fit Ay (no prior) 

Fig. 19. — Comparisons of fitted Ay distributions, using a flat 
prior, for data and simulations. Simulations were generated us- 
ing the underlying Ay distribution derived from the SDSS-II dust 
sample, with Ay > 0. SN samples used in the cosmology analysis 
are indicated on each panel. 



8. FITTING THE HUBBLE DIAGRAM FOR 
COSMOLOGICAL PARAMETERS 

The estimation of cosmological parameters from super- 
novae is based on measurements of the luminosity dis- 
tance, di, as a function of redshift. For a Friedmann- 
Robertson- Walker cosmological model, assuming the 
stress-energy comprises non-relativistic matter (M) and 
dark energy (DE) with constant equation of state param- 
eter w =p/pc^, the luminosity distance depends on four 
parameters: w, the matter density f^M, the dark energy 



density IIde, said the Hubble parameter Hq'- 
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where the curvature density f2i; = 1 ~ ^^m ~ ^^de, and the 
function Sk{x) = sin(a;) for < 0, Sk{x) = sinh(a;) for 
fik > 0, and Sk{x) = x for a flat Universe with il^ = 0. 
The dimensionless expansion rate is given by 



E{z) = H{z)/Ho = 



1/2 



VlMil + zf + f7DE(l + z)3(i+-) + a(l + zf (21) 



For this analysis, we assume that the dark energy equa- 
tion of state parameter w does not evolve in time, mainly 
because current data do not yield precise constraints on 
the time derivative of w. For constraints on time- varying 
w, see Sollerman et al. (2009). Loosely, one can interpret 
the constraints we derive on constant w as constraints 
on the mean value of w{z) over the redshift interval 
< z < 0.4. More precisely, one can interpret these con- 
straints as providing information on one parameter in a 
two-parameter model that describes the evolution of w. 
That is, if the evolution of w is described by a linear two- 
parameter model, e.g., by w{z) = wo + Waz/{l + z), then 
there is a pivot redshift Zp at which the measurements of 
Wo and Wa are uncorrelated and the error in Wp = w{zp) 
reaches a minimum. For the combined SN la data sets 
considered in this paper, Zp « 0.25, and the constraints 
we derive on constant w are equivalent to constraints on 

Wp. 

We will study constraints on two lower-dimensional 
subspaces in the above family of cosmological models. 
For the first, which we denote ACDM, we consider mod- 
els in which the dark energy is vacuum energy, i.e., the 
cosmological constant A, with w = —1, but we allow 
non-zero spatial curvature the parameters of interest 
arc 51m and JIa. In the second case, denoted FwCDM, 
we assume spatial flatness, fik = 0, but allow w to dif- 
fer from —1; here the parameters of interest are Qm 
and w. The rationale for considering the ACDM model 
is that one can consider the cosmological constant as 
a "null hypothesis" for dark energy, so it is worth ex- 
ploring whether it provides a reasonable description of 
the data. The rationale for the FwCDM model is that 
the WMAP data from the cosmic microwave background 
(CMB) anisotropy constrain the spatial curvature to be 
very small. 

In deriving parameter estimates, we will combine the 
SN data with two independently measured constraints. 
The first is from the measurement of baryon acoustic 
oscillation (BAO) feature in the SDSS Luminous Red 
Galaxy (LRG) sample by Eisenstein et al. (2005). The 
BAO measurements constrain several different parame- 
ters (Eisenstein et al. 2005; Percival et al. 2007), depend- 
ing on whether and how information from the CMB is 
used; wc explore this in more detail in Sollerman et al. 
(2009) and Lampeitl et al. (2009). Here, we follow Astier 
et al. (2006) and WV07 in using the combination of angu- 
lar diameter distance, Hubble parameter, and I^m given 
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by Eisenstein et al. (2005), 

A{zi ; til, fiM,f^DE) = 
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where the SDSS LRG BAO constraint is given by 
ii;,l]M,f^DE)-0.469)/0.017]2 , 
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with the effective LRG rcdshift zi ~ 0.35. The second 
constraint is from the five-year results of the Wilkinson 
Microwave Anistropy Probe (Komatsu et al. 2008). We 
use the shift parameter, 



^(^cmb; WjI^Mj^^De) = 



with the constraint 
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,(24) 



XCMB — [(i?(zcMB; «i,r!M,ODE)-1.710)/0.019]2, (25) 

where zcmb is the rcdshift of decoupling of the CMB. 
Although zcMB depends upon J^m and upon the baryon 
density ^I-q at the ~ 2 percent level, we fix this rcdshift to 
the WMAP5 maximum-likelihood value of zcmb = 1090 
(Komatsu et al. 2008). 

Minimizing Xbao + Xcmb without using SNe, the 
best-fit cosmological parameters for the FwCDM model 
are w = -0.80 ± 0.20 and Q.m = 0.30 ± 0.05;'*i for the 
ACDM model, VIm = 0.27 ± 0.02 and f^A = 0.74 ± 0.02. 
Combining these measurements with the SN samples re- 
sults in improved constraints for the FwCDM model but 
has little impact for the ACDM model. 

We have analyzed several different combinations of the 
five SN la data sets mentioned at the beginning of §5: 
SDSS-II, nearby (low-redshift), ESSENCE, SNLS, and 
HST. In § 9 and § 10, we present systematic uncertainties 
and results for the six sample combinations {a-f) shown in 
Table 4. Combination (a) includes only the SDSS-II SN 



TABLE 4 

Sample-combinations used to extract cosmological 
parameters. 

sample-combination 
SDSS-only 
lb) SDSS-I-ESSENCE+SNLS 

(c) Nearby-I-SDSS 

(d) Nearby-|-SDSS-I~ESSENCE+SNLS 

(e) Nearby-l-SDSS-l-ESSENCE+SNLS-fHST 
(/) Nearby-I-ESSENCE+SNLS 



data, without lower- or higher-redshift SNe. For combi- 
nation (b), the nearby SNe la are again excluded, leaving 
the SDSS-II SN sample to serve as the low-redshift an- 
chor for the Hubble diagram. Combination (f) excludes 
the SDSS-II sample, so that we can directly compare with 
previously published results such as those of WV07. In 

'^^ If one marginalizes over f^M, the marginalized value of to is 
0.09 smaller. When SN constraints are included, the difference be- 
tween best-fit and marginalized parameter values is much smaller. 



combination (c), the SDSS-II serves as the 'high-redshift' 
sample. Combination (d) includes all four ground-based 
samples; combination (e) includes all five samples and is 
used for our nominal analysis. 

8.1. Fitting Distances with MLCS2k2 

As described in §5.1, mlcs2k2 provides an estimate of 
the distance modulus, ^, for each supernova. Cosmolog- 
ical parameter estimates are derived by minimizing the 
following statistic (= —2 In of the posterior probabil- 
ity) for the SN la sample over a grid of model parameter 
values, 



E 
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"XbAO + XcMB 
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where is the distance modulus estimated from the 
mlcs2k2 fit for the i'th supernova, Zi is its spectroscop- 
ically determined redshift, and //(z^; w, r^M, ^^de, -f^o) = 
5 log((iL/10 pc) is computed from Eq. 20. The Xbao ^'^^ 
XcMB terms in Eq. 26 incorporate the prior information 
from the SDSS LRG BAO measurement (Eq. 23) and 
the WMAP CMB measurement (Eq. 25). When report- 
ing values and errors on individual cosmological parame- 
ters, we use the values at the x^-minimum, marginalized 
only over iJo (due to the degeneracy between the Hubble 
parameter and the fiducial peak rest-frame model mag- 
nitude). In determining w, marginalizing over shifts 
the value from the X^-minimum by 0.025 for the SDSS- 
only sample and by ^ 0.01 for the other sample combi- 
nations. 

In Eq. 26, the distance-modulus uncertainty is given 

by 



(27) 



where crJJ' is the statistical uncertainty reported by 
mlcs2k2, cr™* = 0.16 is the additional (intrinsic) error 

added so that the x^ per degree of freedom is unity for 
the Hubble diagram constructed from the nearby SN la 
sample (§ 10.1) and cr^j is calculated from the rcdshift 
uncertainty as described below. 

The redshift uncertainty contains two components: 
from spectroscopic measurements and from peculiar mo- 
tions of the host galaxy. The first source of uncertainty, 
c^^spec, is from the uncertain redshift determination, ei- 
ther from the spectrum of the SN or from its host galaxy. 
For SDSS-II SNe we use (iz^spec = 0.0005 for host-galaxy- 
based redshifts, a^.s^ec = 0.005 for SN-based redshifts, 
and the reported SDSS redshift errors for host redshifts 
from the SDSS spectroscopic galaxy survey (§ 3.1). For 
the ESSENCE and SNLS data, we use the estimated red- 
shift errors reported in their public data tables. For the 
nearby sample, redshift measurement errors were usu- 
ally not reported; we take 50 km/sec as a conservative 
estimate (M. Hamuy, private communication). The sec- 
ond source of redshift uncertainty, (Tz^pec, arises from pe- 
culiar velocities of and within host galaxies. We take 
o'z,pec = 0.0012, the quadrature sum of typical galaxy 
peculiar velocities of 300 km/sec and typical internal mo- 
tions of 200 km/sec. For most of the SDSS-II SNe, for 
which the redshift is obtained from the host galaxy, the 
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contribution from internal motions is overestimated be- 
cause these spectra are averaged over the internal galaxy 
motions. The final rcdshift uncertainty is defined to be 
cr^ ~ cr^.spec + '^z.pec- To simplify the treatment of rcd- 
shift errors, we project these uncertainties onto distance 
modulus using the expression for the distance-redshift 
relation for an empty universe, 

Using a different cosmological model to compute cr^ from 
az leads to negligible changes in the Hubble diagram fits. 

Galaxy peculiar velocities are not random, as the above 
treatment assumes, but are spatially correlated, since 
they are induced by the gravitational effects of large-scale 
structure. SN observations are affected by the peculiar 
velocities of both the host galaxies and of the Milky Way. 
We have accounted for the Milky Way peculiar velocity 
by correcting for the CMB dipole: all supernova red- 
shifts have been transformed into the comoving frame of 
the CMB. This is particularly important for the SDSS- 
II, because the equatorial stripe comes within 7° of the 
CMB dipole direction, and the redshift correction from 
the heliocentric frame is negative along the entire stripe. 
We use (1+^;) = {l+ZQ)/{l—Vo-n), where z© is the helio- 
centric redshift as described in § 3.1, n is the unit vector 
pointing from earth to the SN, and vq is the CMB veloc- 
ity vector, with a magnitude of 371 km/sec and direction 
given by Galactic coordinates I = 264.14°, b = -1-48.26° 
(Fixen et al. 1996). The CMB-frame redshifts for the 
non-SDSS-II samples are taken from the literature. 

In transforming to the CMB frame and making no 
other corrections for velocity correlations, we are im- 
plicitly assuming that the peculiar velocities of the host 
galaxies are uncorrelated with that of the Milky Way 
and approximately uncorrelated with each other. That 
assumption is a good approximation for the SDSS-II 
and higher-redshift SN samples, which cover large spa- 
tial volumes and are distant from the Milky Way. It is 
not a good approximation for the nearby SN sample on 
both theoretical (Hui & Greene 2006; Cooray & Cald- 
well 2006) and observational (Neill et al. 2007) grounds. 
Not including velocity correlations means that the low- 
redshift supernovae are over-weighted in the xfi statistic 
and that SN-derived cosmological parameter errors in the 
literature have been underestimated. However, inclusion 
of the SDSS-II SN sample significantly reduces the im- 
pact of the uncertainties due to low-redshift peculiar ve- 
locities. 

8.2. Fitting Distances with SALT— II 

In the first stage of the SALT-ii analysis framework, 
the photometry for each supernova light curve is fit 
to an empirical model (§ 5.2) to determine a shape- 
luminosity parameter (xi), a color parameter (c), and 
an overall flux normalization {xq). For the SALT-II Hub- 
ble diagram analysis, the reference i3-band magnitude 
is m*B = -2.51og[a;o/ dA'Mo(t, A')T^(A')] ■ The flux- 
integral is the same as Eq. 7 with z = c = xi = and 
f = B. The fitted parameter xq depends on the lumi- 
nosity distance and the SN brightness. 

To estimate cosmological parameters, the above pa- 
rameters are related to the distance modulus for each 



supernova by the expression 

Hi = mg,^ — Al + a ■ xi^i — P ■ Ci . (29) 

The global parameters M, a, and /? describe the SN la 
population and are estimated simultaneously with the 
cosmological parameters by carrying out a minimiza- 
tion using an expression analogous to Eq. 26. The min- 
imization and error estimation are performed using the 
program MINUIT.^^ The expression for the distance mod- 
ulus uncertainty, tr^, is similar to that in Eq. 27 for 
mlcs2k2, and the redshift uncertainty is treated iden- 
tically. The intrinsic dispersion (crjf*) is determined sep- 
arately for each sample combination by setting the global 
best-fit for that sample combination to unity, in con- 
trast to the mlcs2k2 method for which we determine 
crj,"* solely from the nearby sample. In addition, the 

SALT-II expression for xfi includes a covariance matrix 
to account for correlations between the parameters xi, 
c, and Xq (§ 9.2). In determining cosmological parame- 
ters, the Hubble parameter is marginalized over, but the 
parameters a and /? are not. 

If the color corrections were due only to extinction 
by host-galaxy dust, the SALT-ii parameter f3 would be 
equal to the Cardelli et al. (1989) extinction parameter, 
Rb (Conley et al. 2007). Further, if host-galaxy extinc- 
tion were similar to the mean extinction in the Milky 
Way, one would expect (3 ~ Rb = Rv + 1 — 4. However, 
in SALT-II, /3 is left as a free parameter that is determined 
in the global fit to the Hubble diagram. 

To account for selection bias, we have applied the 
SALT-II fitting method to simulated sample combinations 
a-/. The Monte Carlo simulations are generated from the 
SALT-II model using parameters based on our analysis of 
the data (§ 10.2). In particular, a ~ 0.13 and (5 = 2.56, 
and the simulated distributions of c and xi are described 
by Gaussians with c = 0.04, CTc = 0.13, xi = —0.13, 
(Jx^ = 1.24. These distributions were estimated directly 
from the light-curve fits rather than from simulating un- 
derlying c and xi distributions such that the observed 
and simulated distributions match. The intrinsic lumi- 
nosity is randomly varied using a coherent luminosity 
smearing factor of 0.12 mag, and the spectroscopic effi- 
ciency is based on the SALT-ii parameters described in 
§6.2.1. 

To isolate a potential SALT-ii bias from a shift due 
to the CMB and BAO priors, the simulated SN sample 
is generated with the cosmological parameter values ob- 
tained from fitting without any SNe, = —0.80 and 
Vlyi = 0.30. The bias in the cosmological parameters is 
estimated to be w—w^ and JIm— ^^M' where w and IIm are 
obtained from the analysis of the simulated SN sample 
combined with the BAO-I-CMB priors. To increase the 
significance of the measured bias, 500 data-sized samples 
were simulated and analyzed; the resulting statistical un- 
certainty on the w-bias is typically 0.004. The average 
bias for each sample combination and cosmological model 
is shown in Table 5. For the results shown below in Ta- 
bles 16 and 17, these bias corrections have been added 
to the cosmological parameters. 

9. SYSTEMATIC UNCERTAINTIES 
http: / /wwwasdoc .web . cern. ch/wwwasdoc/minuit/minmaiii.htinl 
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The likelihood analysis described in the previous sec- 
tion accounts for the impact of statistical errors on the 
determination of cosmological parameters. Here we give 
a detailed description of systematic uncertainties on the 
cosmological parameters w and Qm within the context 
of the FiuCDM model. For the ACDM model, the sys- 
tematic uncertainties are evaluated in a similar fashion 
and are summarized along with the results in § 10. For 
both mlcs2k2 and SALT-II and for each of the six sam- 
ple combinations {a-f) in Table 4, we have carried out 
several dozen systematic tests by varying analysis pa- 
rameters and methods. The resulting variations in the 
Hubble diagram and shifts in the parameter w (denoted 
Sw below) are used to assess the systematic uncertainty. 
To help gauge the significance of systematic shifts ob- 
served in the data, the same systematic parameter and 
method variations have also been applied to the analysis 
of Monte Carlo simulations generated with w ^ —I and 
r^M = 0.3, using the efficiencies determined in § 6. There 
are four main categories of systematic uncertainties: (i) 
uncertainties in SN la model parameters obtained from 
the training procedure for the light-curve fitter; (ii) un- 
certainties in reddening from host-galaxy dust and in in- 
trinsic SN color variations; (iii) errors in survey selection 
efficiencies and in associated selection biases; and (iv) 
uncertainties in relative photometric calibration between 
the nearby sample and the higher-redshift surveys. 

We discuss the impact of these uncertainties in the eon- 
text of the MLCS2k2 method in §9.1 and for the SALT-ii 
method in §9.2. Summaries of the systematic uncertain- 
ties for the six sample combinations listed in Table 4 are 
presented in Tables 6 and 7 for mlcs2k2 and in Tables 
8 and 9 for SALT-ii. For the flat FwCDM model, these 
tables show systematic errors in the dark energy equa- 
tion of state parameter w and matter density parameter 
J7m, including the BAO and CMB priors discussed in §8. 
A more detailed discussion of one major source of sys- 
tematic uncertainty, associated with data-model discrep- 
ancies in the rest-frame [/-band, is postponed to §10.1.3 
and 10.2.4. 

9.1. Systematic Uncertainties with mlcs2k2 
Rest frame [/-band 

As discussed in detail in §10.1.3, there are systematic 
discrepancies between the mlcs2k2 rest-frame [/-band 
model and the light-curve data for all but the nearby 
SN la sample. We have therefore carried out a test 
in which the observer-frame filter corresponding to rest- 
frame [/-band is excluded from the light-curve fits and 

TABLE 5 

SALT-II CORRECTIONS TO COSMOLOGICAL PARAMETERS BASED ON 
SIMULATIONS. 



Sample 


FtoCDM 


ACDM 


combination^ 


5w 




SUa 


SQm 


(a) 


-0.04 


-0.010 


0.001 


-0.001 


(6) 


-0.02 


-0.007 


-0.003 


0.004 




-0.02 


-0.005 


0.001 


-0.001 


id) 


0.00 


0.006 


-0.006 


0.007 


(e) 


0.00 


0.001 


-0.008 


0.005 


if) 


0.02 


0.006 


-0.004 


0.005 



^ a-f arc defined in Table 4. 



the resulting distance estimates. Figure 30 shows that for 
the SDSS-only sample, the exclusion of [/-band causes a 
mean systematic shift in estimated distance modulus of 
0.12 ± 0.02 mag for z > 0.21. The resulting tilt in the 
Hubble diagram leads to a systematic change of —0.310 in 
the equation of state parameter w for both the SDSS-only 
(a) and Nearby-|-SDSS (c) sample combinations; we in- 
clude this shift as an asymmetric systematic uncertainty, 
as indicated by the minus signs in the first row of Ta- 
bles 6-7. 

For the other sample combinations {b,d,e,]), which in- 
clude the ESSENCE & SNLS samples, the exclusion of 
[/-band results in w-shifts of 0.04 to 0.08. Based on tests 
with simulations, we cannot distinguish these shifts from 
random fluctuations; we therefore include the largest 
shift, Sw = 0.08, as a symmetric systematic uncertainty 
for these four sample combinations. 

Minimum Redshift for Nearby SN sample 

As discussed in §4, the choice of minimum redshift z„iin 
for the Nearby sample is complicated by the so-called 
"Hubble bubble," a jump in estimated SN la distance 
modulus around z ~ 0.025 within the Nearby sample 
(JRK07). This shift is shown in the residual Hubble di- 
agram in Fig. 20, which compares the fltted SN distance 
modulus in redshift bins, /ifit(z), with the calculated dis- 
tance modulus (/icaic(-z)) for a cosmological model with 
w ~ —\ and f^M = 0.3. Since /icaicl-z) is a smooth 
function of redshift, the jump between the two left-most 
points in the plot reflects a ~ 0.12 mag discontinuity in 
SN distance modulus {^f^t{z)) at z ~ 0.025. The signif- 
icance of the shift in fi between SNe at z < 0.025 and 
z > 0.025 is 2.4(7. JRK07 found a comparable shift, cor- 
responding to 0.14 mag. To further explore this issue, 
we also show the same residuals for the lower-redshift 
portion of the SDSS-H sample. The fitted SDSS-II SN 
distances are consistent with those of the z > 0.025 sub- 
set of the nearby SNe la and less consistent with the 
z < 0.025 subset 2.5a discrepancy). 

The Hubble bubble may be a real feature due to a 
local, large-scale void, an artifact of selection biases in 
the nearby SN sample, or an artifact of the light-curve 
fitter assumptions about host-galaxy dust and color vari- 
ations. Regardless of its interpretation, this feature has 
been used as an argument to discard SN la measure- 
ments at z < 0.025 from cosmological fits (e.g., Riess 
et al. 2007). We have decided to make a more agnos- 
tic choice for z,nin. We carry out cosmological fits us- 
ing two sample combinations. Nearby -fSDSS (c) and 
Nearby-hSDSS+ESSENCE-fSNLS (rf), and vary z^in 
from 0.015 to 0.03; for mlcs2k2, the resulting varia- 
tion in w is shown in the first two panels of Fig. 21. For 
both sample combinations, the approximate midpoint of 
the w-variation occurs at Zmin = 0.02, which we there- 
fore take as the nominal choice. As z„iin is varied around 
this value, the w- variations are approximately ±0.060 for 
sample combination (c) and ±0.040 for (d); we include 
these variations as systematic uncertainties in Table 6, 
with associated results for in Table 7. 

mlcs2k2 SN la Model Parameters 

The mlcs2k2 model vectors M, p, and q in Eq. 1 were 
determined by training the light-curve fitter on a nearby 
SN la sample, as described in JRK07. As discussed in 
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§5.1, our analysis uses a set of model vectors that includes 
the adjustments to the M'^'^ values given in Eq. 2. To 
estimate the sensitivity of the results to these adjust- 
ments, we have also carried out the cosmology analysis 
using the mlcs2k2 model vectors without those adjust- 
ments. These two sets of model vectors result in values 
for w that differ by 5w = 0.01 to 0.04, depending on the 
sample combination. 

Galactic Extinction 

The wavelength-dependence of the Milky Way Galactic 
extinction is expressed as ^(A) = Ay x {a{X) + b[X) / Ry), 
where Ry = 3.1, Ay = i?v x E{B - V), E{B - V) is 
the color excess, and the functions a(A) and 6(A) are de- 
fined by CardeUi, Clayton, & Mathis (1989). To estimate 
the systematic uncertainty for each SN, we coherently de- 
crease the color excess by 0.01 +0.16 X i?(i?—y) mag (but 
requiring non-negative color excess) relative to the nom- 
inal value in Schlegel, Finkbeiner, & Davis (1998). The 
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Fig. 20. — Difference between fitted mlcs2k2 SN distance modu- 
lus for low-redshift SNe, ii{i^(z), and the calculated distance mod- 
ulus /^calcC-^) for an ui = — 1, Hm = 0.3 model. The residuals are 
averaged in redshift bins of width 0.02 for the Nearby sample (black 
squares) and in bins of width 0.037 for the nearer portion of the 
SDSS-II sample (red circles). Errors are computed from O.lQ/y/N , 
where 0.16 is the typical magnitude dispersion per SN and N is 
the number of SNe in a given redshift bin. The dotted lines in- 
dicate the uncertainty in /icalc(^) resulting from an uncertainty in 
w of 0.15. The vertical arrow a,t z = 0.025 indicates the redshift 
associated with the Hubble bubble in JRK07. 



<N1 0.06 



I 

-0.02 



MLCS 

Nearby+SDSS 




SALT- 1 1 

Nearby+SDSS 

+ESSE+SNLS 




0.02 
Z,., 



0.02 



0.02 



0.02 



Fig. 21. — Variation in dark energy equation of state parameter 
w with minimum redshift 2„iin as the latter is varied around the 
nominal value of ^min = 0.02 (shown by dotted vertical line). The 
fitting method (mlcs2k2 or SALT-Il) and sample combination are 
indicated in each panel. 



color-excess uncertainty of 0.01 mag is based on optical- 
versus-IR discrepancies (Burstein 2003). The resulting 
lu-uncertainties are 0.01-0.02. Note that we have not 
varied the functions a(A) and 6(A), and therefore this 
uncertainty may be underestimated. The mean Galactic 
extinctions in the r-band are 0.20 (nearby), 0.14 (SDSS), 
0.07 (ESSENCE), 0.05 (SNLS) and 0.01 mag (HST). 

Dust Parameter Ry 

As described in § 7.1, the SDSS dust sample was used to 
determine the dust reddening parameter, Ry = 2.18 ± 
0.50 (§ 7.2), and the extinction distribution exponential 
slope, TV = 0.334 ±0.088 (§ 7.1). Propagating the corre- 
lated uncertainties on Ry and tv , the uncertainty on w is 
Sw = 0.036 for the SDSS-only sample. For the combined 
data samples, Sw < 0.03. This color-related uncertainty 
appears rather small, considering that the issue of color 
variations has been a major source of uncertainty in pre- 
vious studies such as WV07. 



Simulated Spectroscopic Efficiency 

The simulated efficiency is part of the mlcs2k2 fitting 
prior, as indicated by the terms Cscarch and ecuts h^ Eq. 5. 
The main effect of an error in determining the spectro- 
scopic efficiency is to introduce a bias in the mlcs2k2 
prior (Eq. 5). As discussed in §6, the largest uncertainty 
in simulating the efficiency is related to modeling of the 
incompleteness of the spectroscopic observations, Cspec- 

For the nearby, SDSS, and SNLS samples, egpec de- 
pends mainly on the Mdim component, i.e., Am ~ 1 in 
§ 6.2.1. We explore the systematic error due to spectro- 
scopic efficiency modeling for these samples by repeating 
the analysis with Cgpec set to unity, an extreme varia- 
tion from the fiducial efficiency model of §6.2.1. Since 
the ESSENCE sample favors a purely redshift-dependent 



(A. 



M 



in Eq. 14), setting eg 



1 has no im- 



pact in that case; we instead evaluate the systematic er- 
ror for that sample by replacing the purely z-dependent 
Cspcc with a purely dim-dependent espoc, i-e., by setting 
Am = 1. In these tests, the efficiency due to photomet- 
ric selection cuts (ecuts) is included in the fitting prior as 
well as the image-subtraction efficiency (egubtr) for the 
SDSS sample. 

For the above systematic changes in egpoc, the corre- 
sponding changes in distance modulus vs. redshift are 
shown in Fig. 22 for both the data and the Monte Carlo 
simulations. One sees clearly that incorporating a non- 
trivial spectroscopic efficiency model has a significant 
systematic impact on distance estimates for the more dis- 
tant SNe in each sample. Moreover, the data-simulation 
agreement is good, adding confidence in our implemen- 
tation of the spectroscopic efficiency model in the fitting 
prior. 

From Fig. 9, there is no doubt that egpoc is significantly 
less than one and that our model for egpoc is closer to re- 
ality than the extreme assumption of setting egpec = 1- 
To be conservative, we use half of the absolute- value dif- 
ference in w between the fiducial and extreme efficiency 
models as our estimate of the systematic error associated 
with uncertainty in the efficiency model. The largest un- 
certainty related to the simulated efficiency occurs for 
sample combinations (a) and (c); the w-uncertainties 
(half the shifts) are Sw = 0.062 for the SDSS-only sample 
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Fig. 22. — Systematic impact of uncertainty in the spectroscopic 
efficiency model for different samples. Upper panels: data sam- 
ples; lower panels: simulated samples. For each sample, the dif- 
ference in distance modulus between using an extreme spectro- 
scopic efficiency model (/x(syst tspcc)) and using the fiducial effi- 
ciency model (/x) is shown vs. redshift. For the nearby, SDSS, 
and SNLS data sets, the extreme model assumes no spectroscopic 
incompleteness (egpec = 1). For ESSENCE, the extreme model as- 
sumes pure Malmquist-induced incompleteness, and the negative of 
the /.i-difference is plotted instead in order to better compare with 
the results for the other samples. The simulated samples have 2.5 
times more points than the data samples for enhanced statistics. 
The differences in ii appear only for z > 2cut as explained in § 6.2.1. 

and 5w = 0.072 for the Nearby-|-SDSS sample combina- 
tion (c). We also note that including vs. not including 
the simulated efficiency in the nearby SN sample changes 
w by a few hundredths for the other sample combinations 
that include the Nearby sample, half of which is included 
as a systematic uncertainty. 

Calibration of Primary Reference Star, BD+17 
As discussed in Appendix B. the UBVRI magnitudes of 
the primary spectrophotomctric reference star, BD-f 17, 
are used in the K-corrections to relate the flux calibra- 
tion of the nearby SN sample to that of the higher- 
redshift samples. To evaluate the uncertainty in the 
BD-l-17 magnitudes, we compare the consistency of the 
synthetic magnitudes of the HST spectra with the SDSS 
photometric magnitudes. We first compute u,g,r,i,z 
zeropoint offsets as the difference between synthetic 
BD-l-17 magnitudes computed from the HST-measured 
spectrum (Bohlin 2006) and the magnitudes measured 
by Landolt (2007). These offsets are then compared 
to those based on three solar analogs (P3330E, P177D, 
P041C). The differences between the zeropoints off- 
sets are -0.004, -0.013, 0.005, 0.010, 0.012 mag for 
u, g, r, I, z, respectively. We therefore assume a 1% uncer- 
tainty in each of the U, B, V, R, I magnitudes for BD-l-17. 
To propagate this error, we change each BD-l-17 magni- 
tude by 1% independently in each passband; for each 
change, all of the SN light curves are re-fit with the 
mlcs2k2 model and cosmological parameters are ex- 
tracted. The resulting five independent w-shifts are then 
added in quadrature and included as a systematic uncer- 



tainty. The resulting lu-uncertainties are 0.02 ~ 0.03. 

As a crosscheck, we have replaced BD-f 17 with Vega 
as the primary reference star. The resulting changes in w 
are consistent with the changes from applying 1% shifts 
in the U,B,V,R,I magnitudes of BD-hl7. The differ- 
ences between BD-l-17 and Vega are therefore not in- 
cluded as a systematic uncertainty. 

Landolt-Bessell Color Transformations 

As discussed in Appendix B, we use color transforma- 
tions to transform between the Landolt network and the 
filter system defined by Bessell (1990). Table 20 shows 
the color terms (fc^) between Landolt magnitudes and 
synthetic magnitudes using the Bessell (1990) filters. We 
vary each color term independently by one standard de- 
viation, and the corresponding changes in w are added 
in quadrature. The resulting change is dw < 0.01 for all 
sample combinations. 

Shifted Bessell Filters instead of Color Transfor- 
mation 

As an alternative to the Landolt-Bessell color transfor- 
mation, we consider the approach of Astier et al. (2006), 
using a modified set of UBVRI filter responses in which 
the shapes of the Bessell (1990) response curves are held 
fixed, but the central wavelength of each filter passband 
is shifted so that the color terms (Eq. B2) are zero. The 
corresponding wavelength shifts are given in Table 21 
under the column "HST standards." Using these wave- 
length shifts and no color transformations in place of the 
color transformation method of Appendix B, w shifts by 
Sw = 0.007 for the SDSS-only sample. For the other 
sample combinations, Sw ^ 0.01. 

SDSS AB Offsets 

As discussed in §3.2, the uncertainties in the SDSS AB 
offsets are 0.009, 0.009, 0.009 mag for g, r, i, respectively 
(see Table 1). These uncertainties include zeropoint un- 
certainties of 0.004, 0.004, 0.007 mag based solely on the 
consistency of the solar analogs and also account for un- 
certainties in the central wavelengths of the SDSS filters: 
SX ~ 7, 16, 25 A for g, r, z, respectively. We indepen- 
dently vary each AB offset by one standard deviation, 
and the resulting changes in w are added in quadrature. 
The resulting change in w; is a few hundredths for the 
sample combinations that include SDSS-II SNe. Note 
that the w-uncertainty for SDSS-only {5w = 0.004) is 
much smaller than that for the Nearby-|-SDSS combina- 
tion {5w = 0.030). For SDSS-only, the effect of chang- 
ing the AB offsets is mainly to shift all of the distance 
moduli by the same amount; the corresponding change 
in cosmological parameters is small. For Nearby -l-SDSS, 
the SDSS-II distances are shifted relative to those of the 
nearby sample, and the abrupt feature in the Hubble 
diagram results in a larger change in the cosmological 
parameters. 

ESSENCE R-I Color Zeropoint 

WV07 report an uncertainty of 0.02 mag in their R — I 
color zeropoint, resulting in a systematic uncertainty of 
Sw = 0.04 for their analysis of the nearby-f-ESSENCE 
data sets. As a crosscheck, we propagate their R ~ I 
uncertainty into our analysis and find Sw ~ 0.05 for 
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the same nearby+ESSENCE combination, in reasonable 
agreement with WV07. We apply their R — I uncertainty 
to our sample combinations and find that the uncertainty 
in w is much smaller, Sw 0.01 for the combinations 
{b,d,e,f) that include the ESSENCE sample. Note that 
these sample combinations also include the SNLS sam- 
ple, suggesting that the impact of the ESSENCE R~ I 
uncertainty is reduced by the presence of another high- 
redshift data sample. 

SNLS g, r, i, z Zeropoints 

Astier et al. (2006) report zeropoint uncertainties of 0.01, 
0.01, 0.01, 0.03 for their g,r,i,z passbands, respectively. 
In their analysis of the nearby+SNLS combination, they 
vary each zeropoint independently and add the corre- 
sponding w-shifts in quadrature: the resulting system- 
atic is Sw = 0.05. As a crosscheck, we propagate their 
zeropoint uncertainties into our analysis and also find 
Sw = 0.05 for the same nearby+SNLS combination. We 
apply their zeropoint uncertainties to our sample com- 
binations and find slightly smaller uncertainties for the 
combinations {b,d,e,f} that include the SNLS sample. 

HST Zeropoints 

The HST zeropoint uncertainties are 0.02 mag for the 
FllOW and F160W filters^'^ (mean wavelengths are 
1.14 fim and l.Gl/im, respectively), and 0.01 mag for 
the optical filters.''^ For the analysis of the five-sample 
combination that includes HST (e), the resulting w- 
uncertainty is Sw = 0.012. 

The cosmological parameter shifts due to all of these 
sources of systematic error are added in quadrature to 
derive total systematic error estimates Sw and S^m for 
the mlcs2k2 analysis of the FwCDM model; these are 
given in Tables 6 and 7. 

9.2. Systematics Uncertainties with SALT-II 

To examine systematic uncertainties in the context of 
the SALT~ii model, we undertake an analysis similar to 
that carried out for mlcs2k2. We also determine sys- 
tematic uncertainties for the SALT-ii parameters a and 
P that enter Eq. 29. Since the SALT-II training soft- 
ware is not available for public use, with one exception 
we make approximations in cases where re-training of the 
spectral surfaces is needed. In such cases we either use 
the nominal SALT-ii surfaces and propagate changes only 
in the light-curve fits, or we use the uncertainties based 
on the mlcs2k2 analysis. The systematic uncertainties 
in w and flu for the different sample combinations are 
given in Tables 8 and 9. 

Rest frame [/-band 

As discussed in detail in § 10.2.4, there are systematic dis- 
crepancies between the SALT-ii rest-frame [/-band model 
and the observer- frame [/-band light-curve data for the 
nearby SN la sample. As was noted in § 9.1, a related is- 
sue is seen for mlcs2k2. We carry out a test of the salt- 
II fits in which the observer-frame filter corresponding to 
rest-frame [/-band is excluded from the fits. Figure 40 

Sec HST Handbook for NICMOS, http://www.stsci.edu/hst 
/nicmos/documents/handbooks/handbooks/DataHandbookv? 

Sec HST Handbook for ACS, http://www.stsci.edu/hst 
/acs/documents/haiidbooks/DataHandbookv4/ACS_longdhbcover .html 



shows that for the SDSS-only sample, the exclusion of 
[/-band results in a redshift-dependent change in the dis- 
tance modulus. This results in a w-shift of —0.100 for 
the SDSS-only sample and -0.133 for the Nearby-^SDSS 
combination; these shifts are included as asymmetric sys- 
tematic uncertainties. 

For the other sample combinations, which include the 
higher-redshift ESSENCE & SNLS samples (&,d,e,/), the 
exclusion of rest-frame [/-band results in w-shifts of .04 
to 0.09. Based on tests with simulations, we cannot dis- 
tinguish these shifts from random fluctuations; we there- 
fore include the largest shift, Sw = 0.09, as a symmetric 
systematic uncertainty for these four sample combina- 
tions. If observer-frame [/-band is excluded from the 
nearby sample, while rest-frame [/-band is included in 
the higher-redshift samples, the change in w is no more 
than 0.01. 

Minimum Redshift for Nearby SN sample 

The dependence of the SALT-ii results for w on the z„iin 
cut is shown in the third and fourth panels of Fig. 21. For 
■^min ^ 0.018, the inferred value of w is fairly insensitive 
to the value of ^min- However, when Zmin is reduced to 
0.015, w changes by almost 0.09 for the Nearby+SDSS 
combination, and by 0.04 - 0.05 for the combinations 
that include the higher-redshift (ESSENCE & SNLS) 
samples. To account for these variations, we have as- 
signed a systematic uncertainty of Sw ~ 0.05 for the 
Nearby+SDSS combination and Sw = 0.03 for the other 
sample combinations that include the nearby SN la sam- 
ple. 

Galactic Extinction 

The systematic uncertainty from Galactic extinction is 
described in § 9.1, resulting in Sw ^ 0.02. 

SALT-ii Training v^fith SDSS-II Data 

Here we examine the SALT-ii training process that pro- 
duces the spectral surfaces Mo{t,X) and Mi{t,X) in Eq. 
6. Because SDSS-II SN probes a relatively unexplored 
range of SN redshifts, the rest-frame behavior of the 
SDSS-II light curves may not be as well described by 
the SALT-II model as that of other SN samples that were 
used in the training of the model. To quantify this is- 
sue, J. Guy has retrained the SALT-ii spectral surfaces 
twice, first including the light curves of the SDSS spec- 
troscopic sample and second including those of the SDSS 
dust sample (§ 7.1). Evaluating cosmological parame- 
ters obtained with each retrained set of spectral surfaces 
and comparing the results with those from the standard 
SALT-II training, we include the larger of the two w- 
shifts as a systematic uncertainty. For the SDSS-only 
and Nearby+SDSS (a & c) sample combinations, the 
uncertainty is Sw 0.02. For the other combinations, 
Sw - 0.01. 

SALT-II Dispersions 

Recall from § 5.2 that the spectral surfaces, Mo{t, A) and 
Mi(t, A), were retrained using the Bessell filter shifts 
based on HST standards (Table 21). The model dis- 
persions around these surfaces, however, were not deter- 
mined in the retraining, and we therefore use the model 
dispersions from Guy et al. (2007). To allow for the 
resulting uncertainty in the dispersions, we assign a sys- 
tematic uncertainty of half the difference between using 
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TABLE 6 

Systematic uncertainties in dark energy equation of state parameter w for the mlcs2k2 analysis of the FuiCDM model, 
including the bao + cmb prior. negative values indicate asymmetric uncertainties. 



source oi uncertainty 


a, 


uncertainty on w for sample: 
b c d e 


/ 


Zmin cut for Nearby sample 
mlcs2k2 SN la Model Parameters 
Galactic Extinction 


—0.310 
0.000 
0.013 
0.021 


0.080 
0.000 
0.036 
0.012 


-0.310 
0.060 
0.001 
0.004 


0.080 
0.040 
0.040 
0.016 


0.080 
0.040 
0.026 
0.022 


0.080 
0.040 
0.043 
0.023 


FORM OF PRIOR 














correlated la changes Ry and ry 


0.036 


0.023 


0.006 


0.007 


0.008 


0.009 


simulated efficiency for nearby SN la 


0.000 


0.000 


0.014 


0.008 


0.007 


0.013 


spectroscopic efficiency for SF)SS 




0.002 


0.072 


0.007 


0.001 


n nnn 
u.uuu 


spectroscopic efficiency for ESSENCE 


0.000 


0.013 


0.000 


0.012 


0.008 


0.012 


spectroscopic efficiency for SNLS 


0.000 


0.029 


0.000 


0.025 


0.017 


0.026 


CALIBRATION 














0.01 mag errors in U, B, V, R, I 


0.030 


0.029 


0.030 


0.021 


0.019 


0.021 


shifted Bessell filters 


0.007 


0.017 


0.007 


0.016 


0.013 


0.014 


vary ki color terms 


0.003 


0.007 


0.005 


0.006 


0.006 


0.006 


vary SDSS AB offsets for g, r, i 


0.004 


0.028 


0.030 


0.012 


0.013 


0.000 


vary ESSENCE R - I color zeropoint 


0.000 


0.011 


0.000 


0.015 


0.008 


0.010 


vary SNLS g, r, i, z zeropoints 


0.000 


0.040 


0.000 


0.036 


0.024 


0.035 


vary HST zeropoints 


0.000 


0.000 


0.000 


0.000 


0.012 


0.000 


Total 


+0.08 
-0.32 


0.11 


+0.10 
-0.33 


0.11 


0.11 


0.12 



" SDSS-only 

' SDSS+ESSENCE+SNLS 

" Nearby+SDSS 
Nearby+SDSS+ESSENCE+SNLS 
Nearby+SDSS+ESSENCE+SNLS+HST 

f Nearby+ESSENCE+SNLS 



TABLE 7 

Systematic uncertainties in matter density parameter Qm for the mlcs2k2 analysis of the FuiCDM model, including the 
BAO+CMB PRIOR. Negative values indicate asymmetric uncertainties. 



source of uncertainty 


a 


uncertainty on O 
b c 


M for sample; 
d e 


/ 


Rest frame [/-band 


-0.051 


0.016 


-0.051 


0.016 


0.016 


0.016 


2min cut for Nearby sample 


0.000 


0.000 


0.014 


0.009 


0.009 


0.009 


mlcs2k2 SN la Model Parameters 


0.003 


0.009 


0.000 


0.010 


0.006 


0.011 


Galactic Extinction 


0.005 


0.003 


0.001 


0.004 


0.006 


0.006 


FORM OF PRIOR 














correlated la changes Ry and tv 


0.007 


0.006 


0.002 


0.002 


0.003 


0.002 


simulated efficiency for nearby SN la 


0.000 


0.000 


0.003 


0.002 


0.002 


0.003 


spectroscopic efficiency for SDSS 


0.014 


0.001 


0.015 


0.002 


0.000 


0.000 


spectroscopic efficiency for ESSENCE 


0.000 


0.004 


0.000 


0.002 


0.002 


0.002 


spectroscopic efficiency for SNLS 


0.000 


0.008 


0.000 


0.007 


0.004 


0.007 


CALIBRATION 














0.01 mag errors in U, B, V, R, I 


0.006 


0.008 


0.006 


0.005 


0.005 


0.005 


shifted Bessell filters 


0.001 


0.004 


0.001 


0.004 


0.003 


0.003 


vary ki color terms 


0.001 


0.002 


0.001 


0.002 


0.002 


0.001 


vary SDSS AB offsets for g, r, i 


0.001 


0.007 


0.006 


0.003 


0.003 


0.000 


vary ESSENCE R - I color zeropoint 


0.000 


0.003 


0.000 


0.004 


0.002 


0.002 


vary SNLS g, r, i, z zeropoints 


0.000 


0.010 


0.000 


0.009 


0.005 


0.008 


vary HST zeropoints 


0.000 


0.000 


0.000 


0.000 


0.003 


0.000 


Total 


+0.019 
-0.054 


0.027 


+0.023 
-0.056 


0.026 


0.023 


0.026 
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and ignoring the dispersions. The resulting uncertain- 
ties arc Sw ^ 0.01 to 0.02 for combinations that include 
the higher-redshift ESSENCE & SNLS samples. For the 
SDSS-only and Nearby+SDSS sample combinations, the 
w-uncertainty is negligible. 

/3- Variation with Redshift 

If the SALT-ii SN parameters (a,/3,M) are allowed to 
vary independently in redshift bins, while the cosmo- 
logical parameters are fixed, we find a strong redshift- 
dependence of (3 for z > 0.6 (see § 10.2.3). To estimate 
the corresponding systematic uncertainty, the Hubble di- 
agram fits have been redone allowing a, /3, and M to vary 
with redshift using a simple model in which each SN 
parameter is constant for z < 0.6, and is then allowed 
to vary linearly with redshift for z > 0.6. Compared 
to the nominal SALT-ii model with redshift-independent 
parameters, the largest change, Sw = -1-0.073, occurs 
for SDSS-f-ESSENCE+SNLS (6) in which the nearby 
sample is excluded. For sample-combinations d and /, 
5w ^ 0.04, and for the full SN set (e) that includes the 
HST, Sw ~ 0.01. These w-shifts are included as asym- 
metric systematic uncertainties. 

Simulated Bias Correction 

For SALT-II, we have determined bias corrections from 
simulations, as described in § 8.2 (sec Table 5). We in- 
clude half the w-shift as a systematic uncertainty. The 
largest uncertainty is Sw = 0.020 for SDSS-only. 

Calibration of Primary Reference Star, Vega 

We assume uncertainties of 0.01 magnitudes in the cal- 
ibration of U, B, V, R, I for the primary reference, Vega. 
Since a full accounting of this effect would require an- 
other retraining of the SALT-ii surfaces, we instead adopt 
the uncertainties derived from the MLCS2k2 analysis (Ta- 
ble 6). In Astier et al. (2006), the corresponding uncer- 
tainty is Sw = 0.024 for the nearby-|-SNLS combination; 
as a crosscheck, we have evaluated this uncertainty for 
the same sample combination and find good agreement, 
Sw ~ 0.021. For the sample combinations analyzed here, 
the resulting iw-uncertainties are 0.02 - 0.03. 

Calibration: Shifted Bessell Filters for Nearby 
Data 

As discussed in § 5.2, we use the Bessell (1990) filter re- 
sponses with wavelength shifts given in Table 21 of Ap- 
pendix B. Since these shifts differ from those in Astier 
et al. (2006), we use the difference in cosmological results 
derived from both sets of wavelength shifts to define an 
additional systematic uncertainty on w. This uncertainty 
is Sw ~ 0.02 for sample combinations that include the 
nearby sample. 

Calibration: Zeropoint offsets for SDSS, 
ESSENCE, SNLS, HST 

Zeropoint uncertainties for the SDSS, ESSENCE, 
SNLS, and HST bandpasses are propagated in the same 
manner as for the mlcs2k2 method (§ 9.1). The SNLS 
zeropoints are varied in the fit, but not in the training of 
the spectral surface, and therefore these w-uncertainties 
might be slightly overestimated. The other survey 
samples were not used in the training, and therefore 



varying the zeropoints in the fit is sufficient to estimate 
the systematic error. Note that the ly-uncertainty from 
the SDSS AB offsets is smaller for SDSS-only than for 
Nearby-|-SDSS, as explained in § 9.1. 

The parameter shifts due to all of these systematic er- 
rors are added in quadrature to derive total systematic 
error estimates in Tables 8 and 9. 

10. RESULTS 

Here we present the Hubble diagram and inferred cos- 
mological parameters using the framework described in 
§8. Results based on the mlcs2k2 and SALT-II methods 
are presented in §10.1 and 10.2 respectively. We compare 
the results from the two methods in § 11. 

The Hubble diagram for the five samples considered in 
this analysis is shown in Fig. 23. The distance moduli 
here are obtained from the mlcs2k2 method described 
above; the Hubble diagram based on the SALT-ii method 
looks quite similar. Detailed Hubble-residual plots are 
given for each method in §10.1 and 10.2. 
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Fig. 23. — Fitted distance modulus (from mlcs2k2) versus red- 
shift for the 288 SNe la from the five samples indicated on the 
plot. 



10.1. Results with mlcs2k2 

Using the mlcs2k2 method, we present cosmological 
results for the six sample combinations {a-f) of Table 4. 
Table 10 gives the spectroscopically determined redshift 
and marginalized mlcs2k2 fit parameters for SNe that 
pass the selection cuts described in §4. We use the en- 
semble of redshifts Zi and estimated distance moduli 
for each sample combination to fit cosmological model 
parameters, as explained in §8.1. 

10.1.1. Goodness of Fit and Hubble scatter 

Before considering the cosmological parameter results 
for the various combined samples, we examine several 
measures of fit quality and dispersion for each SN la sam- 
ple treated independently, since they provide diagnostic 
information that is useful to consider before combining 
the samples. Table 11 displays these measures: (i) the 
statistic for the best-fit FwCDM model for that sample 
from Eq. 26 — in goodness-of-fit tests, this statistic is 
usually compared to the number of degrees of freedom. 
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TABLE 8 

Systematic uncertainties in w for the salt-ii analysis of the FuiCDM model, including the BAO+CMB prior. +/— values 

indicate asymmetric uncertainties. 







Uncertainty 


on w for Sample: 






Source of Uncertainty 


a 


b 


c 


d 


e 


/ 


Rest frame f7-band 


-0.100 


0.104 


-0.133 


0.104 


0.104 


0.104 


•^min cut for Nearby sample 


0.050 


0.030 


0.050 


0.030 


0.030 


0.030 


Galactic Extinction 


0.021 


0.012 


0.004 


0.016 


0.022 


0.023 


SALT-II SN la MODEL PARAMETERS 














retraining : include SDSS data 


0.008 


0.005 


0.017 


0.011 


0.005 


0.005 


dispersions of SALT-II surfaces 


0.001 


0.003 


0.002 


0.006 


0.006 


0.004 


/3-variation with redshift 


0.000 


+0.073 


0.000 


+0.045 


+0.013 


+0.036 


SELECTION EFFICIENCY 














simulated bias 


0.020 


0.011 


0.009 


0.002 


0.001 


0.012 


CALIBRATION 














0.01 mag errors in U, B, V, R, I 


0.029 


0.030 


0.027 


0.022 


0.020 


0.022 


shifted Bessel90 filters 


0.000 


0.000 


0.015 


0.010 


0.008 


0.013 


vary SDSS AB offsets for g, r, i 


0.018 


0.037 


0.031 


0.015 


0.016 


0.000 


vary ESSENCE R - I color zeropoint 


0.000 


0.035 


0.000 


0.036 


0.021 


0.025 


vary SNLS g, r, i, z zeropoints 


0.000 


0.057 


0.000 


0.046 


0.030 


0.043 


vary HST zeropoints 


0.000 


0.000 


0.000 


0.000 


0.015 


0.000 


Total 


+0.06 


+0.15 


+0.07 


+0.13 


+0.12 


+0.13 


-0.12 


-0.14 


-0.15 


-0.13 


-0.12 


-0.12 



" SDSS-only 
SDSS+ESSENCE+SNLS 
Nearby+SDSS 

Nearby+SDSS+ESSENCE+SNLS 
"= Nearby+SDSS+ESSENCE+SNLS+HST 
f Nearby+ESSENCE+SNLS 



Systematic uncertainties in Qj, 



TABLE 9 

for the SALT-II ANALYSIS OF THE FuiCDM MODEL, INCLUDING THE BAOH 
INDICATE ASYMMETRIC UNCERTAINTIES. 



CMB PRIOR. 



VALUES 







Uncertainty 


on 


for Sample: 






Source of Uncertainty 


a 


b 


c 


d 


e 


/ 


Rest frame f/-band 


-0.020 


0.022 


-0.024 


0.022 


0.022 


0.022 


^min cut for Nearby sample 


0.012 


0.007 


0.012 


0.007 


0.007 


0.007 


Galactic Extinction 


0.005 


0.003 


0.001 


0.004 


0.006 


0.006 


SALT-II SN la MODEL PARAMETERS 














retraining : include SDSS data 


0.002 


0.001 


0.003 


0.001 


0.002 


0.001 


dispersions of SALT-II surfaces 


0.000 


0.001 


0.000 


0.001 


0.001 


0.001 


/3-variation with redshift 


0.000 


+0.016 


0.000 


+0.010 


+0.002 


+0.007 


SELECTION EFFICIENCY 














simulated bias 


0.005 


0.003 


0.002 


0.003 


0.000 


0.003 


CALIBRATION 














0.01 mag errors in U, B, V, R, I 


0.006 


0.008 


0.005 


0.005 


0.005 


0.005 


shifted Besscl90 filters 


0.000 


0.000 


0.003 


0.002 


0.001 


0.002 


vary SDSS AB offsets for g,r,i 


0.004 


0.007 


0.006 


0.003 


0.003 


0.000 


vary ESSENCE R - I color zeropoint 


0.000 


0.006 


0.000 


0.006 


0.003 


0.004 


vary SNLS g, r, i, z zeropoints 


0.000 


0.011 


0.000 


0.009 


0.005 


0.008 


vary HST zeropoints 


0.000 


0.000 


0.000 


0.000 


0.003 


0.000 


Total 


+0.015 


+0.033 - 


f0.015 


+0.028 


+0.025 


+0.027 


-0.025 


-0.029 - 


-0.029 


-0.027 


-0.025 


-0.026 



given by A'dof = Nsn — 1;*^ (ii) the root-mean-square 
measure of Hubble scatter, RMS^ = ^ 

where is the estimated distance modulus from the 
light-curve fit for the i'th SN, and /if is the best-fit 
FwCDM model distance modulus at the corresponding 

Af(jof = number of SNc minus the number of cosmology pa- 
rameters {Hq, w, Qm) + the number of priors (BAO+CMB). 



redshift zt; and (iii) the vahie of a™*^ that would be 
required in Eq. 27 to make xfi = ^dof in Eq. 26 for 
that sample. In computing the first two measures, we 
adopt ctJj"* = 0.16, the value that yields = -^dof 
for the nearby SN sample and that we use in analyz- 
ing the combined samples, as explained in §8.1. The 
bottom row of Table 11 shows the statistic from 
the SNe in each sample when the best-fit FwCDM 
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model parameters for the global sample combination 
Nearby+SDSS+ESSENCE+SNLS+HST (e) are used to 
determine fif*". Compared to the values for the indepen- 
dent fits to each sample, the xfi values for the global 
fit are only shghtly larger. For the nearby, ESSENCE, 
SNLS, and HST samples, the reduced statistic, xfi/^doi, 
is close to unity, and RMS^ = 0.f9 to 0.28 mag. For 
the SDSS-II sample, the reduced x^i/-^dof ~ 0.5, and 
RMSp = 0.f5, both substantially smaller than for the 
other samples. 

We attribute the smaller scatter and xfi of the SDSS-II 
sample largely to spectroscopic selection effects. As de- 
scribed in § 2, when prioritizing candidate SNe for spec- 
troscopic observations, preference was given to those that 
were far from host-galaxy cores and/or that were hosted 
by redder (and presumably less dusty) galaxies. In Ap- 
pendix E, this explanation is quantified by comparing 
the Hubble scatter and xfi statistic for the spectroscopic 
SN la sample to those for the host-z sample of photomet- 
rically identified SNe la that have spectroscopic redshifts 
from subsequent host-galaxy observations. The cosmol- 
ogy analysis accounts for this selection effect via the 
model for the search efficiency (§ 6.2.1) in the mlcs2k2 
fitting prior and by including a systematic error that re- 
flects uncertainties in the search efficiency. 

10.1.2. mlcs2k2 Hubble Diagrams and Cosmological 
Parameters 

Figure 24 shows the differences between the estimated 
SN distance moduli fii and those for an open CDM model 
with no dark energy (f2M = 0.3, IIde = 0) as a func- 
tion of redshift, for each of the six SN sample combina- 
tions; the large square (pink) points show weighted av- 
erages of these residuals in redshift bins. Also shown is 
the Hubble distance-residual curve between the best-fit 
FwCDM model for each sample combination (including 
the BAO-I-CMB prior) and that for the open model (solid 
curves). In each panel, the Hubble parameter for the 
open CDM model has been adjusted to agree with that 

TABLE 10 

Parameters from mlcs2k2 light curve fits (uncertainties 
IN parentheses). The complete table for all 288 SNe is 

GIVEN in electronic FORM IN THE .lOURNAL, AND ALSO AT 
http : //das . sdss . org/va/SNcosmology/siicosm09_f its . tar .gz. 



SNID 


redshift''^ 




Av 


A 




762 


0.1904(.0001) 


40.05(0.10) 


0.20(0.08) 


-0.22(0.07) 


53624.4(0.4) 


1032 


0.1291(.0002) 


38.80(0.10) 


0.07(0.06) 


0.88(0.09) 


53626.9(0.2) 


1112 


0.2565(.0002) 


40.84(0.18) 


0.13(0.09) 


0.05(0.16) 


53629.3(1.0) 


1166 


0.3813(.0005) 


41.51(0.18) 


0.16(0.11) 


-0.25(0.13) 


53630.1(1.1) 


1241 


0.0858(.0050) 


38.10(0.09) 


0.40(0.06) 


0.02(0.06) 


53634.7(0.2) 


1253 


0.2609(-0050) 


40.65(0.14) 


0.06(0.05) 


0.08(0.12) 


53634.2(0.5) 



for the best-fit FwCDM model, so that the FwCDM ver- 
sus open CDM residuals vanish at z = 0. Since the Hub- 
ble parameter is not determined in this analysis, a con- 
stant vertical offset in Fig. 24 is irrelevant: what is sig- 
nificant are the slope and curvature of the points and the 
best-fit (solid) curves vs. redshift. Figure 25 shows the 
distributions of normalized residuals, (/li — /-iFiuCDM)/^/^, 
where ^ywCHM is the distance modulus from the best-fit 
FwCDM model for sample combination (e) including the 
BAO+CMB prior, and is the total uncertainty defined 
in Eq. 27. The bulk of the distribution of all 288 nor- 
malized residuals (upper left panel of Fig. 25) is well fit 
by a Gaussian with a = 0.77; outliers increase the rms 
to 0.90. For the nearby sample, the rms is one as ex- 
pected, because a™*" = 0.16 mag is determined such that 

X^/A^dof = 1. 





redshift redshift 



Fig. 24. — Hubble residuals for the mlcs2k2 method: differences 
between measured SN distance moduli and those for an open CDM 
model (Om = 0.3, Hde = 0) vs. redshift for the six SN sample 
combinations. Large, square (pink) points show weighted averages 
in redshift bins: within each redshift bin, the points are plotted 
at the weighted mean redshift given by 2 = (X) •^i/<''?)/(5Z l/""?)! 
where cri is the distance-modulus uncertainty. Solid curves show 
the difference between the best-fit FuiCDM model distance mod- 
ulus and that for the open model, normalized to have the same 
value of the Hubble parameter. The error bars on the data points 
correspond to the distance modulus error crjj* from the mlos2k2 
light-curve fit (Eq. 27), i.e., they do not include the intrinsic scatter 
or the effects of redshift and peculiar velocity errors. The vertical 
bar in panel (a) shows the intrinsic uncertainty, crjf* = 0.16, in- 
cluded in the cosmology fits so that the P^r degree of freedom is 
unity for the Hubble diagram constructed from the nearby SN la 
sample alone. 

Figs. 26 and 27 show the mlcs2k2 statistical- 
uncertainty contours for the FwCDM and ACDM mod- 
els: for each of the six sample combinations, the SN la, 
BAO, and CMB contours are shown along with the com- 
bined constraints. For the combined SN-I-BAO-I-CMB 
results, we include systematic uncertainties to derive to- 
tal (statistical plus systematic) uncertainty contours, as 
explained in Appendix F, with results shown in Figs. 28 
and 29 for the FwCDM and ACDM models (note the 



^Spectroscopic redshift in CMB frame. 



TABLE 11 

Hubble diagram fit-quality parameters using mlcs2k2 
distances. 



fit-quality 
parameter 



A^dof 
RMS„ 



Result for sample: 
Nearby SDSS ESSENCE SNLS HST 



^ A^dof) 



31.9 


55.3 


46.8 


63.0 


32.5 


32 


102 


55 


61 


33 


0.19 


0.15 


0.23 


0.24 


0.28 


0.16 


0.09 


0.13 


0.16 


0.15 


32.9 


56.6 


48.4 


64.6 


32.4 



Xi (global fit) 
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Fig. 25. — Distribution 
residuals (pull) for the mlcs2k2 method, for the 
Nearby+SDSS+ESSENCE+SNLS+HST sample combination 
e (upper left, along with Gaussian fit) and for each SN sample 
indicated on the plots, fj, is the measured SN distance modulus, 
MFujCDM is the distance modulus from the best-fit FtoCDM model 
(for sample combination e) at the same redshift, and a-^ is the 
total uncertainty (Eq. 27). Mean and RMS for each distribution 
are indicated on each panel. 



zoomed axis scales compared to the previous figures). 
The best-fit cosmological parameter values and uncer- 
tainties, including the BAO and CMB priors, are given 
in Tables 12 and 13 for the FwCDM and ACDM models. 
The distance modulus vs. redshift curve for the best-fit 
FwCDM cosmological parameter values (relative to that 
of an open CDM model) are shown as the solid curves in 
the Hubble residual plots in Fig. 24. 

Among the six SN sample combinations, the best-fit 
values of the dark energy equation of state parameter 
w fall roughly into two groups. In the first group, the 
highest-rcdshift sample is from SDSS-II: for the SDSS- 
only sample (a) and Nearby -|~SDSS sample combination 
(c), w ~ —0.84 and —0.92. The agreement between 
these values is consistent with the expected RMS spread 
of 0.07 based on simulations. The second group com- 
prises the other four sample combinations, which in- 
clude the higher-redshift ESSENCE & SNLS samples: 
w = —0.71, —0.76, —0.76, —0.78 for sample combinations 
b,d,e,f. Simulations predict an RMS spread in w of ~ 0.1 
between the results from these two groups; the observed 
difference is therefore not statistically significant but may 
nevertheless be an indicator of systematic effects. 

Table 12 also shows that the statistical and system- 
atic errors in w and Qm for sample combinations (6) and 
(/) are very similar. Since these two sample combina- 
tions differ only in the substitution of the nearby SN la 
sample with the SDSS-II sample, this indicates that the 
first-season SDSS-II SN sample anchors the Hubble dia- 
gram with comparable constraining power to the nearby 
sample. 

Using the Nearby-hSDSS-f ESSENCE-hSNLS+HST 
sample combination ( e) , which covers the widest redshift 
range, we obtain w = — 0.76 ±0.07(stat) ±0.11 (syst) and 
Qm = 0.307 ± 0.019(stat) ± 0.023(syst). Although this 
value for w is higher than that obtained from other re- 
cent SN measurements, we stress that the difference is 
not due to inclusion of the SDSS-II SN data: as Table 12 



shows, we infer nearly identical parameter values for sam- 
ple combination (/), which excludes the SDSS-II data. 
By contrast, WV07 inferred w = -1.07 ± 0.09(stat) ± 
0.13(syst), nu = 0.267t{5;^i^(stat) using mlcs2k2 for a 
sample combination nearly identical to (f) and including 
the BAO (but not CMB) constraints. In § 10.1.4, we 
trace the differences between the WV07 result and ours 
to changes in mlcs2k2 model parameters and assump- 
tions. 

For the ACDM model, in comparison with the FwCDM 
model, the SN results carry less weight relative to the 
combined BAO and CMB results in constraining the pa- 
rameters. In particular. Fig. 27 and Table 13 show that 
the SDSS-only (a) and Nearby-hSDSS (c) SN samples 
have almost no impact on the maximum likelihood pa- 
rameter values and uncertainties. For the other four SN 
sample combinations, there is some tension between the 
SN and BAO results: the SN results pull the maximum 
likelihood parameter values along the CMB contour, 
away from the BAO contours. Since the BAO contours 
arc in all cases narrower than those for the SNe, how- 
ever, this shift is small, corresponding to SH^e — —0.03, 
SriM — 0.04 or less. 




Fig. 26. — For the FuiCDM model, mlcs2k2 statistical- 
uncertainty contours in the CM-f plane for each of the six SN 
sample combinations indicated on the plots. Shaded regions: 68% 
and 95% confidence level regions for the SN data alone; green: 
corresponding CL contours for SDSS BAO (Eisenstein ct al. 2005); 
blue: CL contours for WMAP5 CMB (Komatsu et al. 2008); closed, 
black contours: combined constraints from SN-|-BAO-|-CMB. 



10.1.3. U-Band Anomaly with mlcs2k2 

As noted in §9.1, the largest single contribution to the 
systematic error budget comes from consideration of the 
rest-frame [/-band. Within the mlcs2k2 framework, the 
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— MLCS total-error contour 



stat-error contour 



Fig. 27. — For the ACDM model, mlcs2k2 statistieal-uneertainty 
contours in the Qm-^A plane for each of the six SN sample 
combinations indicated on the plots. Shaded regions: 68% and 
95% confidence level regions for the SN data alone; green: cor- 
responding CL contours for SDSS BAO; blue: CL contours for 
WMAP5 CMB; closed, black contours: combined constraints from 
SN-I-BAO+CMB. Grey region indicates models with no Big Bang, 
i.e., with a bounce at finite value of the FRW scale factor. Solid 
diagonal line indicates a spatially flat Universe, with f2M + f^A = 1- 
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Fig. 28.— For mlcs2k2 and the F«;CDM model, 68% CL con- 
tours in the Q,y\-w plane for each of the six SN sample combina- 
tions, using the combined SN-^BAO-f CMB constraints. Solid eon- 
tours are total (statistical+systematic) uncertainty; dashed con- 
tours are statistical only. Systematic errors have been included 
using the prescription in Appendix F. 
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Fig. 29.— For mlcs2k2 and the ACDM model, 68% CL contours 
in the Qm-^^A plane for each of the six SN sample combinations, 
using the combined SN-I-BAO+CMB constraints. Solid contours 
are total (statistical+systematic) uncertainty; dashed contours are 
statistical only. Systematic errors have been included using the 
prescription in Appendix F. 

issue is manifest as a difference between the nearby sam- 
ple (JRK07) and the other SN la samples. We have 
carried out a series of tests in which the observer-frame 
passband corresponding to rcst-framc U -band is excluded 
from the mlcs2k2 light-curve fits. We compare the re- 
sulting distance modulus estimates, fJ-noU-, with those 



TABLE 12 

For the FwCDM model, constraints on w and Qm 
mlcs2k2 SN distances combined with SDSS BAO 
WMAP5 CMB RESULTS. 



from 

AND 





Result for sample combination 








a 


b 


c 


d 


e 


/ 




170.9 


406.7 


279.5 


517.8 


568.1 


341.9 


^dof 


102 


220 


135 


253 


287 


150 


RMS^ 


0.15 


0.20 


0.16 


0.20 


0.21 


0.23 




-0.84 


-0.71 


-0.92 


-0.76 


-0.76 


-0.78 


fJu, (stat) 


0.15 


0.09 


0.13 


0.08 


0.07 


0.08 


aw (syst) 


+0.08 
-0.32 


0.11 


+0.10 
-0.33 


0.11 


0.11 


0.12 


am (tot) 


+0.17 
-0.35 


0.14 


+0.16 
-0.35 


0.14 


0.13 


0.14 




0.289 


0.319 


0.273 


0.306 


0.307 


0.302 


o-Qm (stat) 


0.033 


0.025 


0.028 


0.021 


0.019 


0.022 


o"nM (syst) 


+0.019 
-0.054 


0.027 


+0.023 
-0.056 


0.026 


0.023 


0.026 


o"nM (tot) 


+0.038 
-0.064 


0.036 


+0.036 
-0.062 


0.033 


0.030 


0.034 



SDSS-only 

SDSS+ESSENCE+SNLS 
Nearby+SDSS 

Nearby+SDSS+ESSENCE+SNLS 
<= Nearby+SDSS+ESSENCE+SNLS+HST 
f Nearby+ESSENCE+SNLS 
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for which rcst-framc J7-band is included, and define 
AfinoU = fJ-noU " Z^- If thc mlcs2k2 model is a good de- 
scription of the data, we would expect AfinoU to scatter 
about zero. 

The left panel of Fig. 30 shows the resulting change 
in the SDSS-II SN Hubble diagram, A^noU versus red- 
shift. For z > 0.21, observer-frame g corresponds to rest- 
frame U ; excluding g-band in this redshift range results 
in a shift in the distance modulus of AjinoU = (0.12 ± 
0.02) mag. In the redshift interval 0.21 < z < 0.285, the 
remaining observer-frame SDSS passbands r and i corre- 
spond to rest-frame V and R; for z > 0.285, they corre- 
spond to B and V. In both of those redshift intervals, we 
see a distance modulus shift of A/i„o(7 — 0.1 — 0.15 mag. 
For z < 0.21, the gri filters do not map onto rest-frame 
C/-band, therefore AnnoU = in this redshift range. 
Since A^noU is determined from strongly correlated sam- 
ples, i.e., the r- and z-band data are the same in fits with 
and without [/-band, its uncertainty is estimated to be 
ims/^/N, where rms is the root-mean-square and TV is 
the number of SNe la. Applying the same {/-exclusion 
test to simulations shows that the observed shift has a 
significance of 6cr, consistent with our estimate of the 
uncertainty. The large tilt in the Hubble diagram that 
results from excluding g-hand at z > 0.21 results in a 
shift of 6w - 0.3 for the SDSS-only and Nearby+SDSS 
sample combinations, by far the dominant contribution 
to the systematic error budget for those samples. 

We further investigate the J7-band anomaly by com- 
paring the fitted distance moduli with and without rest- 
frame U for three subsamples: (i) observer- frame UBV 
vs. BV for 17 nearby SNe la; (ii) observer-frame ugr vs. 
gr for 9 low-redshift SDSS-H SNe for which the u-band 
signal to noise is sufficient (z < 0.1), and (iii) observer- 
frame gri vs. ri for 13 SNLS SNe with 0.21 < z < 0.5 and 
that have at least one g-haiid measurement within ±10 
days of maximum brightness. In each case, the subsam- 
ples are chosen such that the light curves pass the selec- 
tion criteria of § 4 and include three observer passbands, 
one of which maps onto rest-frame [/-band. The dif- 
ferences in distance modulus (Aiinou) between the two- 
and three-band fits (without and with rest-frame U) are 
shown in the right panel of Fig. 30. The SDSS-II & SNLS 
subsamples show a consistent shift of about 0.1 mag. For 

TABLE 13 

For the ACDM model, constraints on Qm and Qa from 

MLCS2k2 SN DISTANCES COMBINED WITH SDSS BAO AND 

WMAP5 CMB RESULTS. 
Result for sample combination: 
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b 


c 


d 


e 


/ 




55.2 


171.7 


87.1 


203.9 


237.9 


145.8 


^dof 


102 


220 


135 


253 


287 


150 


RMS^ 


0.15 


0.20 


0.16 


0.20 


0.21 


0.23 


Qa 


0.735 


0.714 


0.735 


0.713 


0.705 


0.718 


atij^ (stat) 


0.019 


0.019 


0.019 


0.019 


0.018 


0.019 


o-cij^ (syst) 


0.000 


0.000 


0.006 


0.004 


0.004 


0.004 


CflA (*°*) 


0.019 


0.019 


0.019 


0.019 


0.018 


0.019 


flu 


0.274 


0.300 


0.274 


0.302 


0.312 


0.294 


o-Om (f^tat) 


0.023 


0.023 


0.023 


0.023 


0.022 


0.023 




0.000 


0.000 


0.001 


0.001 


0.001 


0.001 


<^nM (tot) 


0.023 


0.023 


0.023 


0.023 


0.022 


0.024 



comparison, the right panel of Fig. 30 also shows the av- 
erage shift for the points in the left-panel test at z > 0.21, 
again showing consistency. By contrast, the nearby sub- 
sample is consistent with no shift or a slightly negative 
shift. Since the MLCS2k2 model is trained on a superset 
of the nearby SN la data, we would expect no signifi- 
cant shift for the nearby subsample. The redshift range 
z < 0.1 in the left panel of Fig. 30 is not the same as the 
SDSS-wijr (z < 0.1) test in the right panel: the former 
is based on observer-frame gri and does not map onto 
rest-frame /7-band, while the latter is based on ugr vs. 
gr in order to test excluding rest-frame t/-band. 

Since the [/-band is particularly sensitive to host- 
galaxy extinction, it is worth exploring whether this 
anomaly might be an artifact of the assumed extinc- 
tion law. We have repeated the test above, replacing 
the CCM89 color law nominally used in MLCS2k2 with 
the empirically determined SALT-II color law, CL{\) in 
Eq. 6. The SALT-ii color law results in a [/-band ex- 
tinction that is 0.07 mag larger than that from using 
CCM89 with Ry = 2.18 and the mean extinction value. 
Ay = 0.35. The differences in the extinction for the other 
passbands are much smaller: 0.014, 0, 0.007, -0.002 for 
B,V,R,I. Repeating the [/-band exclusion tests with 
the SALT-II color law results in distance-modulus off- 
sets for the subsamples in the right panel of Fig. 30 that 
are ~ 20% smaller than for the nominal test using the 
CCM89 color law. 

Although excluding rest-frame [/-band reveals a prob- 
lem, it does not definitively indicate that rest-frame [/, 
as opposed to one of the other passbands, is the source of 
the problem. To study this in more detail, we carry out 
a similar test in which observations in passbands corre- 
sponding to rest-frame B-hand are excluded. For SDSS- 
II SNe with redshifts z < 0.21, this B-exclusion test cor- 
responds to comparing distance moduli from observer- 
frame gri (rest-frame BV R) with those from just r and i 
(rest- frame VR). The difference in the average distance 
modulus, nvR — I^BVR, is —0.01 ± 0.02 mag, consistent 
with no shift. This test suggests that rest-frame B is not 
the source of the anomaly and strengthens the circum- 
stantial evidence that rest-frame U is the source of the 
anomaly. 

To further diagnose the [/-band anomaly, we compare 
the light-curve data and the MLCS2k2 model light-curve 
fits as a function of epoch for the different rest-frame 
passbands and the five different SN samples. Figure 5 
shows the data- model residuals for the nominal fits. Fig- 
ure 31 shows the residuals when the filter correspond- 
ing to rest-frame [/-band is excluded from the fit; to see 
the [/-band residuals, we must use the nominal mlcs2k2 
model parameters that include [/-band in the training. 
When [/-band is excluded for the nearby SNe la, there 
is a negligible change in the [/-band residuals; this is ex- 
pected because mlcs2k2 is trained on the nearby data. 
For the SDSS-II, excluding [/-band from the fits results 
in a ~ 0.05 mag shift in the [/-band residuals near the 
time of peak brightness (T^st = 0). The ESSENCE sur- 
vey used only two passbands, R and /, and only a few of 
their SNe la probe rest-frame [/-band, so we cannot use 
this sample to probe the problem. For the SNLS sam- 
ple, the residuals have been plotted for those SNe la that 
have observer-frame (7-band measurements mapping into 
rest-frame [/-band (most of the SNe with z < 0.48); i.e.. 



36 



the same subset used in Fig. 30. When rest-frame U- 
band is excluded, the corresponding shift in the U-hand 
residuals is consistent with the shift seen in the SDSS-II. 
We have also examined the higher-redshift SNLS subset 
for which observer-frame r-band maps onto rest-frame 
i7-band; the results of those tests are consistent with 
the tests based on gi-band, but with much larger uncer- 
tainties. Note that the SNLS subset with excluded g- 
band has a spectroscopic efhciency Cspoc = 1 (see §6.2), 
while the SDSS-II spectroscopic efficiency is significantly 
smaller than 1. Although the corresponding MLCS2k2 
priors are quite different, the [/-band anomaly is consis- 
tent between them, indicating that it is not caused by 
errors in espoc- 

Since the mlcs2k2 model is trained with the nearby 
sample, the U-hand anomaly in the other SN la samples 
results in significant systematic uncertainties that limit 
the precision of cosmological parameters obtained with 
the current implementation of the mlcs2k2 method. 
There are a number of possible causes for the t/-band 
discrepancy: (i) the SN la [/-band flux could be redshift 
dependent; (ii) selection effects for the nearby sample 
could result in a [/-band flux distribution that is not 
representative of the true SN la population; (iii) there is 
a problem with the mlcs2k2 model; (iv) the observer- 
frame [/-band flux for the nearby sample is not properly 
translated into the Landolt system; (v) the SN spectral 
energy distribution in the UV region is not adequately 
constrained, leading to errors in the K-corrections. 

The first possibility, a redshift-dependent flux in the 
[/-band, is unlikely, given our test based on the ugr 
passbands for 9 SDSS-II SN la with 0.04 < z < 0.09 
(right panel in Fig. 30). Although this redshift range 
is still slightly higher than that for the nearby sample 
(z ^ 0.03), a very rapid redshift evolution of SN la prop- 
erties would be needed to account for the discrepancy. 
The second possibility is motivated by the very low selec- 
tion efficiency for the nearby sample, as indicated by the 
data-simulation comparison of the redshift distribution 
in the upper-left panel of Fig. 9. However, the nearby 
sample shows good agreement with the other SN la sam- 
ples in the B, V, R passbands, so one would have to pos- 
tulate a selection effect that biases the [/-band more than 
the other bands. The third possibility, of a problem 
with the light-curve model, is difficult to exclude, but the 
residuals in the nearby sample ( "Nearby" column of the 
left panel of Fig. 5) look reasonable; that rules out obvi- 
ous problems, thus narrowing potential problems to the 
extrapolation to higher redshifts and to different pass- 
bands. 

The fourth possibility seems at first sight unlikely, since 
calibration errors are typically quoted at the level of 0.01- 
0.02 mag. but it could be that those errors have been 
underestimated. A mis-calibration of more than 0.1 mag 
in [/-band would be needed to account for the observed 
anomaly. JRK07 report that the light-curve residuals are 
40% larger in [/-band than in the other bands, but that 
does not necessarily point to a calibration problem. The 
[/-band residuals for the nearby sample (upper-left panel 
in Fig. 5) vary with epoch as the SN becomes redder, 
suggesting a problem with the definition of the [/-band 
filter. 

The last possibility, a difference in the UV region of 
the SN spectral energy distribution, has been suggested 



in previous works (Foley et al. 2008; Sullivan et al. 2009). 
Both mlcs2k2 and SALT-ii models assign larger uncer- 
tainties in the UV region compared to the redder bands. 
EUis et al. (2008) found that maximum-light SN la SEDs 
varied significantly more at wavelengths A < 4000 A than 
for redder wavelengths and concluded that the additional 
dispersion is not due to extinction from host-galaxy dust. 

Since the U -band anomaly appears in multiple samples 
(SDSS-II and SNLS) as well as in multiple redshift ranges 
for the SDSS-II sample {z < 0.1 and z > 0.21 in right 
panel of Fig. 30) we believe that the problem lies within 
the nearby SN sample, i.e., with observations in observer- 
frame [/-band. The most likely source of the problem is 
either the training procedure or the translation into the 
Landolt system. We will address this issue in the future 
by retraining light-curve models with SDSS-II SNe; the 
modeling of rest-frame [/-band will make use of u-band 
measurements at z < 0.1 and g-band measurements at 
z > 0.21. 
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Fig. 30. — Left panel: redshift dependence of the average dif- 
ference in /I for SE)SS-II SNe between the nominal mlcs2k2 fits 
and fits in which the observer-frame passband corresponding to 
rest-frame [/-band is excluded (g-band for z > 0.21). Labels on 
the plot indicate the corresponding rest-frame UBV (Bessell 1990) 
passbands. Error bars (rms/\/7V) refiect the statistical uncertainty 
on the mean ^-difference in each redshift bin. Right panel: shift 
in average distance modulus when rest-frame (/-band is excluded 
from the mlcs2k2 fits, for the subsamples discussed in the text. 



10.1.4. Comparison with WV07 Result 

The Nearby-hESSENCE-hSNLS sample combination 
(f) corresponds to one of the combinations analyzed by 
the ESSENCE collaboration in WV07 (with the excep- 
tions of a different minimum redshift cut and different 
light-curve selection criteria) . To compare with those re- 
sults, we repeat the sample (/) analysis using the same 
BAO prior as in WV07, i.e., without the CMB prior. We 
flnd w = -0.75±0.11(stat), which differs from the WV07 
value by 0.32, or about idstat. Since the two analyses 
are based on the same data, the statistical significance of 
the discrepancy is much larger than ScJstat ■ If we add the 
systematic uncertainties in quadrature (Swsyst = 0.18), 
which is clearly an overestimate, the discrepancy is still 
fairly significant (1.8cr). Both analyses are based on the 
mlcs2k2 method, but our changes in mlcs2k2 param- 
eters and priors result in systematic differences that are 
explained below. 

We run the mlcs2k2 fitter in "WV07 mode" (§ 5.1), 
using the WV07 parameter choices and reproducing their 
result for w, and then make cumulative changes in the 
analysis that evolve it toward the parameter inputs and 
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Fig. 31. — Data-model fractional residuals as a function of rest- 
frame epoch in 5-day bins for mlcs2k2 light-curve fits. Same as 
Fig. 5, except fits exclude observer-frame filter corresponding to 
rest-frame fZ-band 

light curve fitting co(ic used in our fiducial analysis. The 
resulting shifts in w relative to the WV07 value are 
shown in Fig. 32. The changes from WV07 that we im- 
plement sequentially are: change the minimum redshift 
for the nearby sample from 0.015 to 0.02 {z > 0.02); 
set off-diagonal model covariances to zero (off-diag); im- 
plement K-correction improvements in item 2 of § 5.1 
(Kcor); replace the Bessell filter shifts introduced by 
Astier et al. (2006) and adopted by WV07 to the color 
transformation method used in our analysis and simul- 
taneously change the primary standard from Vega to 
BD-l-17 (Calibration); use mlcs2k2 model parameters 
M, p, q (Eq. 1) corresponding to Ry = 2.18 (§ 7.2), but 
still use Rv = 3.1 and WV07 priors in the light curve fits; 
use Rv = 2.18 with WV07 priors in the hght curve fits; 
use the Ay prior and efficiency from this analysis and re- 
move the WV07 requirement that each observation has 
SNR> 5 (prior). 

The largest source of change in w (^ 0.25) results from 
our different assumptions about host-galaxy dust and the 
mlcs2k2 fitting prior. Our fitting prior is based on a 
measurement of Ry and of the Ay distribution using 
the SDSS dust sample (§ 7) along with a comprehensive 
model of survey efficiencies using Monte Carlo simula- 
tions. In contrast to our analysis, the fitting prior used 
in WV07 is based on the assumption that Ry =3.1, that 
the Ay distribution is represented by the "galactic line- 
of-sight" (glos) distribution, and that the spectroscopic 
targeting efficiency is unity. 

10.2. Results with SALT-II 

Using the SALT-ii method (§8.2), we present cosmolog- 
ical results for the six sample combinations (a-j) of Table 
4. Table 14 gives the spectroscopic redshift and derived 
SALT-II fit parameters for each supernova that passes the 
selection cuts of §4. Recall that the fit parameters xi and 
c are estimated from the light-curve fits for each object. 
For each sample combination, those fit parameters are 
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Fig. 32. — Changes in w between the WV07 analysis and the 
analysis presented here, for the Nearby+ESSENCE-|-SNLS sample 
combination. Changes are implemented sequentially from left to 
right, and the cumulative differences in w are shown. The changes 
are described in the text. 

used to estimate the cosmological parameters, the global 
parameters a and (3 in Eq. 29 (see Tables 16 and 17), 
and the distance moduli, by minimizing the scatter in 
the Hubble diagram. Since the distance modulus esti- 
mate for each SN depends upon the sample combination 
in which the SN is included, upon the cosmological model 
parametrization, and upon the BAO and CMB priors, we 
provide tables for each of the six sample combinations 
and for both the FwCDM and ACDM models. Ahhough 
the cosmological parameters have been corrected for the 
selection bias using simulations (Table 5), the distance 
moduli in Table 14 do not include any bias corrections. 
The entries in this Table should therefore not be used to 
derive cosmological constraints. 



TABLE 14 

SALT-II LIGHT CURVE FIT PARAMETERS INCLUDING BAO + CMB 
PRIORS, (uncertainties IN PARENTHESES). THE COMPLETE SET 
OF TABLES FOR EACH SAMPLE-COMBINATION (a-f) IS GIVEN IN 
ELECTRONIC FORM IN THE JOURNAL, AND ALSO AT 
http : //das . sdss . org/va/SNcosmology/sncosm09_f its . tar . gz. 



SNID 


redshift^ 










762 


0.1904(.0000) 


39.91(0.09) 


-0.01(0.03) 


0.73(0.33) 


53625.2(0.4) 


1032 


0.1291(.0000) 


38.75(0.17) 


0.14(0.06) 


-3.09(0.39) 


53626.6(0.4) 


1112 


0.2565(.0000) 


40.53(0.16) 


0.00(0.04) 


-1.06(0.69) 


53630.2(0.7) 


1166 


0.3813(.0000) 


41.23(0.25) 


0.03(0.07) 


1.10(1.12) 


53631.9(1.2) 


1241 


0.0858(.0000) 


37.91(0.09) 


0.09(0.03) 


-0.90(0.17) 


53635.3(0.2) 



''Spectroscopic redshift in CMB frame. 



10.2.1. SALT-II Hubble Dispersion 

In the SALT-II method, an intrinsic dispersion ((tJ"*) is 
added in quadrature to the distance modulus uncertainty 
(Eq. 27) such that the resulting Hubble diagram x^/^dof 
is equal to one. Using the FwCDM model parametriza- 
tion and the BAO-I-CMB prior. Table 15 gives the crjf* 
values obtained from fitting each SN sample indepen- 
dently and setting xfi = -^dof- For the nearby, SDSS-II, 
and ESSENCE samples, the SALT-II values are similar to 
those from mlcs2k2 in the fourth Hne of Table 11. For 
the SNLS sample, the SALT-ii dispersion is significantly 
smaller than that from mlcs2k2, while for the HST sam- 
ple the SALT-II dispersion is larger. The smaller disper- 
sion for SNLS may derive in part from the fact that the 
SALT-II model was partially trained on SNLS data. For 
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sample combinations a-/, the a™^ values are in Tables 16 
and 17. 



TABLE 15 

Intrinsic dispersion (o-jf*) required for = Af^of for each 

SN lA SAMPLE FIT SEPARATELY WITH THE SALT-II METHOD. 





Nearby 


SDSS-II 


for sample: 
ESSENCE 


SNLS 


HST 


independent fits 


0.15 


0.08 


0.17 


0.11 


0.23 



10.2.2. SALT-II Hubble Diagrams and Cosmological 
Parameters 

Figure 33 shows the dilferences between the SALT~ii 
estimated SN distance moduli /i^ and those for an open 
CDM model with no dark energy {Q.m = 0.3, Ode = 0) as 
a function of redshift. Figure 34 shows the distribution of 
normalized residuals, (^i — /^FmCDM)/^'^!, where /xfu,cdm 
is the distance modulus from the best-fit FiuCDM model 
for sample combination (e), and cr^ is the total uncer- 
tainty defined in Eq. 27. The bulk of the distribution of 
all 288 normalized residuals (upper left panel of Fig. 34) 
is well fit by a Gaussian with a = 0.75; outliers increase 
the rms to 0.92. The rms for combination (e) is slightly 
smaller than one because covariances from the fit are not 
included in the calculation of a,, . 
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Fig. 33. — Hubble residuals for the SALT-II method: differences 
between measured SN distance moduli and those for an open CDM 
model {Qm = 0.3, Ode = 0) vs. redshift for the six SN sample 
combinations. Large, square (pink) points show weighted averages 
in weighted redshift bins (see Fig. 24 caption). Solid curves show 
the difference between the best-fit FuiCDM model distance mod- 
ulus for that sample combination and that for the open model, 
normalized to have the same value of the Hubble parameter. The 
error bars on the data points correspond to the distance modulus 
error tr^* from the SALT-II light-curve fit (Eq. 27), i.e., they do not 
include the intrinsic scatter or the effects of redshift and peculiar 
velocity errors. The vertical bars show the values of the intrinsic 
uncertainty, crjf * , included in the cosmology fits so that the P^i' 
degree of freedom is unity for each Hubble diagram. 
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Fig. 34. — Distribution of normalized Hubble residuals (pull) for 
the SALT-II method, for sample combination (e) comprising all five 
samples (upper left, along with Gaussian fit), and for each SN sam- 
ple indicated on the panels, fi is the measured SN distance mod- 
ulus, ^FmCDM is the distance modulus from the best-fit FuiCDM 
model for sample combination (e), and cr^ is the total uncertainty 
(Eq. 27). 

Figures 35 and 36 show the SALT-ii statistical- 
uncertainty contours for the FwCDM and ACDM 
models; for the latter, the SALT-ii SN contours are 
more consistent with the BAO-I-CMB constraints than 
the MLCS2k2 contours were. For the combined 
SN-I-BAO-I-CMB results, the total uncertainty contours, 
including systematic errors, are shown in Figs. 37 and 
38. The best-fit cosmological parameter values and un- 
certainties, marginalizing over Hq and incorporating the 
bias corrections of Table 5, are given in Tables 16 and 17. 
The SALT-II statistical errors on cosmological parameters 
are consistent with those from mlcs2k2. The system- 
atic errors for the two methods are similar for sample 
combinations (6) and {d-f), but the SALT-ii systematic 
uncertainty is significantly smaller for the sample combi- 
nations in which SDSS-II is the high-redshift sample (a 
and c) . The large difference in the systematic uncertainty 
is driven by the f/-band anonaly. 

Among the six SN sample combinations, the best-fit 
values of w for the FwCDM model again fall into two 
groups, as for the mlcs2k2 results, but with the op- 
posite trend in w compared to mlcs2k2. For sample 
combinations (a) and (c), in which SDSS-II is the high- 
redshift sample, the SALT~ii method results in w = —0.87 
and —0.92, comparable to the MLCS2k2 results of —0.84 
and —0.92 for the same sample combinations. However, 
for the other four sample combinations, which include 
higher-redshift SNe, SALT-ii yields w = -0.98 to -0.95, 
compared with the mlcs2k2 result oi w = —0.71 to 
—0.76 for the same sample combinations. Based on stud- 
ies with simulated sample combinations, the observed dif- 
ference between these two groups of SALT-ii results is not 
statistically significant. However, the difference between 
the mlcs2k2 and SALT-ii results for the higher-redshift 
samples appears to be significant. We discuss these dif- 
ferences further in § 11. 

10.2.3. Redshift Evolution of SALT-ii Parameters 

In the SALT-II model fits, a, /3, and M (see Eq. 29) are 
global parameters that are assumed to be independent of 



TABLE 16 

For the FtoCDM model, constraints on w and Qm from 

SALT-II SN DISTANCES COMBINED WITH SDSS LRG BAO AND 

WMAP5 CMB RESULTS. 



Result for sample: 





a. 


b 


c 


d 


e 


f 




0.084 


0.124 


0.105 


0.128 


0.140 


0.160 


RMS^ 


0.178 


0.219 


0.170 


0.210 


0.232 


0.231 


w 


-0.87 


-0.98 


-0.92 


-0.98 


-0.96 


-0.95 


(Tuj (stat) 


0.12 


0.08 


0.11 


0.07 


0.06 


0.08 


(Tiu (syst) 


+0.06 

n 1 o 
— U. iZ 


0.15 


+0.07 

— U. 10 


0.13 


0.12 


0.13 


(Tm (tot) 


+0.14 
-0.17 


0.17 


+0.13 
-0.18 


0.15 


0.14 


0.15 




0.281 


0.256 


0.271 


0.264 


0.265 


0.267 


0-S2M (stat) 


0.030 


0.019 


0.025 


0.017 


0.016 


0.019 




+0.015 
-0.025 


0.033 


+0.015 
-0.029 


0.028 


0.025 


0.027 




+0.034 
-0.039 


0.038 


+0.029 
-0.038 


0.033 


0.030 


0.033 


a 


0.127 


0.123 


0.113 


0.107 


0.124 


0.106 


ctq [stat) 


0.017 


0.015 


0.014 


0.013 


0.014 


0.019 


<Ta{syst) 


0.020 


0.021 


0.016 


0.020 


0.023 


0.023 


CTq {tot) 


0.026 


0.026 


0.021 


0.024 


0.027 


0.030 


p 


2.52 


2.62 


2.50 


2.66 


2.64 


2.56 


ap(stat) 


0.16 


0.13 


0.15 


0.12 


0.12 


0.17 


a p (syst) 


0.11 


0.19 


0.11 


0.19 


0.18 


0.25 


o-p{tot) 


0.19 


0.23 


0.19 


0.22 


0.22 


0.30 



" SDSS-only 
SDSS+ESSENCE+SNLS 
Nearby+SDSS 

Nearby+SDSS+ESSENCE+SNLS 
" Nearby+SDSS+ESSENCE+SNLS+HST 
f Nearby+ESSENCE+SNLS 



TABLE 17 

For the ACDM model, constraints on S7m and Qa from 
salt-ii sn distances combined with bao and cmb results. 







Result for sample: 








a 


b 


c 


d 


e 


/ 


int 

M 


0.085 


0.123 


0.105 


0.128 


0.140 


0.160 


RMSp 


0.177 


0.220 


0.170 


0.210 


0.232 


0.231 




0.734 


0.735 


0.734 


0.734 


0.727 


0.734 


o-Qj^ (stat) 


0.019 


0.017 


0.019 


0.017 


0.016 


0.017 




0.019 


0.026 


0.018 


0.020 


0.019 


0.018 


un^ {tot) 


0.027 


0.031 


0.026 


0.026 


0.025 


0.025 




0.275 


0.274 


0.275 


0.275 


0.279 


0.275 




0.023 


0.021 


0.023 


0.020 


0.019 


0.021 




0.014 


0.021 


0.013 


0.020 


0.017 


0.016 


{tot) 


0.027 


0.030 


0.027 


0.029 


0.026 


0.027 


a 


0.126 


0.123 


0.113 


0.113 


0.116 


0.105 


(Ta{stat) 


0.017 


0.015 


0.014 


0.013 


0.014 


0.019 


aa{syst) 


0.020 


0.021 


0.016 


0.019 


0.023 


0.023 


<Ja {tot) 


0.027 


0.026 


0.022 


0.023 


0.027 


0.030 





2.58 


2.64 


2.50 


2.57 


2.65 


2.46 


ap{stat) 


0.15 


0.12 


0.15 


0.12 


0.12 


0.17 


af3{syst) 


0.13 


0.18 


0.09 


0.17 


0.18 


0.19 


ag{tot) 


0.20 


0.22 


0.17 


0.20 


0.22 


0.25 



39 




0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 1.0 



Fig. 35. — For the FuiCDM model, SALT-ii statistical-uncertainty 
contours in the Q^-^ plane for each of the six SN sample combi- 
nations indicated on the plots. Long, black contours: 68%, 95%, 
and 99% confidence level regions for the SN data alone; green con- 
tours: corresponding CL regions for SDSS BAO (Eisenstein et al. 
2005); blue contours: CL regions for WMAP5 CMB (Komatsu 
et al. 2008); closed, black contours: combined constraints from 
SN+BAO+CMB. 

redshift. To test the consistency of this assumption, we 
have carried out SALT-II fits separately in five redshift 
bins for sample combinations (d) and (e). For the fit 
in each redshift bin, the cosmological parameters w and 
flu in the FwCDM model are fixed to the values from 
the sample (e) fit, and the BAO+CMB prior is applied. 
The redshift dependence of the best-fit results for a, f3, 
and M are shown in the left panels of Fig. 39. There is 
evidence for redshift evolution of the parameters, partic- 
ularly for the color parameter /3, which falls with increas- 
ing redshift above z ~ 0.6. This trend is more evident 
without the HST sample, suggesting that the /3- variation 
is driven primarily by the SNLS sample. We have per- 
formed this test with high-statistics simulations (right 
panels of Fig. 39) and find that the fitted SALT-ii pa- 
rameters are consistent across all redshift bins. Since 
the simulation accounts for Malmquist bias, selection ef- 
fects, and measurement errors, the redshift dependence 
favored by the data is likely due to some other effect. 

10.2.4. U-Band Anomaly with SALT-ii 

Here we study how the t/-band anomaly is manifest 
in the SALT-ii method. Fig. 40 shows the change in the 
binned SDSS-II SN Hubble diagram when the SALT-ii 
fit excludes data from the observer-frame passband cor- 
responding to rest-frame [/-band, i.e., excluding g-band 
data for z > 0.21. Although this change only affects the 
light-curve fits for z > 0.21, it can alter the estimated 
distance moduli at all redshifts since they arc derived 
from a global fit to the Hubble diagram. For z < 0.21, 
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Fig. 36. — For the ACDM model, SALT-ii statistical-uncertainty 
contours in the Cm-^^A plane for each of the six SN sample com- 
binations indicated on the plots. Long, black contours: 68%, 95%, 
and 99% confidence level regions for the SN data alone; green 
contours: corresponding CL regions for SDSS BAO (Eisenstein 
et al. 2005); blue contours: CL regions for WMAP5 CMB (Ko- 
matsu et al. 2008); closed, red contours: combined constraints from 
SN-I-BAO+CMB. 
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Fig. 37.— For salt-ii and the FwCDM model, 68% CL contours 
in the ^y^-w plane for each of the six SN sample combinations, 
using the combined SN-I-BAO+CMB constraints. Solid contours 
are total (statistical-l-systematic) uncertainty; dashed contours arc 
statistical only. Systematic errors have been included using the 
prescription in Appendix F. 

the average distance modulus shift is (0.012±0.003) mag; 
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Fig. 38. — For salt-ii and the ACDM model, 68% CL contours 
in the CM-f^A plane for each of the six SN sample combinations, 
using the combined SN-I-BAO+CMB constraints. Solid contours 
are total (statistical+systcmatic) uncertainty; dashed contours are 
statistical only. Systematic errors have been included using the 
prescription in Appendix F. 

for z > 0.21 the mean shift is (0.079 ± 0.028) mag rel- 
ative to the z < 0.21 shift. This relative shift is about 
half of that for mlcs2k2, but it is still significant: for the 
SDSS-only sample, the exclusion of [/-band in the salt- 
ii fits results in a shift in w of —0.100 that is included as 
a systematic uncertainty. 

The data-model residuals for the SALT-ii light-curve 
fits are shown in Fig. 7 for all of the SN la samples 
in all of the rest-frame passbands. Fig. 41 shows the 
residuals when the filter corresponding to rest-frame U- 
band is excluded from the fit. In both cases, there 
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Fig. 39. — salt-ii fitted values and uncertainties for a, 
f}, and M, evaluated independently in redshift bins for the 
Nearby+SDSS+ESSENCE+SNLS+HST sample combination (e) 
(solid dots) and for the Nearby+SDSS+ESSENCE+SNLS combi- 
nation (d) (open circles). Left is for data; right is for a simulation 
with 10 times the number of SNe as in the data. For each redshift 
bin, the cosmological parameters w and SIm are fixed to the val- 
ues from the global combination (d or e) fit. The redshift bins are 

= 0.2 for z < 0.8; the highest-redshift bin includes all SNe la 
with z > 0.8. The dashed lines show the ±1(T statistical error band 
based on the global fit to combination (e). 
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are systematic discrepancies between the model and the 
data in rest-frame L/-band for the Nearby SN la sam- 
ple at all epochs. Since the SALT-ii rest-frame L/-band 
model was trained primarily on higher-rcdshift SNe, 
i.e., downweighting observer- frame f/-band data from 
nearby SNe, this points to a systematic offset between 
the nearby, observer-frame U-hand data and the higher- 
redshift, rest-frame U-hand data, which was also quali- 
tatively seen in the mlcs2k2 fits. 
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Fig. 40. — Redshift dependence of the average difference in dis- 
tance modulus for SDSS-II SNe between the nominal SALT-II light- 
curve fits and fits in which the observer-frame passband corre- 
sponding to rest-frame (7-band is excluded (g-band for z > 0.21). 
Labels on the plot indicate the corresponding rest-frame UBVR 
(Bessell 1990) passbands. Error bars (rms/v^) reflect the statis- 
tical uncertainty on the mean /^-difference in each redshift bin. 
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Fig. 41. — Data-model fractional residuals as a function of rest- 
frame epoch in 5-day bins for SALT-II light-curve fits. Same as 
Fig. 7, except fits exclude observer-frame filter corresponding to 
rest-frame (7-band. 



11. COMPARISON BETWEEN mlcs2k2 AND salt-ii 
RESULTS 

As noted in §10.2, the best-fit values of w for the 
FwCDM model agree between SALT-ii and mlcs2k2 to 
within 0.04 for the SDSS-only and Nearby-^SDSS sam- 
ple combinations. However, as indicated in Tables 12 



and 16, when the ESSENCE and SNLS data are included 
(sample combinations 6,(i,e, and/), the u)-valucs increase 
by ^ 0.1 for the mlcs2k2 method, while decreasing by 
nearly 0.1 for the SALT-II method, leading to a discrep- 
ancy of Aw ~ 0.2 between the two methods. To esti- 
mate the statistical significance of this discrepancy, we 
have run mlcs2k2 and SALT-ii fits on a set of ten simu- 
lated sample combinations generated with the mlcs2k2 
model, each with the same statistics as the data. The 
predicted rms-spread on Aw for combination (e) (all five 
samples) is Awstat — 0.05, so the observed discrepancy 
appears to be statistically significant. 

To help diagnose this discrepancy, we compare the 
mlcs2k2 and SALT-ii fitted parameters in Fig. 42 for 
each of the SN samples. In this figure, the SALT-ii pa- 
rameters are based on the fit to sample combination (e), 
including the BAO-I-CMB prior, and the SALT-ii dis- 
tances have not been corrected for selection bias. The 
left panels show the mean difference in distance modu- 
lus, A^ = HSAUT2 — Mmlcs, as a function of redshift. 
The values of the Hubble parameter have been relatively 
adjusted so that the values of fi for the best-fit mlcs2k2 
and SALT-II models for FwCDM agree at z ^ 0. The 
two methods yield consistent distance estimates for the 
Nearby, SDSS-II, and ESSENCE samples over the red- 
shift ranges they cover. Moreover, the scatter in Afi 
is comparable to the intrinsic scatter: ^(A/i) = 0.10, 
0.15 and 0.16 mag for these three samples, respectively. 
Analyzing the Nearby-hSDSS+ESSENCE combination 
(not one of our standard combinations), w = —0.84 for 
mlcs2k2 and w = —0.89 for SALT-ii, indicating good 
agreement. The situation changes dramatically when we 
include the SNLS data, for which there is a clear trend 
of increasing Afi with redshift: this is the primary cause 
of the difference in w between MLCS2k2 and SALT-ii for 
sample combinations b, d, e, and /. 

The middle and right panels of Fig. 42 show the cor- 
relations between the mlcs2k2 and SALT-ii light-curve 
fit parameters, SALT-ii color c versus mlcs2k2 extinc- 
tion Ay and SALT-ii versus mlcs2k2 shape-luminosity 
parameters xi versus A. The fitted slopes, dc/dAy and 
dxi/dA, are consistent among the SN samples, except 
for ESSENCE which has somewhat larger slopes. The 
latter could be related to the fact that this sample has 
only two observer-frame passbands and hence less reli- 
able color information. 

The final light-curve fit parameter to compare is 
the epoch of maximum light in rest-frame i3-band 
(to)- Fig- 43 shows the distributions of to(SALT- 
ii)— io(MLCS2K2) for the five SN samples. On average, 
the SALT-ii epoch of peak light is about 1 — 1.5 days later 
than that for MLCS2k2, with a dispersion of about one 
day for the ground-based samples. For the HST sample, 
the offset is ~ 2 days and the dispersion is larger (3 days) 
because of a handful of SNe with large to-discrepancies. 
Fitting simulated SN samples shows that for the poorly 
sampled HST light curves, the fitted to from SALT-ii is 
expected to be more discrepant from the input value than 
for mlcs2k2. 

To gain some insight into the discrepancy between the 
two fitting methods, wc have investigated modifications 
to the mlcs2k2 model that are designed to partially 
replicate the features of the SALT-ii model. The mod- 
ifications are: (i) change the mlcs2k2 model vectors so 
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Fig. 42. — Comparisons between mlcs2k2 and SALT-II light-curve 
fit parameters for the five SN samples: difference in distance mod- 
ulus vs. redshift (left), color c vs. Ay (middle), and xi vs. A 
(right). Crosses in left panels show average & uncertainty in red- 
shift bins. The solid straight linos (middle & right columns) are 
the same within a column to guide the eye, and they are derived 
from a fit to all of the samples except for ESSENCE, as explained 
in the text. 
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Fig. 43. — Comparison of SALT-II and mlcs2k2 fitted epoch of 
maximum B-band light for the five SN samples: shown are the 
distributions of fo(sALT-ll)— to(MLCS2K2). Insets show mean and 
rms for each sample. 

that the UBVRI hght-curve templates match those of 
the synthesized SALT-ii model, as described in Appendix 
G; (ii) use a flat prior. Ppi.ior(-z, ^v, A) = 1 in Eq. 5, al- 
lowing negative values of Ay; (iii) use the SALT-II color 
law CL{\) in place of the CCM89 model in the extinction 
term in Eq. 1; and (iv) exclude measurements for which 
the mean filter wavelength exceeds 7000 ■ {1 + z) A. Us- 
ing this modified mlcs2k2 model to fit the light curves, 
the resulting values of w decrease by 0.2 for sample com- 
binations (d) and (e) and are in good agreement with 



the SALT-II results. Fig. 44 compares the SALT-ii pa- 
rameters to those from the nominal mlcs2k2 fits and to 
those from the modified mlcs2k2 model. For the latter, 
the distance moduli show good agreement over the entire 
redshift range, and the correlations between the shape- 
color parameters (c vs. Ay and xi vs. A) are stronger 
than for the nominal mlcs2k2 fits. 

Of the four modifications to MLCS2k2, the last two 
(use of the SALT-ii color law and a 7000 A cutoff) 
have a negligible effect on the cosmology analysis. The 
change in w is mainly due to matching the template 
light curves and to using a flat Ay prior. Implement- 
ing either of these changes alone results in a w-shift of 
—0.07 or less; both changes are needed to obtain the 
full w-shift of —0.2. This exercise shows that the w- 
discrepancy between MLCS2k2 and SALT-ii is related in 
part to the training procedure that determines the spec- 
tral and light-curve templates and in part to the different 
assumptions about SN color variations. The model V- 
band light curve and color evolution for mlcs2k2 and 
modified MLCS2k2 are compared in Fig. 45 for differ- 
ent values of the shape- luminosity parameter A. The 
main discrepancy is in [/ — B, which differs by about 
0.1-0.2 mag between mlcs2k2 and modified mlcs2k2. 
Since the modified mlcs2k2 colors agree well with the 
SALT-II colors (not shown), the U — B discrepancy in 
Fig. 45 can be interpreted as the discrepancy between 
the mlcs2k2 and SALT-ii models. In particular, since 
the models agree in i3-band, this plot illustrates the dif- 
ferences in rest-frame [/-band. 
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Fig. 44. — For all five samples (combination e), comparison of 
SALT-II fit parameters with those from the nominal mlcs2k2 fits 
(left) and from the modified mlcs2k2 fits (right) as described in 
the text. 



As noted earlier, there is strong evidence of systematic 
discrepancies in rest-frame [/-band between the nearby 
and higher- redshift samples. These discrepancies are 
reflected in the differences between the mlcs2k2 and 
SALT-II [/-band models, differences that account for part 
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A = +0.4 
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Fig. 45. — For mlcs2k2 model, rest-frame V-hand magnitude 
(top row) & colors (lower 3 rows) versus epoch. Solid curves are 
for the nominal mlcs2k2 model, and dashed arc for the modified 
mlcs2k2 model. The latter was constructed to reproduce the SALT- 
II flux and color evolution. Each column corresponds to a different 
value of mlcs2k2 parameter A = -0.4,0,0.4. More than 70% of 
SNe la have A < 0, and 5% have A > 0.4. 

of the cosmological parameter disagreement between the 
two models. The other major contributor to the cos- 
mological disagreement is the differing treatment of SN 
color variation in the two models. There is a trend to- 
ward negative apparent SALT-ii color at high-redshift 
within the SNLS sample. SALT-II and mlcs2k2 with 
a flat- Ay prior assign these blue events large intrinsic 
luminosities and therefore large distance moduli. By 
contrast, mlcs2k2 with the nominal Ay prior identifies 
these events as having Ay ^ and assigns them lower 
luminosities and distances. As illustrated in Fig. 19, the 
nominal MLCS2k2 interpretation of these events is con- 
sistent with the observed color distributions, so it is not 
obvious which model is correct. 

In § 10 we have considered the approach of exclud- 
ing data corresponding to rest-frame [/-band altogether. 
In that case, for sample combinations that include 
ESSENCE and SNLS (6, d, e, /), the salt-ii - mlcs2k2 
discrepancy of Aw ^ 0.2 is reduced by nearly 50%. By 
contrast, the same test applied to the SDSS-only sample 
increases Aw from 0.06 to 0.2 when [/-band is excluded. 
Moreover, in the nominal fits (including [/-band), the 
SALT-II and mlcs2k2 Hubble diagrams appear contin- 
uous around z = 0.2, and there is good agreement be- 
tween them across the SDSS-II redshift range. Removing 
[/-band introduces a noticeable step in both Hubble di- 
agrams at z = 0.2, suggesting that both models display 
unphysical behavior when this information is removed 
completely. Such a sharp feature in the Hubble diagram 
is not well captured by a model with constant w, so the 
changes in Aw quoted above should be interpreted with 
extreme caution. 

12. CONCLUSIONS 

We have presented measurements of the Hubble di- 
agram for type la supernovae discovered during Sept.- 



Nov. 2005 by the SDSS-II Supernova Survey and com- 
bined them with other SN la data and with BAO and 
CMB results to derive cosmological constraints. 

For the SDSS-II sample, based on stringent light-curve 
sampling and signal-to-noise criteria, we selected 103 
SNe la with redshifts 0.04 < z < 0.42 for analysis from a 
parent sample of 130 spectroscopically confirmed SNe la. 
In an effort to make the analysis as homogeneous as pos- 
sible among the SN data sets, we have used similar light- 
curve selection criteria and applied the same analysis 
techniques to all samples. We have estimated distances 
for these SNe using both the mlcs2k2 and SALT-ii meth- 
ods, which differ from each other in a number of respects 
(§ 5.3). The analysis includes significant improvements 
in the implementation of the mlcs2k2 method. 

To determine the efficiency functions used in the 
mlcs2k2 priors and the biases in the SALT-ii results, 
we have carried out detailed Monte Carlo simulations 
for each of the SN la data sets, making use of recorded 
observing conditions. The simulation accurately models 
the light curves (fluxes & errors) and the photometric se- 
lection effects; we have also incorporated a quantitative 
model for the spectroscopic selection efficiency based on 
comparing observed and simulated distributions in SN 
redshift and mean extinction. The simulation provides 
an excellent description of the data, as illustrated by the 
data-simulation comparison of the flux distribution for 
each SN sample (Fig. 14). 

Due to selection effects, supernovae in spectroscopic 
samples tend to be brighter and less extinguished, on 
average, than those of the parent population. For SDSS- 
II, we have augmented the spectroscopic sample with 
a larger sample that includes photometrically identified 
SNe la with host-galaxy redshift measurements. We 
used this SDSS-II "dust" sample in conjunction with the 
Monte Carlo simulations to determine the mean host- 
galaxy reddening parameter Ry and the distributions of 
extinction and light-curve shape parameter, P{Av) and 
P(A), for the underlying SN la population. By match- 
ing observed and predicted SN la colors in MLCS2k2, 
we find Rv = 2.18 ± 0.14stat ± 0.48syst, consistent with 
the recent trend toward values lower than the canon- 
ical Milky Way average of 3.1. For comparison, the 
SALT-II analysis of the combined SN la sample (e) yields 
(3 = 2.64±0.12(stat)±0.18 (syst). If aU SN la color varia- 
tion (beyond that associated with light-curve shape) were 
due to dust extinction, we would expect /3 — Ry = 1 . In 
this analysis we find /? — Ry = 0.5 ± 0.5, consistent with 
the expectation for extinction from dust, but also sug- 
gesting that there may be additional sources of SN color 
variation. We find that the underlying Ay distribution 
is well described by an exponential function; this dis- 
tribution is marginally consistent with both the galactic 
line of sight model used by WV07 and the exponential 
distribution found for the nearby sample by JRK07. 

For both the mlcs2k2 and SALT-ii methods, we have 
carried out extensive studies of systematic errors, vary- 
ing a large number of model parameters and assump- 
tions, and we have included the resulting cosmological 
parameter variations in the systematic error budget. For 
each method, the largest source of systematic uncertainty 
is associated with the rest-frame [/-band anomaly. The 
anomaly is manifest in discrepancies between the light- 
curve data and the models, indicating the need for im- 



44 



proved training of the rest-frame J7-band models. When 
data corresponding to rest-frame C/-band are excluded 
from the light-curve fits, the estimated distance modu- 
lus shifts by 0.12 mag for MLCS2k2 and by 0.08 mag for 
SALT-ii. This shift occurs in particular redshift ranges, 
leading to abrupt features in the Hubble diagram; the 
effect is most severe for the SDSS-II SNe, since the U- 
band shift occurs at the median redshift of the sample. 
For this reason, and because dropping t/-band leads to 
larger uncertainties from the significantly reduced color 
constraints, we have chosen not to exclude f7-band for 
our nominal analysis but to include the corresponding 
changes as part of the systematic uncertainty. 

We have combined the SN Hubble diagram with 
BAO and CMB measurements to estimate the eos- 
mological parameters. For the FwCDM model and 
the combined sample of 288 SNe la from all five sur- 
veys, we find w = -0.76 ± 0.07(stat) ± 0.11 (syst), 
VLm = 0.307 ± 0.019(stat) ± 0.023 (syst) using mlcs2k2 
and w = -0.96 ± 0.06(stat) ± 0.12 (syst), = 
0.265 ± 0.016(stat) ± 0.025 (syst) for SALT-ii. This dis- 
crepancy of 0.2 in w between the two analysis meth- 
ods is not due to inclusion of the SDSS-II data: for 
the Nearby-hESSENCE-^SNLS (/) sample combination, 
which excludes SDSS-II, we find the same difference in w 
between the two methods. Our mlcs2k2 results for this 
sample combination (/) differ substantially from those of 
WV07 (§ 10.1.4), who used an earlier version of mlcs2k2 
to analyze a nearly identical sample. We have traced the 
differences primarily to the different priors used, i.e., to 
the different values of i?y, to the different Ay distri- 
butions, and to the inclusion of spectroscopic targeting 
efficiency in the prior. 

We have traced the mlcs2k2 vs. SALT-ii discrepancy 
to the SN model parameters determined in the train- 
ing, particularly the rest-frame [/-band model (Fig. 45), 
and to the different assumptions about the source of 
color variations in mlcs2k2 and SALT-ii. If we restrict 
the analysis to the 136 SNe la in the Nearby-|-SDSS 
combination, we find much better agreement between 
the two analysis methods, but with larger uncertain- 
ties, w = -0.92 ± 0.13(stat)t°;33 (syst) for mlcs2k2 

and w = -0.92±0.11(stat)lo °5 (syst) for SALT-ii. We 
also note that the cosmological parameter uncertainties 
for the SDSS+ESSENCE+SNLS sample combination are 
similar to those for Nearby+ESSENCE+SNLS, i.e., the 
first-season SDSS-II data sample anchors the Hubble di- 
agram with comparable power to the Nearby sample. 

The mlcs2k2 vs. SALT-ii discrepancy for the higher- 
redshift samples raises the question of which method, 
if either, is more reliable. Since the mlcs2k2 [/-band 
model relies entirely on [/-band measurements at low red- 
shift, while SALT-II uses a combination of low and high- 
redshift data, some at redder observer-frame passbands 
than [/, the SALT-ii results may be less biased by the U- 
band anomaly if the problem is related to the observer- 
frame [/-band measurements. On the other hand, the 
SALT-II f3 parameter has a notable redshift dependence 
(Fig. 39) that is inconsistent with the underyling assump- 
tion that it is constant. Concerning the interpretation 
of SN color variations, both models are consistent with 
the data; in particular, a supernova that is bluer than 
the templates could be extremely bright, as indicated 



by SALT-II, or it could be due to measurement uncer- 
tainty and intrinsic color fluctuations, as interpreted by 
mlcs2k2 (Fig. 19). Since we cannot definitively deter- 
mine from the current data that either method is better 
or incorrect, the overall conclusion is that the cosmolog- 
ical parameter w conservatively lies between —1.1 and 
—0.7 . This result reflects the fact that present SN la 
samples have reached the systematic limits of current SN 
model fitters. Although this result is disappointing, we 
are optimistic that this situation is temporary. The full 
three-season SDSS-II sample, with its homogeneous and 
well-modeled selection function, will serve both to an- 
chor the Hubble diagram and to retrain the light-curve 
models. Since the [/-band anomaly is most likely asso- 
ciated with the published nearby SN sample, use of the 
full SDSS-II data as well as updated low-redshift samples 
from the CfA (Hicken et al. 2009a), the Carnegie Super- 
nova Project, and the Nearby Supernova Factory, should 
significantly reduce this major source of systematic un- 
certainty. 

Finally, our use of photometrically identified SNe la to 
measure host-galaxy dust properties is an important step 
toward including these SNe in the Hubble diagram, which 
will increase the statistical power of the data. This will 
be a growing trend in the future, as large surveys, includ- 
ing PanSTARRS, the Dark Energy Survey, and LSST, 
will discover vastly more SNe than can be confirmed with 
available spectroscopic resources. 
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APPENDIX 

A. SPECTRAL WARPING FOR K-CORRECTIONS IN THE mlcs2k2 LIGHT CURVE MODEL 

The MLCS2k2 model makes predictions about rest-frame light curves in U BRVI; K-corrections are used to translate 
the model light curves so that they can be used to fit measurements at non-zero redshift in a variety of passbands. 
This translation requires knowledge of the supernova spectral energy distribution (SED) time sequence. Since the 
SN la population exhibits intrinsic color variations, a single time sequence will not provide an adequate model. The 
standard practice is therefore to "warp" the model SEDs so that they match the colors of the photometric model for 
a given object. Here we describe our procedure for the spectral warping used in determining K-corrections within the 
mlcs2k2 framework (§ 5.1). 

The procedure begins by computing the rest-frame MLCS2k2 model magnitudes for the assumed values of extinction 
(Av), time of peak brightness {to), and shape-luminosity parameter (A). These assumed values are typically deter- 
mined iteratively as part of the ^^-minimization (Eq. 3). For each iteration, an epoch-dependent SN la SED from 
Hsiao ct al. (2007) is warped so that the synthetic colors of the warped SED match the mlcs2k2 model colors at the 
corresponding rest-frame epoch. 

In detail, for each epoch and passband, the rest-frame model magnitude is first computed as the sum of the first four 
terms on the right-hand side of Eq. 1, where the host-galaxy extinction Xhost has been determined from the unwarped 
SED using the values of Ay and Ry- Although the K-correction depends strongly on how the SED is warped, the 
value of Xhost is only weakly dependent on the warping, so Xhost can be determined to good approximation from the 
unwarped SED. As an illustration, when the SED is severely warped to modify the V ~ B color (at peak brightness) 
by 0.5 mag, the estimate of the F-band extinction changes by less than 1%. Next, the SED is warped by multiplying 
it with the CCM89 galactic extinction law, following JRK07. This usage of the CCM89 law to warp is purely a 
mathematical convenience — it has nothing to do with physical extinction. The "Ay-warp" is the value of "Ay" in the 
warp factor for which the synthetic SN la color matches the color of the rest-frame model magnitude computed above. 
The redshift-dcpendcnt K-correction to observer-frame magnitude is then determined from this Ay- warped SED, using 
the redshift and knowledge of the rest-frame and observer-frame passbands. This method is model-independent and 
can therefore be applied to any light-curve model. The SED is warped independently at each epoch and locally in 
wavelength near the passband of interest; we do not do a global fit to match all colors simultaneously. There are 
two potential limitations in this procedure. First, brightness-dependent spectral features are ignored, i.e., we use a 
single, A-independent composite SED at each epoch. Second, for rest-frame U-hand (/-band) there is no constraint 
for warping the SED blueward (redward) of the central wavelength. 

Compared to the treatment in JRK07, we have made a slight improvement to the spectral warping used for K- 
corrections. The rest-frame filter, /', is chosen as the one that covers the equivalent rest-frame wavelength, Arest = 
Aobs/(l + z), where Aobs is the central wavelength of the observed passband /. For spectral warping, JRK07 used a 
fixed rest-frame color for a given rest-frame filter; e.g., B — V color was used when /' = B. In this example, the new 
code uses cither B ~ V or B U for the warping, depending on whether the value of Ajost is closer to the central 
wavelength of F or t/. 

B. FILTER SET FOR THE NEARBY SN lA SAMPLE 

The Nearby SN la sample serves as the low-redshift anchor for the Hubble diagram. A superset of the nearby sample 
was also used to derive ("train") model parameters for the mlcs2k2 method. Within the framework of mlcs2k2, here 
we discuss the filter-response functions and color terms needed for the K-corrections that transform from rest-frame 
model magnitudes {U BVRI) to observer-frame magnitudes. 

The nearby SN la sample is a heterogeneous data set for which the filter response functions vary. While some of the 
filter curves have been published, the applicability of the filter curves to precision photometry is uncertain. However, 
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all of the nearby SN la magnitudes have been transformed into the Landolt system (Landolt & Umoto 2007), even if 
the eolor-terms arc not always available. Although the Landolt magnitude system is well defined, there are no standard 
filter responses for this system and hence one cannot compute K-corrections. In order to compute K-corrections, we 
determine color transformations from the Landolt system to the standard UBVRI filter response functions defined 
by Bessell (1990). While previous analyses with SNe la used Vega as the primary reference, we note that Vega has 
not been measured by Landolt. Instead of using Vega, we define the primary reference from the Landolt network, 
BD+17°4708 (hereafter BD+17), which happens to be the primary reference for SDSS photometry and also has a 
precisely measured HST spectrum. 

In order to determine the Bessell-Landolt transformation, we use Landolt standards that have excellent spectropho- 
tometric data and compare the observed Landolt magnitudes to synthetic Bessell magnitudes based upon the spec- 
trophotometric data and knowledge of the Bessell response functions. We use spectra from HST CALSPEC 2006'*^ 
(Bohlin 2007) because of its availability, high quality, and consistent calibration. For the Landolt standards that 
overlap with HST CALSPEC, the Landolt measurements are given in Table 18 and the synthetic magnitudes in the 
Bessell system in Table 19. 

To proceed, we define a synthetic Landolt-magnitude in passband "X" by 

^synth - ^synth + ^-^synth , 

where A^^^ft^ is the synthetic magnitude constructed from the source spectrum and the filter response from Bessell 
(1990), and AX 

synth is a correction based on color transformations as follows. 
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(B2) 

The subscript "obs" refers to observed (instrumental) magnitudes for a standard star or supernova, and the BD+17 
subscript indicates a Landolt measurement. These color transformations are defined so that there is no correction 
for the reference star, BD+17. The color coefficients (fc,-o.4) are determined by fitting Eq. B2 with the Landolt 
measurements and synthetic Bessell magnitudes in Tables 18-19, with Vega excluded from the fit. The resulting ki 
values are shown in Table 20; typical values are in the few percent range. We have also calculated the color terms 
determined with Vega as the primary reference, with results given in the last column of Table 20. The use of color 
terms is an approximation that can lead to significant errors, but the small size of the color terms indicate that this 
error should be negligible. 

As an alternative to using color transformations, we follow Astier et al. (2006) and define a modified set of UBVRI 
Bessell filter response functions in which the central wavelength of each filter passband is shifted from that of Bessell 
(1990) but the shape of the response curve is unchanged. The shifts arc defined such that the color terms, as defined 
in Eq. B2, are zero. The corresponding wavelength shifts relative to the filter responses defined by Bessell (1990) are 
given in Table 21 under the column "HST standards." The shifts used in Astier et al. (2006) are also given in Table 21 
for comparison. The differences in wavelength shifts are likely due to the different choices of spectral standards: we 
use HST standards (Table 18), while Astier et al. (2006) used ground-based spectra from Hamuy et al. (1992, 1994). 

The recipe for K-corrections is as follows. The UBVRI model magnitudes for mlcs2k2 are assumed to be in the 
Landolt system. The inverse of Eq. B2 is used to convert these Landolt magnitudes into magnitudes in the Bessell 
(1990) system, and then a K-correction is applied in the usual manner. For the SDSS, ESSENCE, and SNLS surveys, 
the filter response functions are well understood, and therefore the K-corrections are well-defined. For the nearby 
SN sample, the K-correction transforms into a hypothetical observer-frame UBVRI system with a filter response 
described by Bessell (1990); in this case, Eq. B2 is applied again to transform back to Landolt magnitudes. 

C. PRIMARY MAGNITUDES FOR BD+17 AND VEGA 

For the mlcs2k2 method, the primary magnitudes for each filter system are given in Table 22 for BD+17 and Vega. 
The magnitudes for UBVRI are taken from Landolt measurements, and the SDSS gri magnitudes are given in the 
AB system. The primary magnitudes in the other filter systems are obtained by the interpolation method described 
below. 

We compute the BD+17 magnitudes by first interpolating UBVRI magnitudes from Vega and using these mag- 
nitudes to determine offsets that are used to correct the synthetic magnitudes (from HST spectra) for BD+17. The 
Vega interpolation is based solely on the central wavelength of each passband and docs not depend on the detailed 
shape of the transmission curve. This approach is reasonable since the Vega magnitudes and colors are all small, and 
the Vega spectrum is smooth. 

Since we use the published magnitudes for the SN data, using our own analysis of the BD+17 magnitudes means 



'-' http: //www. stsci . edu/hst/observatory/ cdbs/ calspec .html 
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that wc have effectively adjusted the photometry of the pubhshcd data. The imphcit assumption is that the original 
photometry is correct relative to Vega, but that Vega itself has a slightly different value than was assumed previously. 
The Vega and BD+17 magnitudes do not agree exactly with the difference expected for the synthetic magnitudes 
computed from the HST spectra, so Vega and BD+17 define slightly different photometric systems. While we believe 
that our use of BD+17 provides an accurate and more consistent description of all the photometric data, it should be 
emphasized that our assumed magnitudes for Vega differ by amounts that are small compared to the zeropoint errors 
quoted in the original publications and the differences are well within our quoted systematic errors. 
As an example, we illustrate the determination of the primary BD+17 magnitudes for the SNLS griz passbands. 

TABLE 18 

Measurements from Landolt 2007 that overlap with HST CALSPEC data. N^bs is the number of Landolt observations. 
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-1.236 


-0.149 


-0.201 


-0.350 






±0.0022 


±0.0016 


±0.0033 


±0.0022 


±0.0043 


±0.0049 


HZ44 <= 


40 


11.673 


-0.291 


-1.196 


-0.141 


-0.181 


-0.322 






±0.0016 


±0.0011 


±0.0027 


±0.0009 


±0.0011 


±0.0014 


HZ4 " 


51 


14.506 


0.086 


-0.675 


-0.074 


-0.060 


-0.136 






±0.0027 


±0.0017 


±0.0022 


±0.0013 


±0.0017 


±0.0020 



" Vega measurements are not from Landolt. The ^-magnitude is from the HST analysis of Bohlin (2006), and the colors are from Bessell 
et al. (1998). Vega error estimates are not given in the references. 

^ The reported V magnitude was adjusted by Bohlin (2000). The magnitude reported by Landolt (2007) is 11.781. The errors and the 
colors are from Landolt loc cit. 

"^Landolt (2006) 

Private communication from Landolt (2006) reports more observations but the same values as published previously (Landolt 1992). 
^Landolt (2007). 



TABLE 19 

For Landolt stars in Table 18, synthetic magnitudes are computed using HST CALSPEC spectra and filters defined by 

Bessell (1990). 



star 




V 


B 




U 


- B 


V 


~ R 


R-I 


BD+17°4708 


9, 


,4640 


0, 


,4430 


-0, 


.1830 


0. 


2980 


0.3200 


Vega 


0, 


,0169 


0, 


,0048 


0, 


.0213 


-0. 


0114 


-0.0086 


G191B2B 


11, 


,7777 


-0, 


,3021 


-1, 


.2475 


-0. 


1535 


-0.2016 


GD71 


13, 


,0371 


-0, 


,2265 


-1, 


.1476 


-0. 


1438 


-0.1840 


GD153 


13, 


,3509 


-0, 


,2583 


-1, 


.1933 


-0. 


1477 


-0.1916 


AGK+81°4211 


11, 


,9226 


-0, 


,3221 


-1, 


.2573 


-0. 


1583 


-0.2051 


BD+28°4211 


10, 


.5076 


-0, 


.3202 


-1, 


.2735 


-0. 


1581 


-0.2122 


BD+75°325 


9, 


.5301 


-0, 


.3136 


-1, 


.2672 


-0. 


1509 


0.1987 


Feige 110 


11, 


.8295 


-0, 


.2969 


-1, 


.2064 


-0. 


1475 


-0.1697 


Feige 34 


11, 


.1731 


-0, 


.3312 


-1, 


.2726 


-0. 


1334 


-0.1563 


GRW+70°5824 


12, 


,7515 


-0, 


,0605 


-0, 


.8913 


-0. 


1118 


-0.1287 


HZ21 


14, 


,6880 


-0, 


,3270 


-1, 


.2677 


-0. 


1346 


-0.1969 


HZ44 


11, 


,6606 


-0, 


,2643 


-1, 


.2213 


-0. 


1385 


-0.2001 


HZ4 


14, 


,4818 


0, 


1210 


-0, 


.6627 


-0. 


0879 


-0.0758 
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For the interpolation, we need Landolt UBVRI magnitudes referenced to BD+17, My^g^ = ^^Vogl'' + (^^bdi7 ~ 
M^^lj). For F-band, M^^'^;^{V) = 0.003 + (9.464 - 9.450) = 0.017 mag. The results for aU passbands are 
My^g^iUBVRT) =0.026, 0.022, 0.017, 0.028, 0.037 mag. We interpolate these UBVRI magnitudes to the central 
wavelengths of the SNLS griz filters to get interpolated Vega magnitudes iMy^°^^), with the results shown in the 

first row of Table 23. The zeropoint offests for synthetic magnitudes (ZP''-^"'^) are the differences between the in- 
terpolated and synthetic magnitudes for Vega. The BD+17 magnitudes in the Landolt system are then given by 
M^^'f^ = M^^l^ + ZP'y"*'', and the results are given in the seventh row of Table 22. 

D. DETERMINING THE UNDERLYING DISTRIBUTION OF EXTINCTION {Av) AND LIGHT-CURVE SHAPE (A) 

Within the framework of the mlcs2k2 model (§ 5.1), we extract the underlying Ay and A distributions from the 
SDSS dust sample (§ 7.1) by making use of simulated light curves processed in exactly the same way as the observed 
light curves. The underlying distributions of Ay and A are defined such that when these underlying distributions 
are input into the simulation, the fitted distributions from the simulated light curves match the data distributions. 
We assume that these underlying distributions are independent of redshift, and we fix the CCM89 extinction law 
parameter to Ry = 2.18, following the analysis of § 7.2. 

Here we use the Bayesian unfolding method of D'Agostini (1995), and tests with simulated samples shows that one 
iteration is adequate with our statistics. For this discussion, we use an asterisk superscript to indicate the underlying 
true value for a parameter; the lack of an asterisk indicates a measured value obtained from fitting with the mlcs2k2 
method. Simulated light curves are generated with a flat distribution in Ay (over the range to 4), a flat distribution 
in A* (from —0.6 to -1-2), and a flat distribution in redshift (from to 0.4). For each simulated light curve that passes 
the selection criteria of § 4, the mlcs2k2 model is used with a flat, non-informative prior on Ay and A to extract a 
fltted value of Ay and A. These fitted values, in general, are different from the underlying Ay and A*. The accuracy 
of the inferred extinction. Ay — Ay, has a typical root-mean-square (rms) of 0.2, with little dependence on redshift. 
The accuracy of the fitted shape-luminosity parameter, A — A*, has an rms of a few hundredths for z < 0.1, and the 
rms increases to about 0.2 at z ~ 0.3. 

We use the simulation to calculate the conditional probability of extracting fitted values of Ay and A from a light 
curve generated with values Ay and A*, as a function of redshift. Since the generated, fitted, and observed values are 
drawn from continuous distributions, we bin these quantities with bin-sizes 0.2 for Ay, 0.1 for A, and 0.05 for redshift. 
Each bin is identified by an index: x* = ,*ai** f^'^'' generated quantities, and x = Mv,*A,iz for fitted quantities. 
Since redshifts are spectroscopically determined with high precision, we always have iz = i%. From the simulation, we 
calculate the conditional probability distribution, 

PsiAS\x*) = n{x)/N{x*) , (Dl) 

where n{x) is the number of fits producing values of {Ay , A, z) that lie in the specified bins, and N{x*) is the number 
of light curves generated with {Ay, A*, z*) that lie in the specified bins. The underlying two-dimensional distribution 
of Ay and A* is then obtained from the SDSS-II data by 

P{A*y,A*)^Y. 

X 

TABLE 20 

Calculated color coefficients for BD-i-17 and for Vega. 









value for 


value for 


Band 


coefficient 


BD+17°4708 


Vega 


V 




fco 


-0.010 ± 0.004 


-0.015 ±0.004 


B - 


V 


ki 


+0.027 ± 0.005 


+0.033 ± 0.005 


U- 


B 


k2 


-0.035 ± 0.004 


-0.050 ± 0.004 


V - 


B. 


ka 


-0.010 ± 0.005 


+0.004 ± 0.005 


R- 


I 


ki 


-0.029 ± 0.007 


-0.010 ± 0.007 



esSr.{x*)Y.APsSr.{x\x*)] 



(D2) 



TABLE 21 

Wavelength shifts for the Bessell (1990) filters. 

filter shift in A for: 
Bessell HST Astier et al. 



filter 


standards 


(2006) 


U 


+13 ±4 




B 


-15±4 


-41 


V 


+12 ±6 


-27 


R 


+ 7 + 9 


-21 


I 


-45 ± 21 


-25 
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TABLE 22 

Primary Landolt magnitudes used for K-corrections. 



filter system 


primary 


magnitudes 








* Landolt U, B, V, R, I 


BD+17 


9.724, 9.907, 


9.464, 


9.166, 


8.846 




Vega 


0.017, 0.021, 


0.023, 


0.030, 


0.026 


SDSS 2.5m g,r,i 


BD+17 


9.644, 9.350, 


9.256 






(AB system) 


Vega 


-0.106, 0.142, 0.356 




1^ CTIO(ESSENCE) R, I 


BD+17 


9.152, 8.855 










Vega 


0.024, 0.020 








CFHT(SNLS) g,r,i,z 


BD+17 


9.720, 9.222, 


8.911, 


8.756 






Vega 


0.016, 0.022, 


0.021, 


0.016 




HST (FllOW, F160W, 


BD+17 


8.558, 8.141, 


9.337, 


8.898, 


8.746 


F606W, F775W, F850LP) 


Vega 


0.003, 0.000, 


0.020, 


0.021, 


0.015 



''Measured by Landolt. 

'^Interpolated as described in the text. 

where esim = fsubtr x Ccuts is the simulated efficiency that includes the combined effects of the image-subtraction pipeline 
(§ 6.2) and selection cuts (§ 4), N^ataix) is the number of observed SNe la in rcdshift bin with fitted Ay, A that 
lie in the specified bins, and the summation is over the three-dimensional grid of observed Ay, A, and rcdshift. The 
corresponding one-dimensional distributions, P{Ay) and P(A*), are obtained by marginalizing the two-dimensional 
distribution over the other variable. We have extensively tested this procedure for extracting the true distributions 
of Ay and A* on simulated data samples, and the technique gives excellent agreement between the true (generated) 
distributions and the distributions extracted from the simulated data. 

We apply this technique to the SDSS-II dust sample; the extracted Ay and A* distributions are shown in Figure 46. 
The error bars reflect the statistics of our data sample as well as the statistical uncertainties in the simulated efficiency 
and Psini distributions. We fit these distributions to analytic functions that can be easily computed for the mlcs2k2 
fitting prior. The underlying Ay distribution is accurately fit by an exponential function, 

dN/dA*y = exp(-At./Tv) . (D3) 

The A* distribution is described by a bifurcated Gaussian with peak position Aq and different positive side and 
negative side widths, (j+ and cr— respectively. We find 

TV - 0.334 ± 0.050,tat ± 0.072,y,t (D4) 

Ao = -0.24 ± 0.029stat ± 0.013,y,t (D5) 

a- = 0.24±0.046stat± 0.022,y,t (D6) 

a+ = 0.48 ± 0.029stat ± 0.015,yst , (D7) 

where the uncertainties are statistical and systematic respectively. These distributions are the basis for the mlcs2k2 
priors shown in Figure 17. 

In addition to the statistical errors, we have considered a number of systematic effects in the determination of these 
distributions. We vary Ry by ±1ct (see § 7.2); this variation has the largest effect on the inferred Ay distribution. 
We vary the maximum rcdshift for the SN la light curves included in the data sample from the nominal 0.30 to 0.25. 
We also consider the difference between the two different methods for modeling SN la intrinsic luminosity variations 
(see § 6.1). Numerous minor variations in the analysis procedure, such as changing bin sizes and the use of binned 
or smoothed efficiencies, have a negligible effect. The uncertainties are summarized in Table 24 and are summed in 
quadrature to obtain the total uncertainty. 

As noted above, an exponential shape for the distribution of Ay agrees very well with the data, and we adopt 
this shape in the mlcs2k2 prior. Nevertheless, it is useful to check how well other proposed distributions match the 
data. In particular. WV07 considered a "galaxy line of sight" (glos) distribution for Ay. The glos model includes an 
exponential Ay distribution plus a narrow Gaussian with zero mean and small width that is meant to represent SNe la 

TABLE 23 

Magnitudes used to compute Landolt BD+17 magnitudes for SNLS griz. 



type of 




Magnitude for: 




magnitude 


g 


r 


i 


z 


^intcrp 


0.020 


0.025 


0.035 


0.030 


^synth 
-'"Voga 


-0.099 


0.149 


0.376 


0.513 


2psyiith 


0.119 


-0.124 


-0.341 


-0.483 


, ,synth 


9.601 


9.346 


9.253 


9.240 
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TABLE 24 

Uncertainties in the parameters that describe the distributions of Ay (tv) and A* (Ao, 0-+, cr— ). 



source of 




Uncertainty on: 




uncertainty 


TV 


Ao 


0-+ 


fj— 


statistical 


0.050 


0.054 


0.051 


0.033 


Rv =2.18 ±0.50 


0.050 


0.005 


0.004 


0.008 


redshift range 


0.040 


0.012 


0.014 


0.004 


simulated color smearing 


0.030 


0.003 


0.001 


0.006 


analysis details 


0.012 


0.001 


0.001 


0.002 


Total Systematic 


0.072 


0.013 


0.015 


0.011 


Total Uncertainty 


0.088 


0.056 


0.053 


0.035 



that would be expected to be observed with neghgible host galaxy dust extinction on half of the observed lines of 
sight, given a random distribution of observed host-galaxy orientations. When we fit the inferred Ay distribution with 
the glos model, wc find poor agreement between the best-fit model and the data, as shown in Fig. 46. Even allowing 
the most extreme of the systematic variations (e.g., varying Ry), the best agreement between the data and the glos 
model gives a confidence level of only 2%, calculated from the between the data and the model. If we allow the 
relative amplitude of the narrow Gaussian component of the glos model to vary, the fit returns an amplitude consistent 
with zero. The uncertainty on ry, the slope of the exponential, therefore accurately describes the uncertainty in the 
functional shape of the Ay distribution. 

E. DISCUSSION OF HUBBLE SCATTER FOR THE SDSS SAMPLE 

As noted in §10.1, the Hubble scatter and xfi statistic are significantly smaller for the spectroscopically confirmed 
SDSS-II SN la sample than for the other SN samples (sec Table 11). Here we investigate this anomaly by analyzing 
the SDSS-II host-z sample, described in §7.1, comprising 110 photometrically identified SNe la with spectroscopically 
determined host-galaxy redshifts. Recall that we required z < 0.3 for the SN sample used to measure host-galaxy dust 
properties; for the discussion below, we do not impose a redshift cut, thereby adding 31 host-z SNe la with z > 0.3. 

For the host-z sample, the rest-frame magnitudes at peak brightness are nearly 0.2 mag fainter on average than 
those for the spectroscopically confirmed SNe at the same redshift. The mean inferred extinction. Ay, is nearly 0.2 
mag larger for the host-z sample as well, indicating that the spectroscopic sample is not complete for intrinsically 
underluminous or extinguished events, as already inferred in §6. Performing a cosmological fit to the host-z sample 
alone results in xfi/^doi = 107/102 and RMS^ ^ 0.22 mag (where, as before, we set cr™' = 0.16), both of which are 
significantly larger than the corresponding values for the spectroscopically confirmed SNe la, xfi/^doi = 55.3/102 and 
RMS^ = 0.15 mag. A cosmological fit to the confirmed plus host-z SDSS-II sample results in xfi/^doi — 166/205 
(~ 98% probability) and RMS^ ~ 0.19 mag, consistent with the fit-quality parameters for the other SN la samples 
in Table 11. We therefore conclude that the lower xfi and Hubble scatter for the spectroscopically confirmed SDSS-II 
sample are largely caused by spectroscopic selection effects. As a crosscheck on the simulation, wc have also analyzed 
simulated spectroscopic and host-z samples, and find the corresponding RMS^ values to be in excellent agreement 
with the data. 

For completeness, we consider and exclude several other possibilities for the smaller scatter in the SDSS-II sample. 




Fig. 46. — The distributions of Ay and A* based on the SDSS-II dust sample and the procedure described in the text. The Ay 
distribution includes overlays of the best-fit exponential function (solid curve) and the glos prior (dashed). 
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1. One possibility is that the Scene Model Photometry method (SMP, §3.2) overestimates the flux uncertainties. 
Another is that SMP somehow provides a dramatic improvement in accuracy compared to the photometry 
methods used in other surveys. To test these possibilities, we have processed the SDSS-II spectroscopic sample 
photometry with the image-subtraction pipeline that was used for preliminary photometric measurements during 
the SN survey. The resulting cosmological fit yields very similar xfi and Hubble scatter, indicating that SMP is 
not the cause. 

2. To test if the mlcs2k2 fitter overestimates distance-modulus errors (cr^' in Eq. 27) for the SDSS-II sample, 
we have compared the average cr^' values for each of the SN samples at the mean redshifts of the samples. For 
the Nearby, SDSS-II, ESSENCE, and SNLS samples, the mean redshifts are 0.035, 0.22, 0.42, and 0.63. For SN 
sub-samples within small redshift windows centered on the mean redshifts, the average trj]* values are 0.07, 0.11, 
0.19, and 0.18 mag, with an uncertainty of ^ 0.01 on the average. For the SDSS-II, the corresponding average 
cr^' — 0.11, smaller than for the ESSENCE & SNLS surveys. Therefore, it appears that an overestimate of the 

distance-modulus uncertainty is not the cause of the smaller xfi for the SDSS-II sample. 

3. We have split the SDSS-II spectroscopic sample into lower- and higher-redshift halves at the median redshift 
z — 0.22. Both halves have consistent values of xfi- 

We conclude that the small xfi and Hubble scatter for the spectroscopically confirmed SDSS-II sample are primarily 
effects of the survey selection function, particularly the spectroscopic follow-up. To reduce biases in the cosmological 
analysis, it is important to either include the host-z sample in the analysis or to model the selection effects for the 
spectroscopic sample. In this paper, we have followed the second course. 

F. TOTAL-UNCERTAINTY CONTOURS IN THE W-Qm PLANE 

For the combined SN+BAO+CMB cosmology results, we describe a simple method to generate total-uncertainty 
contours, i.e, contours that include statistical and systematic errors, in the plane of w versus flu for the FwCDM 
model, and in the plane of versus 57de for the ACDM model. Figures 28-29 show the statistical and total error 
contours for the FwCDM and ACDM models using the mlcs2k2 method, and Figures 37-38 show the analogous 
contours using the SALT-ii method. 

A first-principles treatment of systematic errors would include all the systematic-error parameters and variations 
as nuisance parameters, evaluating the likelihood function on a multi-dimensional grid and then marginalizing over 
the nuisance parameters to obtain the likelihood for the cosmological parameters. For the large number of systematic 
effects we have considered, this approach would be computationally expensive. Instead, we take advantage of the 
empirical fact that the best-fit cosmological parameter results from the numerous systematic tests described in § 9 lie 
very close to a straight line defined by the BAO-hCMB prior, with slope dw/dflu — 5 for the FwCDM model, and 
diljyE/dilM — —0.8 for the ACDM model. For FwCDM systematic tests in which the w-variations are within 0.1, 
the root-mean-square w-scatter about this line is ~ 0.002. For larger w-variations, the curvature of the BAO+CMB 
prior becomes more noticable; for the largest w-variation of —0.3, the value of w lies 0.06 away from the straight 
line approximation. For the purposes of illustrating contours, this linear approximation is adequate. To incorporate 
the systematic errors, we stretch the statistical uncertainty contour along the line defined by the BAO+CMB prior, 
where the stretch factor is given by the ratio of total-to-statistical uncertainties on the cosmological parameter w: 
aw (tot) / (Ttu (stat) . 

Note that this approach is not valid in general: it depends on the relative precision of the SN results and the 
BAO+CMB constraints. For example, using a simulated SN sample with three times the data statistics of sample e, 
the scatter about the BAO+CMB line increases by a factor of several. 

G. MODIFICATION OF mlcs2k2 LIGHT-CURVE TEMPLATES TO MATCH salt-ii 

For the comparison in §11, we modified the mlcs2k2 model vectors so that the light-curve templates match synthetic 
light curves derived from the SALT-ii spectral surfaces. This translation of the SALT-ii model begins with the spectral 
surface as a function of rest-frame epoch and xi, with c = (Eq. 6). Synthetic SALT-ii UBVRI magnitudes are then 
calculated on a grid of epochs and xi values, and a ■ xi is added to each synthetic magnitude (with a = 0.12). The 
xi-grid is now relabeled with the mlcs2k2 parameter A using the relation obtained from the data samples, 

A = -0.1799 - 0.19022^1 + 0.038447a;? - 0.0043656a;? . (Gl) 

Next, an overall magnitude adjustment is made so that the peak F-band magnitude for A = matches that of 
the nominal MLCS2k2 model. For the final step, a quadratic fit of magnitude versus A is done for each epoch and 
each UBVRI filter. The resulting quadratic parameters M'''-^ , ■p'^'^ , and (f^^ , where e is an epoch index and /' 
is a filter index, define the SALT~-ii-modified mlcs2k2 model. The A-vs-magnitude fits were done in the interval 
—0.55 < A < +1.1, and the rms-scatter varies between 0.03 and 0.08 mag. 

The UBVRI model-magnitude errors are estimated from Fig. 6 of Guy et al. (2007). Compared to the mlcs2k2 
model errors, the SALT-ii errors are smaller near peak brightness and larger at later epochs. 
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